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3-1.  The  RFA  for  (M,N)  as  (8,7)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  amplitude  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  amplitude  data  (solid  line).  The 
corresponding  phase  plot  Is  given  In  Figure  3-2  .  22 


3-2.  The  RFA  for  (M,N)  as  (8,7)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  phase  plot  of  the  RFA  (dashed  line)  gives  a  good  fit 
to  the  original  phase  data  (solid  line).  The  corresponding 
amplitude  plot  Is  given  In  Figure  3-1 . 23 

3.3.  The  RFA  for  (M,N)  as  (20,19)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  amplitude  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  amplitude  data  (solid  line).  The 
corresponding  phase  plot  Is  given  In  Figure  3-2  .  24 

3-4.  The  RFA  for  (M,N)  as  (20,19)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  phase  plot  of  the  RFA  (dashed  line)  gives  a 
good  fit  to  the  original  amplitude  data  (solid  line). 

The  corresponding  amplitude  plot  Is  given  In  Figure 

3-2 . 25 

3-5.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  simulated 

data  using  the  first  six  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA.  The 
amplitude  plot  of  the  RFA  gives  a  good  fit  to 
the  original  amplitude  data  (solid  line).  The 
corresponding  phase  plot  Is  given  In  Figure  3-6  .  33 


3-6.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  simulated 

data  using  the  first  six  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  phase  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  data  (solid  line).  The  corresponding 
amplitude  plot  Is  given  In  Figure  3-5 . 34 

3-7.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  calculated  data 
for  the  sphere.  The  x's  are  the  data  points  used  for  the 
RFA.  The  amplitude  plot  of  the  RFA  (dashed  line)  gives  a 
good  fit  to  the  original  amplitude  data  (solid  line). 

The  corresponding  phase  plot  Is  given  in  Figure  3-8  ...  37 

3-8.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  calculated 

data  for  the  sphere.  The  x's  are  the  data  points  used  for 
the  RFA.  The  phase  plot  of  the  RFA  (dashed  line)  gives  a 
good  fit  to  the  original  amplitude  data  (solid  line). 

The  correspond!' ng  amplitude  plot  Is  given  In 

Figure  3-7 . 38 

3-9.  The  RFA  (M,N)  as  (16,15)  fit  to  the  calculated  data  for 
the  sphere.  The  x's  are  the  data  points  used  for  the 
RFA.  The  amplitude  plot  of  the  RFA  (dashed  line)  gives 
a  good  fit  to  the  original  amplitude  data  (solid  line). 

The  corresponding  phase  data  Is  given  In  Figure  3-10.  .  .  39 

3-10.  The  RFA  for  (M,N)  as  (16,15)  fit  to  the  calculated  data 
for  the  sphere.  The  x's  are  the  data  points  used  for  the 
RFA.  The  phase  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  phase  data  (solid  line) .  40 
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within  the  frequency  of  Interest  (window  II).  The  x's  are 
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3-13.  The  RFA  fit  to  the  calculated  data  for  the  sphere 

within  the  frequency  of  Interest  (window  III).  The  x's  are 
the  data  points  used  for  the  RFA.  The  amplitude  plot  of 
the  RFA  (dashed  line)  gives  a  good  fit  to  the  original 
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3-14.  The  RFA  applied  to  a  set  of  noisy  data,  where  the  Gaussian 
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3-19.  The  sue  operator  and  the  RFA  applied  simultaneously  to 
the  noisy  data  (solid  line)  for  the  sphere,  where  two 
Gausslanly  distributed  random  noises  have  been 
added  to  the  calculated  data  for  the  sphere.  The  x's  are 
the  data  points  used  In  the  RFA.  The  si gnal -to-nol se 
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4-1.  The  RFA  applied  to  the  calculated  data  for  the  sphere  In 
the  sense  of  least  squared  error  without  weight.  Note 
twenty  complex  data  points  are  used  (amplitude  plot).  .  .  85 
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4- 2.  The  RFA  applied  to  the  calculated  data  for  the  sphere  In 
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6-20.  Backscattered  Impulse  response  waveform  of  commercla 
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horizontal  polarization  . 
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CHAPTER  I 


INTRODUCTION 


In  future  generations  of  radar  systems,  the  natural  electromagnetic 
oscillations  of  a  radar  target  may  be  used  to  determine  the  physical 
properties  (size,  shape,  and  composition)  of  the  target.  The  key  to 
such  a  possible  advance  In  radar  systems  Is  broadband  signaling 
waveforms.  Signaling  waveforms  with  low-frequency  content  give 
Information  about  the  gross  size  and  shape  of  a  target,  and  waveforms 
with  high-frequency  content  give  a  detailed  description  of  Isolated 
parts  of  the  target.  The  Interaction  of  electromagnetic  waves  and 
material  objects  Is  best  summarized  by  the  canonical  response  waveforms 
of  the  object,  l.e.,  the  response  waveforms  for  Interrogating  signals 
with  Impulse,  step,  and  ramp  time  dependence. 

Such  a  linear  system  analysis  was  first  Introduced  for  finite 
objects  by  Kennaugh  and  Cosgrlff  [1]  and  was  later  formalized  by 
Kennaugh  and  Moffatt  [2].  The  concept  of  approximating  the  low- 
frequency  scattering  characteristics  of  an  object  by  rational  functions, 
l.e.,  complex  natural  resonances,  was  first  suggested  by  Kennaugh  and 
Moffatt  [2]  as 


"Keeping  In  mind  that  a  distributed-constant  representation 
Is  exact.  It  should  still  be  possible  to  represent  the 


low-frequency  properties  of  s-2Fr(s)  and  the  corresponding  time 
response  [which  would  be  Fo(t)3  by  that  of  a  lumped  constant 
network.  In  this  case,  FpCt)  would  be  approximated  by  a  series 
of  exponential  terms,  whereas  s"2Fr(s)  would  be  the  ratio  of 
polynomials  rather  than  the  series  (8)*.M 

A  rigorous  solution  to  the  complex  natural  resonance  li'^a  In 
analysis  known  as  the  singularity  expansion  method  (SEM)  has  been 
suggested  by  Baum  [3].  There  Is  much  that  Is  correct  about  this  complex 
natural  resonance  approach  [4,5];  however,  Kennaugh  [6]  has  suggested 
that  the  one  thing  that  Is  missing  from  the  SEM  Is  the  K-pulse  --a 
time-limited  Input  waveform  which  produces  special  time-limited  output 
waveforms.  While  the  main  purpose  of  this  dissertation  Is  not  to  settle 
this  question,  a  K-pulse  response  Is  approximated  for  simple  backscatter 
(a  circular  loop  for  edge-on  Incidence)  by  using  the  complex  natural 
resonances  extracted  In  this  dissertation. 

Accepting  the  fact  that  In  many  useful  situations  the  transient 
response  waveforms  of  a  scatterer  or  Its  corresponding  frequency 
response  can  be  well  approximated  by  lumped  parameter  models,  the 
purpose  of  this  dissertation  Is  to  present  various  methods  for  obtaining 
the  complex  natural  resonances  of  an  object  using  the  model.  Of 
particular  Interest  are  methods  useful  with  measured  data  which  are 
contaminated  by  noise. 

*  The  series  (8)  referred  to  above  was  a  Rayleigh- type  power  series 
In  s. 


The  complex  natural  resonances  of  an  object  can  be  obtained 
analytically  If  the  vector  wave  equation  Is  separable.  For  finite 
objects  therefore  there  are  two  geometries,  spherical,  which  was 
treated  by  Stratton  [7],  and  the  thin  circular  disc  [8],  which  has  not 
yet  been  studied  rigorously  from  a  complex  natural  resonance  viewpoint. 
The  complex  natural  resonances  of  all  other  finite  objects  must  be  found 
by  approximate  methods.  Perhaps  the  best  of  these  Is  a  quasi -rigorous 
approach  where  an  Integral  equation  for  the  object  as  a  scatterer  or 
radiator  Is  developed  and  then  converted  to  a  matrix  equation.  A 
numerical  search  of  the  complex  frequency  plane  Is  then  used  to  find  the 
complex  natural  resonances.  Mains  and  Moffatt  [9]  used  this  method  to 
find  the  complex  natural  frequencies  of  bent  wires  and  simple  wire 
airplane  models.  Unfortunately,  complex  geometries  such  as  aircraft  or 
ships  can  only  be  handled  by  this  approach  for  small  electrical  size. 

For  simple  geometries,  asymptotic  theories  such  as  the  Uniform 
Geometrical  Theory  of  Diffraction  (UTD)  can  be  used,  as  by  Kennaugh  [6], 
to  extract  all  but  the  lowest  frequency  resonances. 

For  complex  geometries  In  the  resonance  region  the  only  recourse  at 
this  time  Is  measured  experimental  data  at  real  frequencies.  Useful 
methods  for  extracting  complex  natural  resonances,  therefore,  must  be 
applicable  to  experimental  data.  Such  data  are  Inevitably  contaminated 
by  some  noise  and  clutter,  which  further  complicates  the  extraction 
problem.  The  experimental  data  utilized  In  this  dissertation  are 
frequency  domain  data.  To  test  time  or  other  domain  methods  the 
appropriate  data  are  generated  synthetically. 


Reviewed  briefly  In  Chapter  II  of  this  dissertation  are  those 
methods  for  extracting  the  complex  natural  resonances  of  a  scatterer  or 
a  radiator  which  have  been  suggested  and  shown  to  have  some  success  when 
applied  to  real  or  synthetic  data.  A  rational  function  approximation 
model  Is  developed  In  Chapter  III.  Application  of  this  model  to 
extraction  of  the  poles  of  sets  of  simulated  data  for  a  conducting 
sphere  are  presented.  A  window  technique  is  developed  for  the  rational 
fuctlon  approximation.  The  effects  of  uncorrelated  Gaussian  white  noise 
added  to  the  calculated  data  on  the  location  of  extracted  poles  are 
demonstrated.  A  sum  operator  and  a  zero-phase-shift  digital  filter  are 
used  to  combat  noise.  An  estimate  Is  made  of  the  required  signal -to- 
nolse  ratio  for  which  the  complex  natural  resonances  are  still 
extractable.  In  Chapter  IV,  an  elgenanalysls  model  Is  reviewed,  the 
emphasis  being  the  application  of  the  model  using  real  frequency  data. 
The  application  of  singular  value  decomposition  to  elgenanalysls  Is 
studied.  Presented  In  Chapter  V  are  comparisons  of  the  complex  natural 
resonances  extracted  by  the  models  derived  In  this  dissertation  with 
models  and  methods  suggested  by  others. 

A  generalized  method  for  extracting  a  set  of  complex  natural 
resonances  from  multiple  frequency  complex  scattering  data  Is  described 
and  Illustrated  In  Chapter  VI.  The  scattering  data  are  real  data  taken 
from  measurements  made  at  the  ElectroSclence  Laboratory  and  are  a  true 
test  of  the  methods  and  techniques  developed.  Stressed  are  the 
techniques  which  can  be  utilized  to  assure  that  the  correct  approximate 


order  for  the  system  has  been  obtained  as  well  as  actual  (as  opposed  to 
pattern-fitting)  complex  natural  resonances  have  been  extracted.  The 
desirability  of  obtaining  the  same  or  very  similar  results  using  two 
somewhat  different  approaches  Is  discussed.  It  Is  assumed  that  the 
available  experimental  data  are  contaminated  by  noise.  For  the 
noise-free  case  there  are  a  number  of  techniques  by  which  the  proper 
complex  natural  resonances  can  successfully  be  extracted.  An  estimate 
Is  made,  based  on  the  material  of  this  dissertation,  as  to  the  required 
slgnal-to-nolse  ratio  which  must  be  achieved  by  the  experimental  data 
before  extraction  of  complex  natural  resonances  Is  feasible. 

Conclusions  and  recommendations  based  on  the  research  reported  In 
this  dissertation  are  given  In  Chapter  VII.  New  analytical  results  for 
complex  natural  resonance  prediction  using  asymptotic  scattering 
estimates  and  first  estimates  of  certain  K-pulse  waveforms  and  response 
are  given  In  the  appendices. 


CHAPTER  II 


AN  OVERVIEW  OF  RESEARCH  ON  THE  EXTRACTION  OF 
COMPLEX  NATURAL  RESONANCES 


The  complex  natural  resonances  (poles)  correspond  to  the  exterior 
resonant  modes  of  an  object  [4].  The  exterior  resonant  frequencies  are 
complex,  where  the  real  parts  account  for  ray  divergence  and  bending 
losses,  and  the  Imaginary  parts  accrunt  for  resonances.  These  poles  are 
of  particular  Interest  In  the  waveform  feature  extraction  because  they 
are  related  to  the  physical  characteristics  of  an  object.  Furthermore, 
the  complex  natural  resonances  are  Independent  of  the  target  aspect  and 
polarization  excitation  [9].  These  poles  may  form  a  minimum  set  of 
parameters  needed  to  characterize  a  target.  Research  on  the  extraction 
of  the  complex  natural  resonances  has  been  ongoing  for  a  decade,  and 
many  engineering  approaches  to  the  complex  natural  resonances  have  been 
established.  This  chapter  summarizes  some  of  the  approximate  methods 
discussed  In  the  literature  which  have  been  used  to  extract  the  complex 
natural  resonances  of  a  radiator  or  scatterer  and  have  shown  some 
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A.  SIGNAL  FLOW  GRAPH 

Signal  flow  graph  methods  have  been  used  extensively  In  network 
analysis,  for  example.  Circuit  and  System  Theory  written  by  Lago  and 
Bennlngfleld.  However,  the  application  of  the  signal -flow-graph  concept 
to  electromagnetic  scatterers  was  first  suggested  by  Bagby  and  Kennaugh 
[10].  In  general  there  are  two  kinds  of  scatterers  being  studied:  thin 
wire  structures  such  as  straight  wires,  cross  wires,  wire  arcs,  circular 
loops,  etc.  and  solid  structures  such  as  spheres,  circular  cylinders, 
prolate  spheroids,  etc.  Descriptions  of  the  scattering  mechanisms  of 
these  two  basic  structures  follow  [10]. 

1.  Thin  Wire  Structures 

It  Is  assumed  that  there  are  four  operations  on  the  current  flowing 
on  the  wire  surface.  These  are:  a)  wire  path  distortion  and  phase 
delay,  b)  wire  end  and  wire  junction  reflections,  c)  wire  junction 
transmission  and  d)  wire  end  radiation  and  radiation  coupling.  The 
parameters  associated  with  these  operations  are  used  In  the  signal -flow 
graphs  to  predict  the  complex  natural  resonances.  For  a  straight  wire, 
the  accuracy  of  predicted  complex  natural  resonances  Is  very  good  for 
very  thin  wire  structure  [10],  l.e.,  wires  whose  length-to-dlameter 
ratios  are  greater  than  1000.  However,  the  accuracy  Is  decreased  If  the 
thickness  of  the  wire  Is  Increased.  Obviously,  the  formula  used  for  the 
current  distribution  should  be  further  modified  for  thick  wires.  It  Is 
very  hard  to  predict  numerically  the  four  scattering  mechanisms  for 
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different  wire  configurations.  The  application  of  the  signal  flow-graph 
to  wire  arcs t  crossed  wires  and  other  complex  wire  shapes  Is  limited  by 
the  accurate  values  of  the  scattering  mechanisms  mentioned  above.  The 
complex  natural  resonances  extracted  using  the  signal -flow-graph  method 
may  be  lacking  In  accuracy  for  the  reason  mentioned  above. 

Nevertheless,  the  signal  graph  method  yields  significant  physical 
Insight  Into  the  complex  natural  resonances  of  an  object. 

2.  Solid  Structures 

In  recent  years,  the  Uniform  Geometrical  Theory  of  Diffraction 
(UTD)  has  become  a  very  powerful  tool  for  solving  such  problems  as 
diffraction  from  a  straight  or  curved  wedge,  vertex,  and/or  curved 
surface.  The  UTD  edge  diffraction  field  and  curved  surface  waves 
derived  by  Pathak  and  Kouyoumjlan  [11,121  and  others  have  been  used  to 
evaluate  the  complex  natural  resonances  of  solid  structures  such  as 
spheres,  prolate  spheroids,  circular  discs,  etc.,  for  both  transverse 
electric  (TE)  and  transverse  magnetic  (TM)  modes  [10].  The  results  are 
fairly  close  to  those  obtained  from  the  exact  solution  (sphere)  and  from 
other  approximate  methods.  It  Is  noted  that  UTD  Is  a  high  frequency 
asymptotic  approximation  and  the  errors  In  the  poles  predicted  at  low 
frequencies  are  due  to  this  high  frequency  approximation. 

B.  POLES  EXTRACTION  FROM  THE  TRANSIENT  RESPONSE-PRONY'S  METHOD 

Prony's  method  was  originally  developed  to  solve  equations  which 
express  the  relationship  between  the  temperature  and  the  expansion  of 

8 


water  vapor  or  phenomena  defining  elastic  fluids  using  experimental 
data*  In  1795  [13].  However,  the  application  of  Prony's  method  to  the 
transient  response  of  electromagnetic  radiation  or  scattering  data**  to 
extract  poles  and  residues  was  first  suggested  by  Van  Blarlcum  In  1976 
[14].  In  the  first  application  of  Prony's  method,  numerically  generated 
data  for  the  transient  current  on  a  thin  dipole  were  used.  The 
resultant  poles  of  the  thin  dipole  were  very  close  to  those  extracted 
via  a  Integral  equation  formulation  and  numerical  search  by  Tesche  [15]. 
Inspired  by  this  Idea,  several  researchers  began  studying  Prony's  method 
for  other  scatterers  and  looking  also  for  solutions  to  the  problems 
Inherent  In  the  Prony  process.  E.M.  Kennaugh  suggested  the  application 
of  an  elgenanalysls  solution  (elgenanalysls  was  mentioned  In  [14]  but 
was  not  exploited)  to  replace  Prony's  method;  this  solution  was 
Illustrated  by  Moffatt,  Young,  Kslenskl  et  al.  [16].  Although  the 

*  The  equations  are  of  the  form; 

1 AX.  1 AX  1  AX  . 

Ti  «  Vl  +  v2p2  +  u3p3  +  .  .  .  ♦  Vn 

1  *  0,  1 . M-l ,  where  uj  and  pj  are  undetermined  constants,  ax 

Is  an  Increment  of  the  variable,  M  Is  the  total  number  of 
observation,  Is  an  observation  result  of  variable  at  Iax. 

**  The  equations  dealing  with  the  Impulse  response  of  electromagnetic 
radiation  or  scattering  data  are  of  the  form: 

p.lAt  p91 At  p,lAt  p  lAt 

Fi  *  Rie  1  +  R2e  2  +  R3e  3  R„e 

1  *  0,  1,  2,  ...,  M-l,  where  Aj,  pi  are  the  residues  and  poles,  M  Is 
the  total  number  of  the  sampled  data.  At  Is  the  time  Increment, 

Is  the  Impulse  response  at  time  lAt. 


complex  natural  resonances  of  an  object  or  class  of  objects  can  be 
determined  using  elgenanalysls,  the  accuracy  Is  seriously  affected  by 
even  moderate  amounts  of  noise.  Prony's  method  and  many  Prony's  related 
methods  are  summarized  In  a  recent  publication  [17].  Although  Prony's 
method  Is  a  useful  tool  for  extracting  the  complex  natural  resonances 
from  transient  data,  two  major  problems  which  occur  In  the  application 
of  Prony's  method  are:  1)  the  determination  of  the  system  order  and  2) 
the  bias  problem  with  noisy  data.  Many  researchers  have  focused  on 
these  two  problems  [17,18,19]  In  the  past  few  years,  but  there  Is  not 
yet  one  method  which  Is  completely  successful. 

C.  INTEGRAL  EQUATION  AND  NUMERICAL  SEARCH 

The  singularity  expansion  method  [SEM]  was  first  Introduced  by  Baum 
In  1971  as  a  technique  for  solving  transient  electromagnetic  scattering 
from  conducting  bodies  [3].  Although,  It  requires  a  search  procedure  In 
conjunction  with  a  space-frequency  Integral  equation  formulation,  a 
recent  publication  showed  that  the  Idea  of  the  singularity  expansion 
method  can  also  be  applied  to  a  space-time  Integral  equation  formulation 
[20].  In  the  application  of  Integral  equation  and  numerical  search  to 
determine  the  complex  natural  resonances,  an  Integral  equation  for  the 
Induced  current  Is  derived  and  reduced  to  a  finite  matrix  equation  using 
the  method  of  moments.  The  complex  natural  resonances  are  found  using 
an  Iterative  search  in  the  complex  wavenumber  plane.  The  complex 
natural  resonances  of  some  simple  scatterers  were  studied  using  this 
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method.  A  thin  wire  structure  was  studied  by  Tesche  In  1973  [15];  the 
prolate  spheroid  was  Investigated  by  Marin  In  1974  [21],  and  the 
circular  loop  was  done  by  Blackburn  and  Wilton  In  1978  [22].  The 
complex  natural  resonances  of  some  stick  models  of  the  M1gl9  and  F104 
were  studied  by  Mains  and  Moffatt  In  1974  [9].  At  the  Ohio  State 
University,  Professor  J.  H.  Richmond's  programs  for  backscatterlng 
calculations  for  the  aforementioned  airplane  models  were  used.  The 
natural  resonances  are  those  complex  frequencies  which  make  the 
admittances  In  a  moment  method  formulation  very  large  In  magnitude.  All 
of  these  studies  showed  some  success  in  the  extraction  of  complex 
natural  resonances.  However,  the  search  procedure  used  has  the 
following  weaknesses. 

1.  The  search  procedure  cannot  be  used  for  extraction  of  poles  of 
complex  structures  for  more  than  moderate  electrical  sizes. 

2.  The  search  procedure  Is  time  consuming  In  machine  computing. 

3.  The  search  procedure  cannot  be  used  to  process  measured 
scattering  data. 

0.  POLES  EXTRACTED  VIA  THE  RATIONAL  FUNCTION  APPROXIMANTS 

As  mentioned  in  Chapter  I,  the  electromagnetic  scattering  system 
can,  within  some  limitations,  be  modeled  as  lumped  parameter  system. 

The  canonical  responses  of  a  scatterer  are  the  Impulse  response,  step 
response,  and  ramp  response.  In  general  the  response  of  a  distributed 
parameter  system  can  first  be  separated  Into  a  forced  response  as  an 
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aperiodic  excitation  moves  over  the  object  and  later  a  natural  response 
as  the  excitation  moves  beyond  the  body.  In  this  dissertation,  the 
transfer  function  of  the  system  Is  approximated  by  a  rational  function 

E(Xi)  -  1  an  Xin+2/(l  +  l  bmX1 
n»0  m=l 

where  Xj  1$  jk^a  (an  electrical  length),  a  Is  a  linear  dimension  of  the 
object,  and  I  Is  the  number  of  the  sampled  data  points.  Thus  the  phasor 
response  of  a  system  Is  approximated  by  the  the  rational  function  E(Xf). 

The  first  attempt  to  extract  the  complex  natural  resonances  of  a 
scatterer  from  multiple  frequency  scattering  data  was  by  Moffatt  [23] 
using  a  rational  function  approximation  to  model  the  transfer  function 
of  the  target.  This  approach  was  largely  abandoned  with  the  advent  of 
Prony's  method  and  the  Improvement  of  Prony's  method— el genanaly sis.  In 
radar  applications,  with  the  exception  of  subsurface  radar,  the  data  are 
recorded  using  multiple  frequency  phasor  response  data  [16];  time 
response  waveforms  have  to  be  produced  synthetically. 

The  relative  advantages  of  pole  extraction  In  the  time  and  frequency 
domains  are  not  completely  obvious.  Assuming  that  the  Initial  data  are 
measured  frequency  domain  samples,  avoiding  approximate  transforms  Into 
the  time  domain  seems  desirable.  This  Is  clearly  true  If  the  data  are 
severely  bandllmlted  and  weightings  must  be  used  to  avoid  Gibbs  and 
aliasing  difficulties.  Also,  windowing  methods,  l.e.,  examining  various 
portions  of  the  measured  spectrum  individually,  appear  somewhat  more 


direct  In  the  frequency  domain.  A  detailed  discussion  of  the  windowing 
technique  Is  given  In  Chapter  III.  The  techniques  which  seek  to 
sequentially  extract  the  most  dominant  poles  are  generally  not  useful 
because  pole-extraction  methods  are  approximate  and  the  errors  will  be 
accumulative. 

No  clearly  superior  method  for  complex  natural  resonances 
extraction  from  measured  data  has  emerged  in  any  domain.  For  this 
reason  alone,  a  technique  utilizing  rational  function  approxlmants  Is 
deserved  of  some  additional  study. 


CHAPTER  III 


RATIONAL  FUNCTION  APPROXIMATIONS 


A.  INTRODUCTION 

As  discussed  earlier.  It  Is  assumed  that  with  the  possible 
limitations  mentioned  the  frequency-dependent  electromagnetic 
scattering  by  an  object  can  be  modeled  by  a  passive  linear  two-port 
system  with  time-invariant  parameters.  The  (normalized)  transfer 
function,  which  corresponds  to  the  (normalized)  impulse  response  of  the 
system,  can  be  approximated  by  a  rational  function  in  the  frequency 
domain  [2],  With  such  a  model  the  numerator  polynomial  will  be  aspect 
and  polarization  dependent  but  the  denominator  will  be  excitation 
Invariant. 

The  rational  function  model  for  electromagnetic  transfer  functions 
Is  written  Initially  as 

N  M 

En(X. )  =  l  a'  X.n+2  /  l  b'  X™  1  -  1,  2 . I  ,  (3-1) 

1  n=0  n  1  m=0  m  1 

where  En(X^)  (In  a  phasor  form)  Is  a  normalized  electric  or  magnetic 
field  at  frequency  X^ ,  Xi  is  jkjl  (an  electrical  length),  1  Is  a  linear 
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dimension  of  the  object  and  an'  and  bma  are  unknown  coefficients*.  I  Is 
the  number  of  sample  data.  One  of  the  coefficients  In  Equation  (3-1) 

Is  completely  arbitrary  because  the  numerator  and  denominator  can  be 
scaled  by  a  common  factor  without  changing  the  rational  function  model. 
The  zero  order  denominator  coefficient  In  Equation  (3-1)  Is  chosen  to 
be  normalized  and  accordingly  the  rational  function  model  (model  I)  Is 
rewritten  as 


E"(v  - 1 S  j,  6»  *?)  •  ,3-2' 

where  Xj  Is  jkja  and  a  Is  a  linear  dimension  of  the  object. 

It  Is  assumed  that  the  normalized  scattered  field,  En,  is  known  at 
a  total  of  I  real  ( jk-|1 )  frequencies.  The  proper  orders  of  the 
denominator  and  numerator  polynomials  are  unknown,  as  are  the  proper 
frequencies  ( jki 1 )  to  use  where  the  scattered  field  Is  assumed  to  be 
known.  Note  that  Equation  (3-2)  Is  selected  to  yield  a  Rayleigh-type 
dependence  if  the  scattering  data  samples  Include  the  Rayleigh  range. 
If  both  sides  of  Equation  (3-2)  are  multiplied  by  the  denominator 
polynomial,  then  I  complex  or  21  real  linear  simultaneous  equations 
can  be  generated.  Four  cases  must  be  considered  when  both  N  and  M  are 
allowed  to  be  even  or  odd.  There  are  N  +  M  +  1  unknown  coefficients 


*  For  a  real  physical  system,  the  Impulse  response  (In  the  time  domain) 
Is  causal  and  real.  Then  the  transfer  function  (i  the  frequency 
domain)  has  the  property  that  En(X  )  =  En  (X),  where  *  denotes 
complex  conjugate.  With  this  restriction,  the  coefficients  a',  b'  i 
Equation  (3_u  are  real.  "  m 


(an,  bnO  t0  be  found,  and  neglecting  the  underdetermi ned  case  these 
coefficients  can  be  found  In  an  exact  or  in  a  least  squared  error  sense 
depending  upon  the  selection  of  I.  One  would  anticipate  that  a  least 
squared  error  approximation  would  usually  be  better.  As  will  be  seen, 
however,  an  exact  solution  with  or  without  a  preprocessing  technique  has 
some  advantages  over  a  least-squared  error  approximation  under  certain 
conditions. 

To  this  point,  the  approximating  equations  given  in  Equation  (3-2) 
are  no  different  than  those  recently  given  by  Brittingham,  Miller  and 
Willows  [24].  In  fact,  they  are  the  same  as  those  originally  suggested 
by  Moffatt  [23].  The  idea  of  using  a  lumped  parameter  model  to 
approximate  a  distributed  parameter  system  goes  back  to  Kennaugh  and 
Moffatt  [2].  They  suggested  such  a  model  for  low  frequencies  In  1965 
(see  Chapter  I),  but  as  is  shown  in  this  dissertation  the  idea  can  be 
extended  to  rather  remarkable  lengths  In  terms  of  higher  frequencies. 

What  makes  the  the  present  approach  unique,  however.  Is  the  use  of 
non-harmonlcally  related  sampling  frequencies.  Also,  when  noisy  data 
are  considered,  the  methods  suggested  and  demonstrated  for  minimizing 
the  deleterious  effects  of  the  noise  are  felt,  to  some  extent,  to  be 
new.  The  final  step  in  the  application  of  Equation  (3-2)  is  routine 
once  the  coefficients  in  Equation  (3-2)  are  known.  The  zeros  of  the 
denominator  polynominal  yields  the  poles  of  the  system,  and  using  the 
numerator  coefficients.  Equation  (3-2)  can  be  rewritten  as  a 
residue  series,  i.e.. 


En(X)  =  l  Rj/(X-P.) 

j=l  1  1 


(3-3) 


V  .  vH+Z/  r, 


En(X)  =  l  a xn+z/  n  (X-P.) 
n=0  n  j=l  J 


(3-4) 


where  Rj  are  residues  and  Pj  are  poles.  Then, 


l  a  Xn+Z/  n  (X-P.)  I 
n*0  n  i=l  1  ' 


(3-5) 


assuming  all  poles  are  simple. 

Numerous  computer  programs  are  available  for  the  purpose  of  root 
extraction.  It  was  found,  however,  that  good  results  require  that 
double  precision  must  be  used  with  these  routines.  The  residues  and 
poles  are  in  complex  conjugate  pairs  for  the  transfer  function  of  a 
linear  causal  system  to  yield  real  time  domain  responses.  It  is  noted 
that  the  complex  natural  resonances  are  aspect  and  polarization 
Independent.  However,  the  residues  are  aspect  and  polarization 
dependent.  The  complex  data  EP(Xi)  can  be  expressed  in  a  phasor 


form  as 


En(Xi )  =  A{Xj)^  P(Xj) 


(3-6) 


where  A( X-f )  Is  the  amplitude  and  P(Xf)  Is  the  phase  at  frequency  of  Xi . 
To  measure  the  error  between  the  rational  function  approximation  and  the 
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original  data,  the  following  definitions  of  amplitude  errors  (norms)  are 
given: 


C-'t,!  lAljk,')  -  »RFA I2!1 


(3-7) 


EaS£  *  MAX  H^i1*  -  ARFA,(jkini  <sl>2 . 1}  ,  (3-8) 


where  A(jkfl)  and  arfa( 1 )  are  the  original  and  fitted  amplitude  at 

frequency  |^i,  and  E*^p  Is  the  root  mean  squared  error  (Euclidean  norm) 

max 

between  the  fitted  amplitude  and  the  original  one.  The  quantity  Eamp  Is 
defined  as  the  maximum  error  (Infinite  norm)  between  the  fitted  amplitude 
and  the  original  one.  Similar  definitions  are  used  for  the  phase  errors. 
The  dimension  of  the  phase  used  here  Is  degree. 

A  "best  approximate  solution"  of  a  system  generally  will  have 
different  solutions  for  different  choices  of  the  norm.  The  root  mean 
squared  error  and  the  maxlmun  error  are  calculated  In  each  run  using 
different  data  sets.  From  various  trials  and  errors,  the  least  error  In 
the  sense  of  Euclidean  norm  or  Infinite  norm  Is  selected  and  the  "best" 
fit  Is  obtained. 

In  applying  the  rational  function  approxlmants  to  a  set  of 
simulated  data,  there  are  a  number  of  parameters  which  must  be  chosen 
before  running  the  program,  the  orders  of  the  system  (M  and  N)  and  the 
selected  sample  data  points.  Therefore,  to  better  understand  the 
limitations  of  the  unknown  system  orders  and  to  gain  experience  In  the 


18 


n 


selection  of  the  data  points,  it  Is  helpful  to  consider  first  the 
simulated  real -frequency  backscatterl ng  data  of  a  conducting  sphere 
where  the  complex  natural  resonances  and  the  corresponding  residues  are 
well  known  [7,  25].  Some  of  the  lower  order  true  poles  and  residues  of 
a  conducting  sphere  (see  Table  3-1)  [25]  are  used  in  Equation  (3-9)  to 
generate  sets  of  simulated  data. 

ES(X)  =  l  tJ*1  ♦  )  ,  (3-9) 

1=1  *-Pl  X-Pf* 

where  ES(X)  are  the  simulated  data  and  R^,  P^,  are  poles  and 
residues  respectively.  X  Is  jkl  and  X*  denotes  complex  conjugate. 

In  generating  simulated  data  to  test  pole  extracting  procedures, 
actual  poles  and  actual  residues  of  the  sphere  were  used.  Much  too 
often  researchers  have  assumed  equal  excitation  of  all  poles  (a  very 
unrealistic  assumption)  and  correspondingly  shown  misleading  results  for 
particular  pole  extraction  techniques.  In  principle,  the  generated 
simulated  waveforms  could  be  exact  if  many  more  pole-pairs  were  used. 
However  even  here,  great  care  must  be  exercised  In  using  residue  series 
summations  [26]. 

B.  EXAMPLE  1:  SIMULATED  SPHERE  DATA  USING  FOUR  POLE-PAIRS 

In  this  example,  only  the  first  four  pairs  of  poles  and  the 
corresponding  residues  In  Table  3-1  are  used  in  Equation  (3-9)  to 
generate  a  set  of  simulated  scattering  data.  The  data  are  generated  In 


TABLE  3-1 


TRUE  POLES  AND  RESIDUES*  USED  FOR  GENERATING  SETS  OF 
SIMULATED  DATA 


Simulated  Poles 
(true  poles  of  a  sphere) 

-0.500000+/-j0. 866025 

-0. 701964+/- jl. 80740 

-0.842862+/-J2. 75786 

-0.954230+/-J3. 71478 

-1.04764  +/-J4. 57641 

-1.12891  +/-J5. 64163 


Simulated  Residues 
(true  residues  of  a  sphere) 

-0.0946447-/+0. 516674 

-0.633323  -/+0. 0853256 

-0.0802221+/-J0. 733736 

-0.822075  +/-J0. 0767481 

-0.0741270-/+J0. 901805 

-0. 0664705+/- jO. 223154 


*  The  units  of  the  poles  and  the  residues  are  In  lea,  where  k  Is 
the  wavenumber  and  a  Is  the  radius  of  the  sphere. 


the  real  frequency  domain  from  ka  of  0.2  to  ka  of  4  In  the  Increment  of 
0.02  ka,  l.e.,  191  data  points  are  generated.  Then,  a  rational  function 
approxlmant  of  order  (M,N)  as  (8,7)*  Is  applied  to  the  above  set  of 
data.  Since  there  Is  no  Rayleigh  region  In  this  set  of  simulated  data, 
the  dependence  Is  taken  out  from  Equation  (3-2).  Thus,  model  II  Is 
defined  as 

En(Xi)  *  l  anXj  /(l  ♦  J  bJC*  ]1  -  1,  2 . I  .  (3-10) 

1  n-0  n  1  m»l  m  1  / 

The  data  points  used  In  the  rational  function  approxlmant  are 
arbitrary  and  need  not  be  equally-spaced  In  frequency.  A  typical  plot  of 
a  rational  function  approximation  to  the  above  simulated  data  Is  In 
Figures  3-1  (amplitude)  and  3-2  (phase)  for  (M,  N)  as  (8,7)  and  In 
Figures  3-3  (amplitude)  and  3-4  (phase)  for  (M,N)  as  (20,19)**.  As  can 
be  seen  In  the  plots,  the  rational  function  model  with  the  above  system 
orders  fitted  all  the  data  points  so  well  that  the  original  and  the 
fitted  curves  of  both  amplitude  and  phase  cannot  be  distinguished  In  the 
plots. 

*  H  »  8  Is  the  smallest  system  order  for  the  simulated  data  using  four 
pole-pairs  (eight  poles);  N  of  7,  l.e.,  N  of  (M-l),  Is  the  best  choice 
for  the  RFA  to  fit  the  original  data  but  satisfies  Equation  (3-9), 
l.e.,  the  order  of  the  numerator  Is  one  order  less  than  that  of  the 
denominator. 

**  The  system  order  M  of  20  Is  chosen  because  the  RFA  programs  of  order 
higher  than  20  are  close  to  the  maximum  capacity  of  the  computer 
(VAX  11/780)  used  In  the  calculation. 


AMPLITUDE 


ml  part 


imaginary  part 


raal  part 

-i .e,s:c3  x:o-' 
-1.160983  X 1 0> 
-1.6-.5.93  Xio-' 
-3.808S3S  X 1 0-* 
-3.80353®  XIO-' 
1.398220  X 1 0° 
-1.42-1322  X10° 


ternary  part 

1 . S3t 33S  XIO" 
2. 635535  XIO'7 
-:.?S43&9  x : o8 
3.018S73  X10°, 
-3.018373  X10° 
0.300C00  X10° 

c.oocooo  x:o° 


-S. Si22£'4  X1C-’ 
-9.542384  XIO’1 
-8.425532  XIO'1 
-8.426552  XIO'1 
-7.019529  XIO'1 
-7.019525  XIO'1 
-S.0CC301  XIO'1 
-5.000002  X10'» 


-3.7:4782" X  T  0 3 
2.757552  >..:3 
-2.757335  XIO5 
1 . 6^7402  XIO** 
-1.807400  XIO5 
8.350247  x:0': 
-8.650243  XJO'1 


FREQUENCY  (  ka  ) 


Figure  3-1.  The  RFA  for  (M,N)  as  (8,7)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 
The  amplitude  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  amplitude  data  (solid  line).  The 
corresponding  phase  plot  is  given  In  Figure  3-2. 
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Figure  3-2.  The  RFA  for  (M,N)  as  (8,7)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  phase  plot  of  the  RFA  (dashed  line)  gives  a  good  fit 
to  the  original  phase  data  (solid  line).  The  corresponding 
amplitude  plot  Is  given  In  Figure  3-1. 
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Figure  3-3.  The  RFA  for  (M,N)  as  (20,19)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 
The  amplitude  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  amplitude  data  (solid  line).  The 
corresponding  phase  plot  Is  given  In  Figure  3-4. 
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Figure  3-4.  The  RFA  for  (M,N)  as  (20,19)  fit  to  the  simulated  data 
using  the  first  four  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 
The  phase  plot  of  the  RFA  (dashed  line)  olves  a 
good  fit  to  the  original  phase  data  (solid  line). 

The  corresponding  amplitude  plot  Is  given  In  Figure  3-3 


Besides  the  original  curve  (solid  line)  and  the  RFA  (dotted  line), 
a  third  one  (dashed  line)  is  generated  (see  Equation  (3-9))  using  the 
poles  and  residues  obtained  from  the  RFA.  Actually  these  procedures  are 
routine  In  the  RFA  program  to  check  the  accuracy  of  the  extracted  poles 
and  residues.  The  approximation  Is  such  that  only  the  solid  line  can  be 
seen  In  the  RFA  plots.  The  zeros,  poles,  residues  and  magnitude  of 
residues  are  listed  on  the  amplitude  and  phase  plots.  The  advantage  of 
these  listings  1$  to  see  the  true  poles,  curve  fitting  poles,  and  the 
cancellation  of  poles  and  zeros  before  the  true  poles  and  residues  are 
deduced.  It  Is  noted  that  only  M  complex  data  points  are  used  to  solve 
the  linear  system  equations  in  the  exact  sense.  The  average  errors*  of 
the  RFA  to  the  above  set  of  simulated  data  are  of  the  order  of  l.E-5  In 
amplitude  and  l.E-4  In  the  phase  (In  degrees).  The  numerical  accuracy 
of  the  entire  rational  function  program  was  tested  using  the  data  In 
Example  I  for  various  system  orders  and  various  Input  data  points.  In 
all  cases  the  programs  were  found  to  be  very  accurate.  Therefore,  it  Is 
assumed  that  the  program  can  be  applied  routinely  to  other  scattering 
data  without  worrying  about  numerical  errors  generated  within  the 
program  Itself. 

The  poles  and  the  residues  extracted  via  the  rational  function 
approxlmants  using  different  system  orders  M,  N,  and  different  sets  of 
selected  data  points  are  listed  In  Table  3-2  for  (M,N)  as  (8,7)  and 


*  Since  the  data  are  exact  for  these  poles  (residues),  only  trancatlon 
contributes  to  the  errors.  Obviously  the  error  Is  dependent  also  on 
the  selected  sample  points. 
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TABLE  3-2 


A  COMPARISON  OF  TRUE  POLES  AND  EXTRACTED  POLES,  TRUE  RESIDUES 
AND  EXTRACTED  RESIDUES  FOR  (M,N)  AS  (8,7);  FOUR  PAIRS  OF  POLES 
AND  RESIDUES*  ARE  USED  IN  THE  SIMULATED  DATA 


Poles 


True 

Extracted 

Percentage  of  Error 

Real  Imaginary 

-0.50000  +/-jO. 866025 

-0.50000  +/-J0. 866025 

0.E0 

O.EO 

-0. 701964+/- jl. 80740 

-0. 701 962+/- j 1.80740 

1.04E-4 

2.07E-4 

-0. 842862+/- j2. 75786 

-0. 842859+/- j2. 757866 

1.03E-4 

4.97E-5 

-0.954230+/-J3. 71478 

-0.954238+/-J3. 714782 

2.08e-4 

4.97E-5 

True 

Residues 

Extracted 

Percentage  of  Error 

Real  Imagl nary 

-0,0946447+/- jO. 516674 

-0.094644-/+J0. 516674 

1.33E-4 

O.EO 

-0.633323  -/+jO. 0853256 

0.633320-/+J0. 08533 

4.76E-4 

3.75E-4 

0.0802221+/-J0. 733736 

0.08024  +/-J0. 73373 

2.42E-3 

8.07E-4 

-0.822075  +/-J0. 07674781 

-0. 822088+/- jO. 07676 

1.57E-3 

1.48E-3 

*  The  units  of  the  poles  and  the  residues  are  In  ka,  where  k  Is  the 
wavenumber  and  a  is  the  radius  of  the  sphere. 


Table  3-3  for  (M,N)  as  (20,19).  The  maximum  and  the  average  discrepancy 
of  the  fitted  amplitude  and  the  original  amplitude,  the  fitted  phase  and 
the  original  phase  are  listed  in  Table  3-4.  Obviously,  the  smaller  the 
maximum  amplitude  discrepancy  or  the  smaller  the  average  amplitude 
discrepancy,  the  closer  the  fitted  curve  Is  to  the  original  one 
(slmilarily,  for  the  phase  error).  Correspondingly,  the  extracted  poles 
are  closer  to  the  original  poles.  The  percentage  error  (P.E.)  Is 
defined  for  both  real  part  and  imaginary  part  of  the  extracted  poles  as 

P.  E.  of  real  part  =  ^  P*rt  l^ext-PoUtine)!  ,  (3-11) 


P.  E.  of  Imaginary  part  -  part(Polee,t.-P°letrue)l  .  (3-12) 

"’Oletraol 

where  Poleext.  is  the  extracted  pole  and  Poletrue  the  true  pole. 
Similar  definitions  are  used  for  the  residues.  A  comparison  of 
extracted  poles,  the  percentage  error  of  extracted  poles  and  the  errors 
of  amplitude  and  phase  are  also  shown  in  Table  3-4. 

As  can  be  seen  in  Tables  3-2  to  3-4,  various  system  orders  of 
the  RFA  have  been  applied  to  the  same  data  set  to  see  the  effect  of  the 
system  order  on  the  accuracy  of  the  extracting  poles.  However,  no 
obvious  result  has  been  found  In  this  simple  simulated  example  except 
that  the  system  order  should  not  be  less  than  the  number  of  simulated 
system  poles.  For  real  targets  there  are  an  infinite  number  of 


TABLE  3-3 


A  COMPARISON  OF  TRUE  POLES  AND  EXTRACTED  POLES,  TRUE  RESIDUES 
AND  EXTRACTED  RESIDUES  FOR  (M.N)  AS  (20,19);  FOUR  PAIRS  OF  POLES 
AND  RESIDUES*  ARE  USED  IN  THE  CALCULATION  OF  THE  SIMULATED  DATA 


True  Poles 

Extracted  Poles 

Percentage  Error 

real  Imaginary 

-0.50000  +/-J0. 866025 
-0. 701964+/- j 1.80740 
-0. 842862+/- j2. 75786 
-0.954230+/-J3. 71478 

-0. 50000+/- jO. 86602 
-0.70196+/-jl. 80741 
-0.84283+/-J2. 75791 
-0.95421+/-J3. 71492 

0.E0  4.95E-4 

2.06E-4  5.16E-4 

1 .  HE— 3  1.74E-3 

5.22E-4  3.65E-3 

Curve  fitting  Poles* 

Curve-fitting  Zeros 

0. 19396+/- jO. 78060 
0.23181+/-J1. 51326 
-0. 39652+/- j 1.86546 
0.12284+/-j2. 97965 
0.18009+/-J3. 84572 

0 . 43016+/- j2 . 16078 

-0. 19396+/- jO. 78060 
-0.23181+/-J1. 51327 
-0. 39649+/- jl. 86548 
0.12284+/-J2. 97967 
0.18008+/-J3. 84573 

0. 43016+/- j2. 16079 

Extracted  Zeros 

-0.38086+/-J3. 01887 
-0. 16461+/- jl. 68440 
-1.16102+/-J0. 

1. 29822+/- jO. 

-1. 42429+/- jO. 

True  Residues 

Extracted  Residues** 

Percentage  Error 

-0.0946447-/+J0. 516674  -0.09463-/+J0. 51668 

0.633323  -/+J0. 0853256  0.63332-/+J0.0854 

0. 0802221+/- jO. 733736  0. 08046+/- jO. 7335 
-0.822075  +/-J0. 07674781  0.8219  +/-J0. 07706 

2.79e-3  1.15E3 

4.76E-4  1.16E-2 

3.22E-2  3.197E-2 

2.119E-2  3.781E-2 

*  The  units  of  the  poles  and  the  residues  are  In  ka,  where  k  Is  the 
wavenumber  and  a  Is  the  radius  of  the  sphere. 

**  All  extracted  residues  corresponding  to  the  curve-fitting  poles  are 
In  the  order  of  l.E-5  or  less  for  both  real  and  imaginary  parts. 
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TABLE  3-4 

1 

THE  EXTRACTED  POLES  AND  THE  ERRORS 

BETWEEN 

THE  RATIONAL  FUNCTION  FIT  AND  THE 

ORIGINAL 

- 

DATA  (SIMULATED  DATA  USED  IN  EXAMPLE  1)* 

■ 

True  poles 

RFA  Model  I 

M  =  14,  N  =  13 

Percentage 
Real  Part 

Error 

Imag.  Part 

: — - 

-0.500000+/-J 0.866025 

-0.50000  +/-J0. 866022 

O.EO 

2 . 98e-04 

- 

-0.70196  +/-J1. 80739 

-0.7019523+/-J 1.80743 

6.06E-4 

1.55E-3 

-0 . 842849+/- J  2.75786 

-0.84282  +/-J2. 75741 

1.457E-3 

1.56E-2 

| 

-0.954299+/-J3. 714787 

-0.9648  +/-J3.7253 

0.27 

0.27 

MAX.  AMP.  ERROR  0.19E-5 
AVE.  AMP.  ERROR  0.2E-6 
MAX.  PHA.  ERROR  0.1068E 

-3 

AVE.  PHA.  ERROR  0.155E- 

4 

1 

RFA  Model  II 

M  =  14,  N  *  13 

Percentage  Error 

Real  Part  Imag.  Part 

•  -  , 

-0.499998+/-J0. 86606 

2.E-4 

3.505E-3 

-0.70196  +/-J1. 80738 

2.06E-4 

1.03E-3 

m 

-0.84295  +/-J2. 75778 

3.05E-3 

2.77E-3 

m _ _ 

* 

-0.95515  +/-J3. 71550 

2.4E-02 

1.877E-2 

» *  •» 

[ 

MAX.  AMP.  ERROR  0.5E-6 
AVE.  AMP.  ERROR  0.1E-6 
MAX.  PHA.  ERROR  0.916E-4 

• 

b 

AVE.  PHA.  ERROR  0.148E-4 

• 

RFA  Model  I 

M  =  16,  N  =  15 

Percentage 
Real  Part 

Error 

Imag.  Part 

-0.499997+/-J0. 866023 

3.01E-4 

1.97E-4 

- 

-0.701981+/-jl. 80739 

8.76E-4 

5.16E-4 

- 

-0.84286  +/-J2. 75810 

7.03E-5 

8.33E-3 

-0.95048  +/-J3. 71478 

9.78E-2 

O.EO 

MAX.  AMP.  ERROR  0.19E-5 
AVE.  AMP.  ERROR  0.2E-6 

MAX.  PHA.  ERROR  0.176E-3 

L 

AVE.  PHA.  ERROR  0.187E-4 

— 

*  The  units  of  the  poles  and  the  residues  are  In  ka,  where 

a  Is  the 

radius  of  the  sphere  and  k  Is  the  wavenumber. 

. 
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singularities  In  the  system  transfer  function.  However,  over  a  limited 
spectral  span  only  a  finite  number,  hopefully  small,  will  contribute 
significantly. 

C.  EXAMPLE  2:  SIMULATED  SPHERE  DATA  USING  SIX  POLE-PAIRS 

It  Is  well  known  that  there  are  an  Infinite  number  of  pole  pairs 
for  a  scatterer  [4].  To  generate  a  set  of  simulated  data  similar  to  the 
real  scattering  ones,  six  pairs  of  simulated  poles  and  residues  (Table 
3-1)  are  used  in  Equation  (3-9).  The  simulated  data  (ka  of 
0.2(0. 02), 4.)  are  generated  in  the  same  way  as  that  In  example  1  except 
that  2  more  pole  pairs  and  the  correspondl ng  residues  outside  the 
original  spectrum  are  used  In  Equation  (3-9).  Once  again,  the  RFA 
routine  Is  used  to  extract  poles  from  this  set  of  simulated  data.  The 
original  data  and  fitted  data  are  plotted  In  Figures  3-5  and  3-6.  The 
extracted  poles  are  shown  on  Table  3-5. 

The  two  highest  frequency  pole-pairs  are  not  as  close  to  the 
original  ones  as  as  those  found  In  example  one  because  the  two  pairs  of 
poles  outside  the  original  spectrum  have  some  Influence  on  the  extracted 
poles,  especially  over  the  higher  frequency  portion  of  the  spectrum. 

This  Is  one  of  the  reasons  why  curve  fitting  poles  are  obtained  when  the 
rational  function  approxlmants  are  used  and  the  reason  why  the  system 
order  should  always  be  chosen  higher  than  the  number  of  poles 
actually  desired.  The  number  of  true  poles  Is  known  In  these  simulated 
data  but  the  number  of  true  poles  In  a  set  of  calculated  data  or 


TABLE  3-5 


A  COMPARISON  OF  TRUE  POLES  AND  EXTRACTED  POLES,  TRUE  RESIDUES 
AND  EXTRACTED  RESIDUES  FOR  (M,N)  AS  (14,13);  6  POLE-PAIRS  AND  THEIR 
CORRESPONDING  RESIDUES  ARE  USED  IN  THE  GENERATION  OF 
SIMULATED  DATA 


Poles  (M  =  14,  N  =  13)* 

True  Poles  Extracted  Poles  Percentage  Error 

Real  Imaginary 


-0.50000  +/-j0. 866025 

-0 .499995+/- jO . 866023 

5.E-4 

1.97E-4 

-0. 701964+/- jl. 80740 

-0.70203  +/-J1. 80746 

3.4E-3 

3.1E-3 

-0. 842862+/- j2. 75786 

-0.8427  +/-J2.7572 

5.62E-3 

2.29E-2 

-0. 954230+/- j3. 71478 

-0.9689  +/-J3.7213 

0.382 

0.17 

True  Residues 

Residues  (  M  *  14,  N  * 
Extracted  residues 

13) 

Percentage  Error 

Real  Imaginary 

-0.0946447-/+J0. 516674 
0.633323  -/+J0. 0853256 

0. 0802221+/- jO. 733736 
-0.822075  +/-J0. 07674781 

-0.094641-/+J0. 516647 
0.6337  -/+J0. 08573 
0.07635  +/-J0.7326 
0.924  +/-J0. 07676 

7.05E-4  5.14E-3 

5.8994E-4  6.328E-2 
0.5246  0.154 

12.34  6.63E-2 

*  The  units  of  the  poles  and  the  residues  are  In  ka,  where  k  Is 
the  wavenumber  and  a  Is  the  radius  of  the  sphere.  Also  the 
three  pairs  of  curve  fitting  poles  and  the  corresponding 
residues  are  not  shown  here. 
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Figure  3-5.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  simulated  data 
using  the  first  six  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  amplitude  plot  of  the  RFA  gives  a  good  fit  to  the 
original  data  (solid  line).  The  corresponding  phase  plot 
Is  given  In  Figure  3-6. 
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6.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  simulated 

data  using  the  first  six  pole-pairs  and  residues  for  the 
sphere.  The  x's  are  the  data  points  used  for  the  RFA. 

The  phase  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  data  (solid  line).  The  corresponding 
amplitude  plot  Is  given  In  Figure  3-5. 
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measured  data  Is  unknown.  A  detailed  discussion  on  the  determination 
of  the  system  order  Is  given  In  a  later  section. 

0.  EXAMPLE  3:  CALCULATED  BACKSCATTERING  DATA  FOR  A  CONDUCTING  SPHERE 

The  spectrum  of  a  conducting  sphere  [27]  used  to  test  rational 
function  approxlmants  Is  from  ka  of  0.2  to  ka  of  4  In  the  Increment  of 
0.02  ka(  l.e.,  191  data  points  are  available.  To  test  the  rational 
function  approximation  in  application  to  this  set  of  calculated  data,  M 
Is  set  to  be  the  number  of  poles  where  the  Imaginary  parts  are  within 
the  given  spectrum  first.  As  discussed  before,  the  extracted  poles  are 
not  accurate  If  the  system  order  M  Is  less  than  the  number  of  true 
poles.  However,  a  better  result  can  be  obtained  If  the  system  order  Is 
set  around  twice  the  number  of  the  true  poles  In  the  system.* 

After  the  polynomial  orders  (M  and  N)  are  selected  for  testing, 
different  sets  of  data  points  are  used  and  the  defined  errors  In 
amplitude  and  phase  are  calculated.  Table  3-6  Is  a  comparison  of  the 
extracted  poles  and  their  percent  errors.  Figures  3-7  to  3-10  are  plots 
of  the  rational  function  fit  under  two  different  sets  of  (M,N)  ((14,13) 
and  (16,15)).  Obviously,  the  original  curve  and  the  fitted  curve  agree 
so  well  that  they  cannot  be  distinguished  In  the  plots  of  both  amplitude 
and  phase.  However,  the  results  fit  better  for  (M,N)  as  (14,13).  If 
different  sets  of  selected  data  points  are  used  for  the  RFA,  the  results 

*  Basically,  the  extra  orders  (the  orders  which  are  higher  than  the 
number  of  true  poles)  are  required  for  the  curve  fittings 
although  the  curve-fitting  poles  have  no  physical  meaning. 


TABLE  3-6 


A  COMPARISON  OF  EXACT  POLES  AND  THOSE  EXTRACTED*  FROM 
CALCULATED  DATA  FOR  SPHERE  VIA  RATIONAL  FUNCTION 
APPROXIMATION 


True  Poles 

gxtrjjttd  Po]5s 

Percentage  Error 

Real  Part  Imag.  Part 

-0.50000  +/-J0. 866025 
-0.701964+/-J1. 80740 
-0.842862+/- j2. 75786 
-0.95423  +/-J3. 71478 

-0. 500000+/- jO. 86662 
-0.70193  +/-j 1.8028 
-0.84244  +/-J2. 76213 
-0.9175  +/-J3.6829 

O.EOO 

1.75E-3 

1.463E-2 

0.9576 

5.95E-2 

0.237 

0.148 

0.8312 

M  »  16,  N  «  15 
(exact  fit) 

-0.50087+/-j0. 86575 
-0. 71008+/- jl. 79815 
-0. 84535+/- j2. 76718 
-0.9076  +/-J3. 69581 

8.7E-2 

0.42E0 

8.63E-2 

1.22E0 

2.75E-2 

0.48E0 

0.32E0 

0.49E0 

*  The  units  of  the  poles  are  In  lea,  where  Ic  Is  the  wavenumber, 
a  Is  the  sphere  radius. 
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Figure  3-7.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  calculated  data 
for  the  sphere.  The  x's  are  the  data  points  used  for  the 
RFA.  The  amplitude  plot  of  the  RFA  (dashed  line)  gives 
a  good  fit  to  the  original  amplitude  data  (solid  line). 
The  corresponding  phase  plot  Is  given  In  Figure  3-8. 
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Figure  3-8.  The  RFA  for  (M,N)  as  (14,13)  fit  to  the  calculated  data 
for  the  sphere.  The  x's  are  the  data  points  used  for 
the  RFA.  The  phase  plot  of  the  RFA  (dashed  line)  gives  a 
good  fit  to  the  original  phase  data  (solid  line).  The 
correspondl ng  amplitude  plot  Is  given  In  Figure  3-7. 
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Figure  3-9.  The  RFA  (H,N)  as  (16,15)  fit  to  the  calculated  data  for 

the  sphere.  The  x's  are  the  data  points  used  for  the  RFA. 
The  amplitude  plot  of  the  RFA  (dashed  line)  gives  a  good 
fit  to  the  original  amplitude  data  (solid  line).  The 
corresponding  phase  data  Is  given  In  Figure  3-10. 
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may  differ  somewhat,  but  they  are  still  within  the  same  order  of 
accuracy. 

E.  IDENTIFICATION  OF  THE  TRUE  POLES  AND  THE  CURVE  FITTING  POLES 

With  proper  use  of  rational  function  approxlmants,  the  polynomial 
order  required  Is  always  higher  than  the  number  of  true  poles. 

Therefore,  the  following  rules  are  used  to  help  to  Identify  the  true 
extracted  poles  and  curve  fitting  poles  In  any  test. 

1.  Pole  locations  found  which  are  In  the  right  half  plane  are 
assumed  to  be  curve-fitting  poles.  Clearly,  right-half  plane  poles 
would  make  the  transient  response  unstable. 

2.  Poles  whose  corresponding  residues  are  relatively  small  compared 
to  others  obtained  with  that  test  are  assumed  to  be  curve-fitting  poles. 
The  reason  Is  that  the  pole-pairs  having  very  small  residues  contribute 
little  to  the  time  domain  response  or  to  the  frequency  domain  response. 
Models  with  higher  orders  have  too  much  detail  about  the  data.  The  fine 
detail  of  the  data  may  be  due  to  the  numerical  errors  from  the  machine 
computation  or  from  using  a  less  perfect  model  In  the  original  data 
calculation. 

3.  Poles  which  can  be  approximately  cancelled  by  zeroes  are  curve 
fitting  poles.  This  type  of  cancellation  Is  due  to  the  polynomial 
orders,  M  and  N,  being  overdeterml ned  (too  large)  both  In  the  denominator 
and  numerator. 


4.  Poles  which  have  only  a  real  part  or  whose  real  parts  are  very 
large  compared  with  other  extracted  poles  are  assumed  to  be  curve 
fitting  poles.  Poles  with  large  negative  real  parts  make  the 
corresponding  transient  response  decay  so  quickly  that  they  make  little 
contribution  In  the  total  transient  response. 

Therefore,  all  of  the  above  poles  are  assumed  to  be  due  to  noise  or 
errors  In  the  original  data  and  are  classified  as  curve  fitting  poles. 

It  Is  true  that  poles  with  only  real  parts  exist  for  certain  targets, 
the  conducting  sphere  1$  an  example.  A  real  pole  showing  consistently 
for  numerous  tests  would  be  reexamined. 

F.  THE  SELECTION  OF  DATA  POINTS 

If  the  matrix  equations  are  solved  In  the  exact  sense,  the  number 
of  data  points  used  ((M+N+D/2)  or  ((M+N)/2+l)  Is  limited  by  the  number 
of  unknowns  (H+N+l).  Equally  spaced  Interpolation  may  or  may  not  yield 
the  best  fit  for  a  rational  function  approximation.  On  the  other  hand, 
data  points  sampled  at  relative  maxima  and  minima  of  amplitude  are 
always  a  good  choice  for  the  best  fit.  This  Is  deduced  from 
observations  of  the  applications  of  the  rational  function  approxlmants 
to  numerous  sets  of  scattering  data.  Furthermore,  In  most  cases,  one 
additive  data  point  chosen  between  the  adjacent  relative  maximum  and 
minimum  Is  helpful  In  making  the  error  smaller. 

Although  there  Is  no  general  proof,  the  following  technique  Is  used 
and  works  reasonably  well.  The  amplitude  and  phase  error  at  all  the 


data  points  must  be  computed.  For  points  where  the  error  is  large,  the 
sample  points  are  moved  closer  together;  and  where  the  local  error  Is 
small,  the  samples  are  moved  farther  apart.  With  this  new  choice  of 
sample  points,  the  process  Is  repeated  until  the  error  Is  as  small  as 
possible  and  is  distributed  evenly. 

If  tiie  matrix  equation  Is  solved  in  a  least  squared  error  sense, 
there  Is  no  limitation  on  the  number  of  data  points  used  in  solving  an 
overdetermined  system.  (However,  the  sampling  rate  has  been  shown  to  be 
related  to  the  accuracy  of  the  extracted  poles  using  the  least  squared 
error  solutions  in  Prony's  Method  [28].)  The  application  of  least 
squared  error  solutions  to  pole  estimation  in  the  frequency  domain  using 
the  RFA  will  be  shown  in  Chapter  IV. 

G.  LIMITATION  OF  THE  SYSTEM  ORDER  M 

For  a  large  system  order  (for  example,  M  greater  than  20),  the 
rational  function  approxlmant  program  may  fall  due  to  overflows  in  the 
computer  during  the  computation.  The  machine  used  in  this  application, 
VAX  11/780,  has  a  capacity  In  the  range  of  0.29  E-38  to  1.7  E+38  for 
floating  point  numbers  in  double  precision.  In  our  applications,  the 
overflow  happens  in  the  subroutine  Dpolyrt  (a  zero-searching  routine) 
during  its  calculation  of  zeros.  If  M  Is  roughly  greater  than  20.  It 
seems  that  the  overflow  or  111 -conditioning  never  happens  in  the 
subroutine  Crout  (a  modification  of  Gauss  reduction  to  solve  a  set  of 
linear  system  equations)  as  long  as  M  is  roughly  less  than  20  if  the  VAX 


11/780  machine  Is  used.  Therefore  the  parameter  M  Is  always  taken  to  be 
less  than  20. 

H.  POLES  EXTRACTED  VIA  A  WINDOW  TECHNIQUE 

In  some  situations  the  data  available  may  be,  1)  from  too  large  a 
spectral  span  to  be  covered  with  one  rational  function  approximation,  2) 
without  low-frequency  content  due  to  difficulties  In  measurement,  or  3) 
Inaccurate  over  the  low  frequency  portion  If  calculated  data  are  being 
used  because  of  an  Inherent  limitation  of  the  method,  e.g.,  Uniform 
Geometrical  Theory  of  Diffraction  (UTD).  Thus,  a  window  technique 
{rectangular  window)  has  been  developed  for  extracting  poles  in  the 
above  situations. 

This  rectangular  window  Is  generated  simply  by  taking  sample  points 
over  a  band  limited  region  of  the  spectrum.  Then  the  RFA  Is  fitted 
accurately  over  only  that  region  and  the  extracted  poles  (whose 
Imaginary  component  falls  In  the  appropriate  range)  are  accurate.  Model 
II  given  below  is  used  In  the  window  technique. 

n  N  M 

En(X1)  =  l  anXj  /I  +  l  bm  Xf  1  =  1,  2 . I  .  (3-13) 

n=0  m=l  1 

2 

where  now  an  X^  In  the  numerator  Is  taken  out  since  there  Is  no  Rayleigh 
region.  It  Is  noted  that  due  to  a  truncation  of  data  on  both  sides  of 
the  window,  the  lowest  and  highest  frequency  of  oscillation  poles 
extracted  via  this  technique  are  not  very  accurate.  However,  the  rest 
within  the  window  are  still  very  close  to  the  true  poles. 
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Plots  using  this  window  technique  applied  to  the  backscatterl ng 
data  of  a  conducting  sphere  are  shown  In  Figures  3-11  to  3-13.  These 
plots  show  that  curves  within  the  window  fit  well  both  In  amplitude  and 
phase.  In  Table  3-7,  the  first  two  pairs  of  poles  are  extracted  via 
window  I  (ka  of  0.36  to  2.4),  the  third  and  fourth  pairs  of  poles  are 
extracted  via  window  II  (ka  of  2.1  to  4.2)  and  the  last  two  pairs  of 
poles  are  from  window  III  (ka  of  3.9  to  5.9).  Obviously,  as  the  windows 
move  to  higher  frequencies,  the  percentage  errors  are  still  In  the  same 
order.  The  rational  function  model  Is  Inherently  more  accurate  at  low 
frequencies  than  at  high  frequencies.  This  was  Initially  suggested  by 
Kennaugh  and  Moffatt  [2]  who  noted  that  for  a  distributed  parameter 
system  a  lumped  parameter  approximation  can  best  be  utilized  where  the 
scattering  or  radiation  properties  of  the  object  are  essentially  global, 
l.e.,  not  dominated  by  portions  of  the  object.  A  more  subtle 
distinction  Is  made  by  examination  of  the  K-pulse  [6]  which  notes 
something  missing  from  the  singularity  expansion  method  [SEM]. 

The  width  of  the  window  Is  chosen  such  that  at  least  three  pairs 
of  poles  are  covered  In  the  spectrum.  The  first  and  the  last  pair 
of  poles  are  deleted  because  the  truncation  of  data  makes  these  poles 
Inaccurate.  Also  note  that  the  poles  extracted  using  one  RFA  over  the 
whole  spectral  span  (one  window)  have  almost  the  same  order  of  accuracy 
as  those  extracted  using  three  RFA  and  three  overlapped  windows  over  the 
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TABLE  3-7 

A  COMPARISON  OF  TRUE  POLES  AND  POLES*  EXTRACTED  VIA 
WINDOW  TECHNIQUE 


Exact  Poles  of  Sphere  Poles  Extracted  via  RFA  A 

Windows  {Three)  Technique 


Percentage  Error 
Real  Part  Imag.  Part 


-0. 500000+/- JO. 866025  -0. 5003+/- jO. 866 10  (WI)** 

3.80E-2  7.50E-3 


-0. 701964+/- jl. 80740 

-0. 7086+ /-Jl. 8133 

(WI) 

0.342 

0.304 

-0. 842862+/- j2. 75786 

-0. 8766+/- j2. 7483 

(WII) 

1.17 

0.332 

-0. 954230+/- j3. 71478 

-0. 9277+/- j3. 7388 

(WII) 

0.692 

0.626 

-1.04764  +/-J4. 67641 

-0.9996+/-J4.6191 

(Will) 

1.002 

1.196 

-1.12891  +/-J5. 64163 

-1. 1338+/- j5. 5328 

(Will) 

8.50E-2 

1.891 

*  The  units  of  the  poles  are  In  lea,  where  k  Is  the  wavenumber  and 
a  Is  the  sphere  radius. 

**  WI  denotes  window  I. 
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Figure  3-11. 


The  RFA  fit  to  the  calculated  data  for  the  sphere 
within  the  frequency  of  Interest  (window  I).  The  x's  are 
22  used  for  the  RFA.  The  amplitude  plot  of 

the  RFA  (dashed  line)  gives  a  good  fit  to  the  original 
amplitude  data  (solid  line)  within  the  window. 
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Figure  3-12.  The  RFA  fit  to  the  calculated  data  for  the  sphere 

within  the  frequency  of  Interest  (window  II).  The  x’s  are 
the  data  points  used  for  the  RFA.  The  amplitude  plot  of 
the  RFA  (dashed  line)  gives  a  good  fit  to  the  original 
amplitude  data  (solid  line)  within  the  window. 
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Figure  3-13.  The  RFA  fit  to  the  calculated  data  for  the  sphere 

within  the  frequency  of  Interest  (window  III).  The  x's  are 
the  data  points  used  for  the  RFA.  The  amplitude  plot  of 
the  RFA  (dashed  line)  gives  a  good  fit  to  the  original 
amplitude  data  (solid  line)  within  the  window. 


whole  spectral  span.  The  detail  of  this  application  to  some 
experimental  data  will  be  demonstrated  In  a  later  chapter. 

I.  ADDITIVE  NOISE  AND  ITS  EFFECTS  ON  POLES  EXTRACTED  VIA  RATIONAL 
FUNCTION  APPROXIMATIONS 

Up  to  this  point,  the  rational  function  approximation  has  been 
applied  to  extract  the  complex  natural  resonances  from  sets  of  noiseless 
data.  However,  In  any  measured  system  noise,  clutter  and  uncertainty 
are  Inevitable.  Two  types  of  pseudo-noise  are  considered  In  this 
section.  One  type  Is  simply  the  calculated  data  rounded  to  the  kth 
decimal.  The  other  type  Is  an  uncorrelated  Gaussian*  white  noise  added 
to  the  calculated  data.  For  the  first  test,  both  the  real  and  imaginary 
part  of  the  calculated  sphere  data  are  rounded  to  the  kth  decimal. 

The  rational  function  routine  Is  applied  to  these  sets  of  data  with 
system  order  for  (M,N)  as  (14,13).  This  set  of  system  orders  Is  chosen 
because  the  system  with  these  orders  works  reasonable  well  (see  Figures 
3-7  and  3-8)  for  eight  poles  In  the  spectrum,  l.e.,  about  twice  the 
number  of  poles.  For  a  general  problem,  an  estimate  of  the  number  of 
poles  can  be  obtained  from  the  measured  amplitude  spectrum.  Tables 
3-8  to  3-9  show  comparisons  of  extracted  poles  from  data  rounded  to  the 
2nd,  and  to  the  1st  decimal  respectively.  Obviously,  as  k  decreases, 
the  extracted  poles  move  away  from  the  true  ones. 

*  The  Gaussian  random  number  generating  subroutlne-Gauss  In  IBM 
Scientific  Subroutine  Package  (SSP)  Is  used. 


TABLE  3-8 


A  COMPARISON  OF  TRUE  POLES  AND  EXTRACTED  POLES* 

USING  THE  RATIONAL  FUNCTION  APPROXIMATION  TO  THE  BACKSCATTERING 

DATA  OF  SPHERE** 


True  Poles  Extracted  Percentage  of  Error 


(M  ■  16t  N  *  15) 

Real  Part 

Imag.  Part 

-0.50000  +/- JO. 866025 

-0.5398+/- JO. 80699 

3.98 

5.903 

-0.701964+/-J1. 80740 

-0. 8233+/- jl. 8558 

6.257 

2.496 

-0.842862+/-J2. 75786 

-1.0036+/-J2.8927 

5.574 

4.676 

-0.954230+/-J3. 71478 

-0.9817+/-J3.7629 

0.716 

1.254 

(M  ■  14,  N  *  13) 

Real  part 

Imag.  Part 

-0 . 5196+/- JO . 8159 

1.96 

5.01 

-0. 7721+/- Jl .8882 

3.617 

4.167 

-0. 9394+/- J3. 1347 

3.348 

13.07 

-1.4014+/-J3.5855 

11.65 

3.37 

*  The  units  of  the  poles  are  In  kat  where  k  Is  the  wavenumber 
and  a  Is  the  sphere  radius. 

**  The  data  were  rounded  to  2nd  decimal . 


TABLE  3-9 


A  COMPARISON  OF  TRUE  POLES  AND  POLES  EXTRACTED 
USING  THE  RFA  TO  THE  BACKSCATTERING  DATA  FOR  A  SPHERE* 


True  Pole 


-0. 500000+/- jO. 866025 
-0. 701964+/- j 1.80740 
-0. 842862+/- j2. 75786 
-0. 954230+/- j3. 71478 


Extracted  Poles 
CM  =  16,  N  s  15) 

-0.4525+/-J1.012 
-0 . 4898+/- j  2 . 009 
-0. 3040+/- j3. 2430 
not  located 


Percentage  of  Error 
Real  Part  Imag.  Part 

4.75  14.59 
10.94  10.39 
18.69  16.82 


(M  *  22,  N  *  21)  Real  Part  Imag.  Part 

-0.415  +/-J0.8779  8.5  1.148 
-0. 5548+/- j 1.7920  7.59  0.794 
-0 . 56 74+ /- j 2 . 6454  9.55  3.899 
-0.4200+/-J3.7250  13.93  0.266 


*  The  data  were  rounded  to  1st  decimal 


The  second  type  of  noise  Is  uncorrelated  Gaussian  random  noise  with 


zero  mean  and  a  standard  deviation.  Two  uncorrelated  Gaussian  random 
processes  with  a  zero  mean  and  a  standard  deviation  of  a  which  are 
generated  from  two  different  seeds  are  added  to  the  real  and  Imaginary 
parts  respectively  of  the  calculated  data  for  the  sphere.  The 
slgnal-to-nolse  ratio  Is  defined  as 


S/N(dB) 


(3.14) 


where  A(X^)  Is  the  amplitude  of  signal  at  frequency  Xj,  and  N  Is  the 
total  number  of  samples  used. 

These  noisy  data  are  fitted  by  the  rational  function  In  the  exact 
sense.  Figure  3-14  Is  a  typical  example  of  rational  function  fit  to  the 
noisy  data,  where  the  noise  level  Is  a  of  l.E-1.  The  extracted  poles 
under  different  noise  levels  are  listed  In  Table  3-10.  It  Is  easy  to 
see  that  the  discrepancies  between  the  extracted  poles  and  the  true 
poles  are  Increased  as  the  noise  level  is  Increased.  Furthermore,  the 
rational  function  approxlmants  may  fall  to  extract  the  complex  natural 
resonances  at  all  when  the  S/N  Is  below  20  dB  or  so.  Based  on  the  above 
characteristics  of  noise  effects,  several  preprocessing  techniques  will 
be  used  to  preprocess  the  noisy  data.  These  preprocessing  techniques 
are  described  and  are  applied  to  several  sets  of  noisy  data  for  the 
sphere  In  the  following  section. 


TABLE  3-10 


COMPARISON  OF  EXACT  AND  EXTRACTED  POLES*  OF  SPHERE 
USING  THE  RFA  TO  THE  NOISY  DATA  FOR  THE  SPHERE** 


Poles  Extracted  via  RFA  (In  the  exact  sense) 
True  Poles  o=2.E-2(Run  1)  o®2.E-2(Run  2) 


-0. 500000+/- jO. 866025 
-0. 701964+/- jl. 80740 
-0. 842862+/- j2. 75786 
-0. 954230+/- j3. 71478 


-0.5362  +/-j0.9616 
-0.4745  +/-j 1.9577 
-0. 38518+/- j2. 643 
-0.4745  +/- j3 . 538 


-0. 4386+/- jO. 80129 
-0. 6482+/- jl. 8509 
-0. 9795+/- j2. 7384 
-0. 6285+/- J3. 5081 


o=l. E-l 


-0.439  +/-j0.5497 
-0.295  +/-J1.6033 
-0.2574+/-J2.784 
not  located 


Run  1  and  2  used  different  data  points. 

*  The  units  of  the  poles  are  in  ka,  where  k  Is  the 
wavenumber  and  a  Is  the  sphere  radius. 

**  Two  Gaussian  white  noise  generator  (using  two  different 
seeds)  with  standard  deviation  a  are  added  to  the 
backscatterl ng  data  for  the  sphere.  The  signal  to  noise 
ratio  corresponding  to  o  of  1 .E-l  Is  18.37  dB. 
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Figure  3-14.  The  RFA  applied  to  a  set  of  noisy  data,  where  the  Gaussian 
random  noise  with  zero  mean  and  a  standard  deviation  of 
0.1  Is  added  to  the  calculated  data  for  the  sphere. 

The  x's  are  the  data  points  used  for  the  RFA. 

The  amplitude  plot  of  the  RFA  Is  plotted  as  a  dashed 
line  (S/N  of  18.37  dB). 
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J.  PREPROCESSING  TECHNIQUES 


1.  Preprocessing  With  Zero-Phase-Shift  Digital  Filter 

In  general,  a  digital  filter  Is  defined  In  the  z-domaln  as 
N  M 

H(z)  =  \  a  z-1/  (1  +  )  b.  z-1)  ,  (3-15) 

1=1  n  1=1  1 

where  aj  and  are  coefficients  which  determine  the  behavior  of  the 
filter.  Emphasized  In  this  section  Is  the  design  of  a  high  order  (10) 
lowpass  digital  filter  with  zero-phase-shift,  where  the  3  dB  cutoff 
frequency  is  movable.  The  Butterworth  approximation  Is  used  because  It 
Is  a  maximal ly-fl at  amplitude  approximation.  In  order  to  obtain  a  10th 
order  zero  phase  shift  digital  filter,  a  5th  order  lowpass  Butterworth 

digital  filter  Is  designed  first.  Then,  a  complex  conjugate  filter 

*  T 

H  (z)  Is  cascaded  to  H(z).  Thus,  the  total  filter  H  (z)  Is  a  10th  order 

lowpass  filter  with  zero  phase  shift. 

HT(z)  =  H(z)  .  H*(z)  *  |H(z)|2  0°  .  (3.16) 

The  details  of  the  zero  phase  shift  digital  filter  and  its 

Implementation  are  described  In  Appendix  A.  Comparisons  of  the  filtered 

data  with  a  12.35  dB  S/N  (see  Figure  3-20  for  the  noisy  waveform)  and 

the  original  noiseless  data  are  plotted  In  Figures  3-15  (amplitude)  and 

★ 

3-16  (phase).  The  3  dB  cutoff  frequency  fc  is  set  at  20  Hz,  the  RFA  are 
applied  to  the  filtered  data.  The  extracted  poles  and  the  percent 


*  Corresponding  to  6.29  TAU  (diameter  transit  time)  of  the  Impulse 
response  (see  Appendix  A). 
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Figure  3-15.  A  comparison  of  amplitude  plots  for  the  calculated  data 
(solid  line)  and  the  filtered  data  (dashed  line)  for  the 
sphere.  The  filtered  data  are  obtained  by  filtering 
the  noisy  data  (calculated  data  with  additive  Gaussian 
random  noise)  using  a  10th  order  zero-phase-shift  digital 
filter.  The  3  dB  cutoff  frequency  fc  is  set  to  20  Hz. 
The  signal-to-nolse  ratio  Is  12.35  dB  before  filtering. 


PHASE  COEGI 


Figure  3-16.  A  comparison  of  the  phase  plots  for  the  calculated  data 
(solid  line)  and  the  filtered  data  (dashed  line)  for  the 
sphere.  The  filtered  data  are  obtained  by  filtering 
the  noisy  data  (calculated  data  with  additive  Gaussian 
random  noise)  using  a  10th  order  zero-phase-digital 
filter.  The  slgnal-to-nolse  ratio  Is  12.35  dB  before 
filtering. 
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errors  are  shown  In  Table  3-11.  Note  that  the  extracted  poles  are 
15  percent  In  error  In  the  real  part  and  10  percent  in  error 
In  the  Imaginary  part.  Therefore,  the  minimum  slgnal-to-nolse  ratio 
which  can  be  tolerated  Is  about  15  dB. 

2.  Preprocessing  With  Sum  Operator 

After  some  Initial  manipulation  the  technique  In  using  the  rational 
function  approximation  Is  to  solve  the  matrix  equation  AX  *  B,  where  A 
Is  an  mxn  matrix,  B  Is  an  mxl  matrix,  and  X  is  an  nxl  unknown  column 
matrix.  It  has  been  shown  that  the  least  squared  error  method  Is  very 
sensitive  to  noise  and  has  no  guarantee  of  convergence  If  an  Iterative 
search  Is  used  [25].  As  will  be  shown  in  Chapter  IV  ,  there  are  some 
advantages  to  solving  the  matrix  equation  In  an  exact  sense  rather  than 
solving  It  In  a  least  squared  error  sense.  However,  there  are  many  data 
points  which  contain  useful  Information  and  are  corrupted  by  noise  to  a 
certain  extent  but  are  unused. 

To  extract  better  poles  from  a  set  of  additive  noise  data,  a  major 
effort  Is  to  reduce  the  effect  of  noise  to  a  mlnlmun.  The  Intention 
here  Is  to  use  the  row  echelon  algorithms  (which  are  used  In  Gauss 
elimination)  for  an  overdeterml ned  matrix  equation  to  reduce  N  rows  to 
one,  l.e.,  N  linear  system  equations  can  be  scaled  arbitrarily  (say  1/M ) 
and  added  together.  A  detailed  description  of  the  sum  operator  Is  shown 
In  Appendix  B.  It  Is  expected  that  the  noise  In  each  element  of  the 
matrix  equation  Is  reduced  by  a  factor  of  about  1/M .  Assuming  that  X(t) 
Is  an  arbitrary  stationary  random  process  with  mean  7  and  variance  a. 


TABLE  3-11 


POLES  EXTRACTED  VIA  RATIONAL  FUNCTION  APPROXIMATION  AND  10th 
ORDER  ZERO-PHASE-SHIFT  DIGITAL  FILTER.* 


Extracted  Poles*  Percentage  of  Error 

_ Real _ Imaginary 

Run  1 


-0.415 

±  j 

0.802 

8.50 

6.40 

-0.364 

±  j 

1.526 

17.43 

14.51 

-0.361 

±  j 

2.744 

16.71 

0.48 

-0.307 

±  j 

3.738 

16.88 

0.61 

Run  2 


-0.637  ±  j  0.798 
-0.752  ±  j  1.716 
-0.529  ±  j  2.982 


13.7 

2.58 

10.88 


6.80 

4.71 

7.77 


tverage 


-0.526  ±  j  0.800 
-0.558  ±  j  1.621 
-0.418  ±  j  2.863 


2.6 

7.42 

14.73 


6.60 

9.61 

3.65 


*  The  original  slgnal-to-nolse  ratio:  12.35  dB.  The  3  dB  cut-off 
frequency  Is  set  at  fc  of  20  Hz. 

**  The  units  of  the  extracted  poles  are  In  ka,  where  k  Is  the  the 
wavenumber  and  a  Is  the  sphere  radius. 


(3-17) 


then  the  following  statistics*  are  true: 

N  - 

E  (  l  X-j/N)  =  J 
1=1 

N  - 

V  (  l  Xi/N)=  o2/N  ,  (3-18) 

1=1 

where  E  and  V  are  the  expected  value  and  the  variance  respectively,  N  Is 
the  number  of  samples.  Thus,  the  variance  of  the  noise  Is  reduced  by  a 
factor  of  1/N. 

The  rational  function  approxlmant  with  the  sum  operator  was  defined 
In  this  study  (see  Appendix  B).  The  sum  operations  are  made  In  the 
matrix  equation  formed  by  the  RFA  Instead  of  In  the  sampled  data.  Then 
the  sum  operator  and  RFA  are  applied  simultaneously  to  a  set  of 
calculated  sphere  data,  where  two  Gaussian  white  noise  processes 
generated  with  two  different  seeds  are  added  to  the  real  and  Imaginary 
parts  respectively.  It  is  noted  that  there  Is  a  tradeoff  between  the 
noise  reduction  factor  and  the  number  of  samples  (N)  used  for  the  sum 
operation.  The  noise  term  In  each  element  of  matrix  A'  (see  Appendix  B) 
Is  enlarged  by  a  factor  X^  ,  where  X^  Is  jk^a,  M  is  a  system  order.  If 

the  N  sampled  frequencies  Xj,  X2,  .  .  .  Xjj  used  for  a  group  summation 
are  assumed  to  be  close  to  each  other,  l.e., 

Xi  s  X2  2  X3  .  .  .  2  XN  ,  (3-19) 


*  John  Neter  and  William  Wasserman,  Applied  Linear  Statistical  Models, 
Richard  0.  Irwin,  Inc.  Homewood,  Illinois,  1974. 


where  2  means  for  almost  equal  to,  then  the  noise  reduction  factor  Is 
1/N  based  on  Equations  (3-17)  and  (3-18).  However,  If  the  number  of 
samples  used  In  the  group  summation  N  Is  large  and  Equation  (3-19)  Is  no 
longer  true,  then  the  additive  noise  reduction  factor  Is  no  more  a 
factor  of  1/N.  Nevertheless,  too  many  samples  as  well  as  too  few 
samples  In  the  group  of  sum  operations  still  reduce  the  additive  noise 
significantly.  Note  the  sum  operator  developed  here  does  not  reduce  any 
of  the  signal  Itself.  In  other  words,  the  sum  operator  can  be  applied 
to  a  noiseless  matrix  equation,  the  solutions  with  or  without  a  sum 
operator  are  exactly  the  same. 

The  manipulation  of  a  sum  operator  In  the  elements  of  the  matrix 
equation  Is  definitely  better  than  that  In  the  sampled  data  because  the 
sum  operator  defined  here  does  not  reduce  the  signal  at  all.  It  should 
be  noted  that  there  Is  an  optimal  N  so  that  the  noise  can  be  reduced  as 
much  as  possible.  However,  the  search  of  the  optimal  N  is  beyond  the 
scope  of  this  dissertation.  The  number  of  samples  (N)  used  in  the  group 
summation  Is  set  optionally  at  10  In  this  dissertation. 

The  plots  shown  In  Figures  3-17  and  3-18  (for  S/N  as  18.37  dB), 
Figures  3-19  and  3-20  (for  S/N  as  12.35  dB),  and  Figure  3-21  (for  S/N 
as  6.33  dB)  are  amplitude  plots  of  the  original  noisy  data  (solid  line) 
and  of  the  fitted  curve  via  sum  operator  In  the  rational  function  fit 
(dashed  line).  Various  sets  of  data  points  are  selected  from  the  first 
two  noisy  data  sets  (S/N  of  18.37  dB  and  12.35  dB)  respectively  (see 
Figures  3-17  to  3-20)  and  then  the  rational  function  approximations  and 
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17.  The  sum  operator  and  the  RFA  applied  slmultaneouly  to 
the  noisy  data  (solid  line)  for  the  sphere,  where  two 
Gausslanly  distributed  random  noises  have  been 
added  to  the  calculated  data  for  the  sphere. 

The  x's  are  the  data  points  used  for  the  RAF.  The 
clgnal-to-noise  ratio  Is  18.37  dB  (amplitude  plot). 
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Figure  3-18.  The  sum  operator  and  the  RFA  (dashed  line)  applied 

simultaneously  to  the  noisy  data  (solid  line)  for  the 
sphere,  where  two  Gausslanly  distributed  random  noises 
have  been  added  to  the  calculated  data  for  the  sphere. 
The  data  points  (the  x's)  are  different  from  those  In 
Figure  3-17.  The  slgnal-to-nolse  ratio  Is  18.37  dB 
(amplitude  plot). 
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Figure  3-19.  The  sum  operator  and  the  RFA  applied  simultaneously  to 
the  noisy  data  (solid  line)  for  the  sphere,  where  two 
Gausslanly  distributed  random  noises  have  been 
added  to  the  calculated  data  for  the  sphere.  The  x's  are 
the  data  points  used  in  the  RFA.  The  slgnal-to-nolse 
ratio  Is  Is  12.35  dB  (amplitude  plot). 
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Figure  3-20.  The  sum  operator  and  the  RFA  are  applied  simultaneously 
to  the  noisy  data  (solid  line)  for  the  sphere,  where 
two  Gausslanly  distributed  random  noises  have  been  added 
to  the  calculated  data  for  the  sphere.  The  selected 
data  points  (the  x’s)  used  In  the  RFA  are  different  from 
those  In  Figure  3-19.  The  slgnal-to-nolse  ratio  Is  12.35 
dB  (amplitude  plot). 
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Figure  3-21.  The  sum  operator  and  the  RFA  (dashed  line)  applied  to 
the  noisy  data  (solid  line)  for  the  sphere,  where  two 
Gausslanly  distributed  random  noises  have  been 
added  to  the  calculated  data  for  the  sphere.  Note  the 
selected  data  points  are  marked  by  x's.  The 
slgnal-to-nolse  Is  6.33  dB  (amplitude  plot). 


the  sun  operators  are  simultaneously  used  to  form  system  matrix 
equations.  The  extracted  poles  are  listed  In  Tables  3-12  and  3-13. 

Averaged,  the  real  part  has  13.64  percent  error  and  the  Imaginary 
part  4.81  percent  error  for  the  S/N  at  12.35  dB.  For  a  S/N  of  18.37  <ffi, 
the  real  part  has  10.5  percent  error  and  the  Imaginary  part  2.4  percent 
error.  The  percentage  error  Is  18.25  percent  for  the  real  part  and  9.6 
percent  for  the  Imaginary  part  for  a  S/N  of  6.33  dB.  The  correspondl ng 
amplitude  plot  using  a  sum  operator  Is  shown  In  Figure  3-21  for  S/N 
ratio  of  6.33  dB.  The  Imaginary  parts  have  been  Improved  more  than  the 
real  part  using  the  sum  operator.  The  fact  that  the  real  parts  of  the 
poles  are  very  sensitive  to  noise  Is  Is  well  known  In  the  extraction  of 
poles  from  the  transient  response  waveforms. 

3.  Preprocessing  with  Both  Digital  Filter  and  Sum  Operator 

The  functions  of  both  digital  filter  and  sum  operator  have  been 
discussed  In  the  two  preceedlng  sections.  In  order  to  extract  the 
complex  natural  resonances  Individually  from  a  set  of  noisy  data,  the 
two  techniques  are  used  sequentially.  In  the  beginning,  the  original 
data  are  preprocessed  with  a  10th  order  digital  filter  with  the  3  dB 
cutoff  frequency  fc.  Then,  the  rational  function  approxlmant  with  sum 
operator  Is  applied  to  the  preflltered  data  for  a  fc  of  20  Hz.*  Two 
examples  are  shown  In  Figures  3-22  and  3-23,  where  the  amplitude  of  the 


*  The  same  reason  for  an  fc  of  20  Hz  Is  given  In  Section  J  of  this 


TABLE  3-12 


POLES  EXTRACTED  VIA  RATIONAL  FUNCTION  APPROXIMATION 
AND  SUM  OPERATOR 


Extracted  Poles*  Percentage  Error 

S/N=18.37  dB  (o=l.E-l) _ Real  Imaginary 


Run  1  (Figure  3-17) 

-0.373+/-J0.764 

12.7 

10.2 

-0. 581 7+/-J 1.8582 

6.20 

2.62 

-0.6737+/-J2.7798 

5.87 

0.76 

-0 . 4253+/-J  3 . 7322 

13.79 

0.454 

Run  2  (Figure  3-18) 

-0.4047+/-J0.7592 

9.53 

10.68 

-0.5807+/-J 1.6620 

6.25 

7.499 

-0.4970+/-J2.7006 

11.99 

1.985 

-0 . 3379+/-J  3 . 6641 

16.07 

1.321 

Run  3 

-0.353+/-J0.8167 

14.7 

4.932 

-0.564+/-J 1.873 

7.12 

3.383 

-0 . 6478+/-J2 . 7828 

6.76 

0.865 

-0.3852+/-J3.7031 

14.84 

0.305 

Average  Poles 

-0.3769+/-J0.78 

12.31 

8.6 

-0.5755+/-J1.798 

6.52 

0.485 

-0.6062+/-J2.7562 

8.21 

0.057 

-0 . 3828+/-J  3 . 700 

14.90 

0.385 

The  units  of  the  poles  are  In  ka,  where  k  Is  the 
wavenumber  and  a  Is  the  sphere  radius. 


TABLE  3-13 


POLES  EXTRACTED  VIA  RATIONAL  FUNCTION  APPROXIMATION 
AND  SUM  OPERATOR. 


Extracted  Poles 

S/N=12.35  dB 
( cr=2E-l ) 

Percentage  of  Error 
Real  Imaginary 

Run  1  (Figure  3-19) 

-0 . 2989+/ -JO . 8643 

20.11 

0.17 

-0.2300+/-J 1.6802 

24.34 

6.56 

-0. 1827+/-J2 .2712 

22.89 

16.88 

-0 . 2455+/-J  3 . 6937 

18.48 

0.55 

Run  2  (Figure  3-20) 

-0 . 4695+/-J0 . 9854 

3.05 

11.94 

-0 . 2530+/ -J  2.160 

23.15 

18.19 

-0 . 5323+/-J  2.5157 

10.77 

8.397 

-0. 4474+/ -J 3. 6547 

13.21 

1.57 

Run  3 

-0 .457 1+/-J  0 . 9676 

4.29 

10.16 

-0.4727+/-J2.088 

11.82 

14.47 

-0.2651+/-J3.296 

20.03 

18.66 

-0.1739+/-J3.674 

20.34 

1.063 

Average  Poles 

-0.405  +/-J0.9391 

9.5 

7.307 

-0.318  +/-J 1.976 

19.80 

8.70 

-0 . 3267+/ - J 2 . 6943 

17.90 

2.20 

-0.2889+/-J3.674 

17.35 

1.06 

S/N  *  6.33  dB 

Extracted  Poles 

Percentage  of  Error 

(Figure  3-21) 

Real 

Imagl nary 

-0.418  ±j  0.9177 

16.4 

5.9 

-0.173  ±j  1.903 

27.3 

4.93 

-0.274  ±j  2.153 

19.7 

20.97 

-0.584  ±j  3.974 

9.6 

6.78 

Average 

18.25 

9.65 

The  units  are  ka,  where  a  Is  the  sphere  radius 
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Figure  3-22.  The  sum  operator  and  the  RFA  are  applied  simultaneously 
to  the  preflltered  data  (solid  line)  for  the  sphere.  A 
10th  order  zero-phase-shift  digital  filter  Is  used  to 
prefllter  the  noisy  data,  where  two  Gausslanly  distributed 
random  noises  have  been  added  to  the  calculated  data  for 
the  sphere.  The  x's  are  the  data  points  used  In  the  RFA. 
The  original  signal -to-nolse  ratio  Is  12.35  dB  (amplitude 
plot). 
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Figure  3-23.  The  sum  operator  and  the  RFA  are  applied  simultaneously 
to  the  preflltered  data  (solid  line)  for  the  sphere.  A 
10th  order  zero-phase-shift  digital  filter  Is  used  to 
prefllter  the  noisy  data,  where  two  Gausslanly  distributed 
random  noises  have  been  added  to  the  calculate  data  for 
the  sphere.  The  x's  are  the  data  points  used  In  the  RFA. 
The  original  slgnal-to-nolse  ratio  Is  12.35  dB.  The  data 
points  used  are  different  from  those  In  Figure  3-22 
(amplitude  plot). 
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filtered  data  (solid  line)  and  the  fitted  curve  using  sum  operator 
(dashed  line)  In  the  rational  function  are  plotted.  The  extracted  poles 
are  listed  In  Table  3-14.  It  has  been  shown  (see  Table  3-14)  that  the 
first  three  pole-pairs  are  much  Improved  using  the  digital  filter  and 
sum  operator  sequentially.  However,  the  fourth  pole-pairs  may  not  be 
located  at  all  because  both  digital  filter  and  sum  operator  are  using 
only  one  sided  Information  to  predict  or  to  smooth  the  data  around  the 
high  frequency  end  (due  to  the  truncation  of  data  there).  The  smoothed 
data  near  the  high-frequency  end  are  therefore  not  accurate. 

K.  BRIEF  SUMMARY 

In  this  chapter,  rational  function  approxlmants  have  been  defined, 
developed  and  applied  to  noiseless  and  noisy  data  for  a  conducting 
sphere.  Two  major  problems,  critical  points  In  any  system 
Identification,  are  the  system  order  and  the  noise  problem.  Inherently, 
the  system  order  Is  not  a  very  critical  point  In  the  rational  function 
fit  as  long  as  the  system  order  M  Is  not  less  than  the  true  poles 
desired.  In  fact  the  system  order  can  be  roughly  determined  from  the 
amplitude  curve  resonances,  l.e.,  the  order  of  the  system  Is  roughly 
determined  by  the  number  of  peaks  In  the  amplitude  data.  However,  some 
curve  fitting  poles  and  some  true  poles  should  also  be  counted  although 
they  are  not  seen  clearly  In  the  amplitude  peaks  (If  there  are 
Interferences  from  different  scattering  mechanisms). 


TABLE  3-14 

POLES  EXTRACTED  VIA  RATIONAL  FUNCTION  APPROXIMATION 
USING  THE  DIGITAL  FILTER  AND  THE  SUM  OPERATOR  (S/N  =  12.35  dB) 


Poles*  Extracted  via  Percentage  Error 

Digital  Filter  and  Sum 

Operator  for  (M,N)  as  (16,15)  Real _ Imaginary 

Run  1  (Figure  3.22) 

M*14,  N=13 
-0.531+/-J0.926 
-0.446+/-J1.730 
-0.321+/-J2.784 
-0.246+/-J3.587 


Run  2  (Figure  3.23) 

-0.568+/-j0.915 

6.8 

4.89 

-0. 664+/- j 1.790 

1.96 

0.89 

-0.565+/-J3.001 

9.63 

8.43 

Run  3 

-0.392+/-J0.927 

10.8 

6.10 

-0. 455+/- j 1.927 

12.73 

6.17 

-0.453+/-j3.0o 

13.52 

11.17 

Run  4 

-0. 568+/- JO. 8235 

6.80 

4.25 

-0. 789+/- j 1.538 

4.49 

13.89 

-0.455+/-J2.902 

13.45 

4.99 

Average 

-0.515+/-J0.897 

1.499 

3.097 

-0.589+/-J1.759 

5.826 

2.455 

-0.429+/-J2.923 

14.35 

5.72 

*  The  units  of  the  poles  are  In  ka,  where  k  Is  the  wavenumber 
and  a  Is  the  sphere  radius. 
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3.10 

6.037 

13.21 

3.992 

18.08 

0.893 

18.46 

3.33 

LJ 


As  far  as  the  noise  problem  Is  concerned,  a  high-order  digital 
filter  has  been  developed  to  prefllter  the  additive  noise  and  a  sum 


operator  Is  used  In  the  rational  function  fit  to  reduce  both  the 
additive  and  multiplicative  noise.  As  shown  In  this  chapter,  the 
poles  still  can  be  extracted  and  Identified  for  a  S/N  ratio  at  12.35 
dB*,  if  a  10th  order  zero-phase-shift  digital  filter  and  a  sum  operator 
are  used.  Although  the  stud|y  In  this  chapter  Is  limited  to  the 
computer  generated  additive  Gaussian  white  noise.  Chapter  VI  will 
present  an  extension  of  using  the  10th  order  zero-phase-shift  digital 
filter  and  RFA  to  the  measured  data  using  the  reflectivity  measurement 
facility  at  the  ElectroSclence  Laboratory. 


*  Although  there  Is  an  example  of  S/N  of  6.33  dB  (Figure  3-21); 
however,  the  percentage  error  Is  relatively  high  (see  Table  3-13). 


CHAPTER  IV 


LINEAR  LEAST  SQUARED  ERROR  SOLUTIONS 


A.  INTRODUCTION 

It  was  shown  In  Chapter  III  that  the  RFA  with  preprocessing  can 
extract  the  complex  natural  resonances  reasonably  well  when  the  signal 
Is  contaminated  by  additive  noise.  It  would  appear  that  the  RFA  applied 
In  a  least  squared  error  sense  might  fit  the  data  better.  In  this 
chapter  several  techniques  using  least  squared  error  solutions  are 
Investigated. 

B.  EIGENANALYSIS 

An  Imhomogeneous  system  of  linear  equations 

AX  *  E  ,  (4-1) 

where  A  Is  a  real  m  x  (n+1)  data  matrix  of  rank  (n+1)  and  X  and  E  are 
(n+1)  and  m  rows  column  vectors,  respectively,  Is  an  overdeterml ned 
system  If  m  >  (n+1).  The  components  of  X  are  unknown  coefficients  and  E 
Is  the  error,  assumed  small.  The  squared  error  Is 

XTATAX-ETE  ,  (4-2) 
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where  YT  means  the  transpose  of  Y. 

There  are  N+2  ways  to  solve  the  matrix  equation  [16].  However,  the 
solution  with  the  constraint 

l  lail2  -  1  ,  (4-3) 

1*0 

the  elgenanalysls  method  [16]  ,  Is  the  solution  to  the  min  max  problems 
and  yields  the  minimum  total  squared  error. 

Elgenanlysls  has  been  used  to  minimize  the  total  squared  error  of 
the  matrix  equation  (Equation  (4-1))  and  has  shown  some  success  In 
simple  simulated  cases  [16].  In  addition,  noise  remains  a  serious 
problem.  The  following  procedures  are  used  to  solve  the  matrix  equation 
given  by  Equation  (4-1). 

1.  Find  the  eigenvalues  of  ATA  first. 

2.  Find  the  eigenvector  which  corresponds  to  the  smallest 
eigenvalue,  this  Is  the  minimum  total  squared  error. 

3.  The  eigenvector  corresponding  to  the  smallest  eigenvalue  Is  the 
solution  X  of  the  matrix  equation  In  the  sense  that  the  sum  of 
the  residuals  squared  Is  minimized. 

Several  efforts  have  been  made  to  extract  the  complex  natural 
resonances  from  rational  function  approxlmants  via  elgenanalysls. 

The  results  are  not  good  even  for  a  relatively  low  system  order  for 
(M,N)  as  (6,5)  and  noiseless  data.  With  noisy  data,  the  result  Is  even 
worse.  Assuming  no  mistake  In  theory,  the  error  must  be  due  to 


numerical  calculation.  In  particular,  the  calculation  of  the 
eigenvalues.  This  will  be  clear  In  the  next  section  on  singular 
decomposition  where  It  Is  shown  that  better  results  can  be  obtained 


C.  EIGENANALYSIS  VIA  SINGULAR  VALUE  DECOMPOSITION 

Given  a  real  m  x  (n+l)  matrix  A  with  m  >  (n+l),  the  matrix  A  can  be 
decomposed  Into  the  following  form. 

A  *  UEVT  ,  (4-4.1) 

where 

uTu  -  VV T  »  vTy  -  I(n+1)  ,  (4-4.2) 

and 

I  *  d1ag(oi,  ....  o( n+l)).  (4-4.3) 

The  matrix  U  Is  an  mxn  matrix  consisting  of  n  orthonormal 1 zed 
elgnevectors  associated  with  the  n  largest  eigenvalues  of  AAT.  The 
matrix  V  Is  an  (n+l)  x  (n+l)  matrix  consisting  of  the  (n+l) 
orthonormal  1 zed  eigenvectors  of  ATA  and  the  diagonal  elements  of  z  are 
the  non-negative  square  roots  of  the  eigenvalues  of  ATA  which  are  called 
singular  values.  It  Is  assumed  that 

oi  >  og  >  03  >  ...  >  o(n+l)  >  0  •  (4-4.4) 

If  the  rank  of  A  Is  r  and  r  Is  less  than  (n+l),  then 


°r+l  3°r+2 


•••.  °(n+l)  *  0 


(4-4.5) 


The  decomposition  of  matrix  A  In  Equation  (4-4.1)  Is  called  singular 
value  decomposition  (SVD).  To  compute  the  singular  value  decomposition 
of  a  given  matrix  A,  Householder  transformations  and  a  QR  algorithm  are 
used  In  the  literature.*  The  advantage  In  solving  least  squared  error 
problems  using  SVD**  Is  that  there  Is  no  necessity  to  compute  ATA,  which 
may  Involve  unnecessary  numerical  Inaccuracies. 

Singular  value  decomposition  routines  are  available  In  the 
literature.***  An  application  of  the  singular  value  decomposition  to 
find  the  complex  natural  resonances  from  the  calculated  data  for  a 
sphere  Is  presented  In  Table  4-1  and  Table  4-3.  Table  4-1  Is  a  list  of 
singular  values  for  the  rational  function  fit  to  the  sphere  data  for 
(M,N)  as  (14,13)  and  I  as  20  (20  data  points).  Table  4-2  Is  a  list  of 
the  ratio  of  two  adjacent  singular  values.  Table  4-3  Is  a  list  of  the 
extracted  poles  from  the  sphere  backscatterl ng  data  via  SVD. 

The  smallest  slngualar  value  Is  of  the  order  of  0.01.  The  ratios  of 
two  adjacent  singular  values  are  all  In  the  order  of  100  to  10,  l.e.,  no 
sudden  change  occurs.  Thus  there  Is  no  clue  for  finding  the  system 
order  via  the  ratio  of  two  adjacent  singular  values  although  the 


*  G.H.  Gould  and  C.  Relnsch,  "Singular  Value  Decomposition  and  Least 
Squares  Solutions",  Numerical  Mathematics  14,  p.  403-420,  1970. 

**  The  eigenvector  of  A^A  which  corresponds  to  the  smallest  singular 
value  of  A  Is  the  solution  of  the  unknown  matrix  X  In  Equation  (4-1). 

***C.L.  Lawson  and  R.  J.  Hanson,  Solving  Least  Squares  Problems, 
Englewood  Cliffs,  N.  0.,  Prentice-Hail ,  1974. 
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TABLE  4-2 


THE  RATIO  OF  TWO  ADJACENT  SINGULAR  VALUES* 


THE  1TH  SINGULAR  VALUE  CMR  THE  2  TH  0.286972E+01 
THE  2TH  SINGULAR  VALUE  CVER  THE  3  TH  0.153009E+02 
THE  3TH  SINGULAR  VALUE  CVER  THE  4  TH  0.187740E+01 
THE  4TH  SINGULAR  VALUE  OVER  THE  5  TH  0.773612E+01 
THE  STH  SINGULAR  VALUE  OVER  THE  6  TH  0.209978E+01 
THE  6TH  SINGULAR  VALUE  OVER  THE  7  TH  0.120647E+02 
THE  7TH  SINGULAR  VALUE  (HER  THE  8  TH  0.200734E+01 
THE  STH  SINGULAR  VALUE  OVER  THE  9  TH  0.277767E+C1 
THE  9TH  SINGULAR  VALUE  OVER  THE  10  TH  0.151 40 5E+01 
THE  10TH  SINGULAR  VALUE  CVER  THE  11  TH  0.101432E+02 
THE  11TH  SINGULAR  VALUE  OV®  THE  12  TH  0.226 56 2E+01 
THE  12TH  SINGULAR  VALUE  OVER  THE  13  TH  0.303491E+01 
THE  13TH  SINGULAR  VALUE  OVER  THE  14  TH  0.154349E+01 
THE  14TH  SINGULAR  VALUE  OVER  THE  15  TH  0.188664E+01 
THE  15TH  SINGULAR  VALUE  OVER  THE  16  TH  0.214739E+01 
THE  167H  SINGULAR  VALUE  OVER  THE  17  TH  0.129813B+01 
THE  17TH  SINGULAR  VALUE  OVER  THE  18  TH  0.320356B+Q1 
THE  18TH  SINGULAR  VALUE  (MR  THE  19  TH  0.102635E+01 
THE  19TH  SINGULAR  VALUE  OV®  THE  20  TH  0.197989E401 
THE  20TH  SINGULAR  VALUE  OVER  THE  21  TH  0.178358E+01 
THE  21TH  SINGULAR  VALUE  CVER  THE  22  TH  0.393260E+01 
THE  22TH  SINGULAR  VALUE  (MR  THE  23  TH  0.117150  E+01 
THE  23TH  SINGULAR  VALUE  OVER  THE  24  TH  0.16 220 0E+01 
THE  24TH  SINGULAR  VALUE  OVER  THE  25  TH  0.220013E+01 
THE  25TH  SINGULAR  VALUE  (MR  THE  26  TH  0. 119880 E+01 
THE  26TH  SINGULAR  VALUE  (MR  THE  27  TH  0.494517E+01 
THE  27TH  SINGULAR  VALUE  OVER  THE  28  TH  0.216 41 8E+01 
THE  28TH  SINGULAR  VALUE  OVER  THE  29  TH  0.3017  91 E+01 


*  The  singular  values  In  Table  4-1  were  used. 


The  singular  values  In  Table  4-1  were  used 


TABLE  4-3 


A  COMPARISON  OF  THE  TRUE  POLES  AND  THE  EXTRACTED  POLES* 
USING  RFA  MODEL  AND  THE  SINGULAR  VALUE  DECOMPOSITION  TO 
THE  CALCULATED  DATA  {SPHERE,  BACKSCATTER) 


True  Poles 

-0. 500000+/- JO .8660250 
-0.701964+/-31. 80740 
-0.842862+/-j2. 75786 
-0. 954230+/- j3. 714780 


Extracted  Poles 

-0.501176+/-  jO. 8552016 
-0. 7111590+/- jl. 806740 
-0 . 8363640+/- J  2 . 76490 
-0.9193100+/-33. 67700 


Percentage  of  Error 
real _ Imaginary 

0.1176  1.082337 

0.4742292  0.0034 

0.2253303  0.2441246 

0.9104712  1.166823 


*  The  units  of  the  poles  are  In  lea,  where  k  Is  the  wavenumber 
and  a  Is  the  sphere  radius. 


extracted  poles  are  very  close  to  the  true  ones.  The  application  of 
singular  value  decomposition  to  noisy  data  falls  even  for  small  additive 
noise.  It  was  disappointing  to  learn  that  zero  phase-shift  digital 
filter  effected  no  improvement  In  the  extraction  of  poles  via  singular 
value  decomposition  of  the  filtered  data. 

0.  PSEUDO-INVERSE  METHOD 

Instead  of  forming  a  matrix  equation  such  as  Equation  (4-1), 
where  E  Is  an  error  matrix  for  an  overdetermined  system,  the 
other  simple  way  of  forming  a  matrix  equation  from  the  RFA  is  similar 
to  that  of  exact  determined  systems  in  Chapter  III,  l.e., 

AX  =  B,  (4-5.1) 

where  A  is  an  mxn  matrix  with  m>n,  X  is  a  nxl  unknown  matrix,  B  is  a  mxl 
matrix  and  the  elements  In  B  are  the  real  and  imaginary  parts  of  the 
complex  data  points.  Then  the  easiest  way  to  solve  this  equation  Is  to 
perform  a  a  pseudo-inverse  as 

ATAX  =  ATB  ,  (4-5.2) 

where  AT  Is  the  transpose  of  A.  In  fact  this  method  is  also  known  as  a 
least  squared  error  solution  without  weight  In  mathematics.  The 
application  of  this  solution  to  the  calculated  data  for  a  sphere  are 
given  In  Figures  4-1  (amplitude)  and  4-2  (phase).  Note  that  twenty 
complex  data  points  are  used  In  this  overdetermined  system  for  (M,N)  as 
(14,13)  (28  unknowns).  Table  4-4  Is  a  comparison  of  the  true  poles  and 
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TABLE  4-4 

A  COMPARISON  OF  THE  TRUE  POLES  AND  THE  EXTRACTED  POLES* 
USING  RFA  AND  PSEUDO-INVERSE  METHOD,  i.e.,  LEAST  SQUARED  ERROR 
WITHOUT  WEIGHT  TO  THE  CALCULATED  DATA  (SPHERE,  BACKSCATTER) 


True  Poles 


Extracted  Poles 


Percentage  of  Error 
real _ Imaginary 


-0.500000+/-J0. 8660250 

-0. 5043270+/- JO. 8679600 

0.4327 

0.1935 

-0. 701964+/- jl. 80740 

-0.70266 

+/-J1. 79561 

0.03589 

0.6081 

-0. 842862+/- j2. 75786 

-0.84936 

+/-J2.76335 

0.2253 

0.1903 

-0. 95423+/- j3. 71478 

-0.92042 

+/-j3. 70323 

0.8815 

0.30115 

*  The  units  of  the  poles  are  In  ka,  where  k  Is  the  wavenumber 
and  a  Is  the  sphere  radius. 
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Figure  4-1.  The  RFA  applted  to  the  calculated  data  for  the  sphere  In 
the  sense  of  least  squared  error  without  weight.  Note 
twenty  complex  data  points  are  used  (amplitude  plot). 
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Figure  4-2.  The  RFA  applied  to  the  calculated  data  for  the  sphere  In 
the  sense  of  least  squared  error  without  weight.  Note 
that  twenty  complex  data  points  are  used  (phase  plot). 
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the  extracted  poles  using  this  method.  The  extracted  poles  are  very 
close  to  the  true  ones  In  the  noiseless  situation. 

The  Intentions  to  use  more  data  points  have  been  made  but  the 
results  are  even  worse  due  to  the  numerical  errors  In  the  computation  of 
large  matrix.  An  attempt  to  extract  poles  using  the  pseudo-inverse 
method  from  sets  of  noisy  data  falls  even  for  small  additive  noise, 
e.g.,  signal  to  noise  ratio  at  30  dB.  Without  doubt  the  bias  problem  Is 
serious  In  this  method  as  well  as  the  other  two  methods  mentioned 
earlier. 

E.  SUMMARY 

All  techniques  discussed  In  this  chapter  concerned  methods  for 
solving  the  matrix  equation  In  the  least  squared  error  sense.  They  are 
elgenanalysls  using  the  eigenvalues  or  singular  values  and  a 
pseudo-inverse  method.  The  elgenanalysls  using  singular  values  works 
well  If  there  Is  no  noise;  however,  the  method  generally  falls  even 
for  slightly  noisy  data.  El"  ^analysis  would  appear  to  be  the  best  of 
these  unweighted  solutions.  The  anticipation  of  better  fits  from  least 
squared  error  solutions  using  more  data  points  Is  not  true  due  to  the 
accumulation  of  numerical  errors  and  the  bias  problem. 


CHAPTER  V 


AN  INVESTIGATION  OF  THE  COMPLEX  NATURAL  RESONANCES 
OF  VARIOUS  OBJECTS 


A.  INTRODUCTION 

It  was  shown  In  Chapter  III  that  the  RFA  could  be  applied 
successfully,  for  pole  extraction,  to  the  calculated  data  for  the  sphere 
with  or  without  additive  noise.  Furthermore,  It  was  shown  In  Chapter  IV 
that  the  RFA  used  for  pole  extraction  In  a  least  squared  error  sense 
were  very  sensitive  to  noise.  An  approximate  method  which  can  be 
utilized  to  extract  the  complex  natural  resonances  of  a  particular 
target  may  or  may  not  be  precisely  applicable  to  other  targets. 
Therefore,  various  sets  of  backscatterlng  data  for  different  targets 
should  be  tested  to  see  the  general  applicability  of  the  method.  In 
this  chapter,  the  complex  natural  resonances  of  some  simple*  conducting 
scatterers  such  as  a  circular  disc,  loop,  thin  wire,  finite  circular 
waveguide,  semi -Infinite  circular  wave-guide,  and  one  stick  model  of  a 
F104  aircraft  are  extracted  via  the  rational  function  approximations 

*  Applications  to  geometrically  complex  targets  are  given  In  Chapter  VI 
(measured  data). 
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(RFA)  In  the  exact  sense.  The  results  (poles)  are  compared  In  some 
cases  with  those  obtained  using  different  approximate  or  exact  (sphere) 
methods.  The  approximate  methods  (In  addition  to  RFA)  which  were  used 
in  this  study  are  the  Uniform  Geometrical  Theory  of  Diffraction  (UTD) 
[11]  for  a  disc  and  an  Integral  equation  formulation  and  numerical 
search  for  a  loop. 


B.  CONDUCTING  SPHERE 


1.  Analytical  Method 


For  convenience,  the  first  few  complex  natural  resonances  of  a 
conducting  sphere  as  taken  from  Stratton  are  listed  again  In  Table  5-1 
for  both  the  transverse  electric  (TE)  and  the  transverse  magnetic  (TM) 
modes.  Both  TE  and  TM  modes  are  excited  at  the  same  time.  The 
oscillations  of  the  electric  mode  are  Independent  of  the  magnetic  mode 
oscillations. 


2.  Signal -Flow-Graph  Method 


Using  the  differential  attenuation  and  phase  shift  of  a  surface 
wave  on  a  spherical  surface  given  by  Pathak  and  Kouyoumjlan  [11,12], 
Kennaugh  [10]  derived  two  characteristic  equations  for  a  geodesic  path 
around  the  sphere.  The  roots  of  these  two  equations  are  the  poles  of 
the  sphere.  The  first  4  pole-pairs  for  the  electric  mode  via 
signal -flow  graphs  [10]  and  rational  function  approximation  are  compared 
In  Table  5-2.  The  low  order  poles  extracted  via  the  rational  function 
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TABLE  5-1 

POLES*  OF  SPHERE  EXTRACTED  VIA  ANALYTICAL  METHOD 


Electric  Mode 


Magnetic  Mode 


■1.60+/-J0. 
■0.50+/-j0.86 
■2.17+/-J0.87 
0.70+/-J1.81 
0.83+/- j2. 77 


-1.  +/-J0. 

-1.5  +/-J0.86 

-2.26+/-j0 

-1.87+/-J1.75 


The  units  of  the  poles  are  In  ka,  where  k  Is  the  wavenumber 
and  a  Is  the  sphere  radius. 


TABLE  5-2 


A  COMPARISON  OF  THE  POLE-PAIRS*  FOR  A 
CONDUCTING  SPHERE  EXTRACTED  USING  DIFFERENT  METHODS 


True  Poles  [25] 
U. 50000  +/-J0. 86025 
0. 701964+/- jl. 80740 
0.842862+/-J2. 75786 
0.954230+/-J3. 71478 


Signal  Flow  Graphs 
■m;TO58+7-jOVB75T' 
-0. 7008+/- j 1.8130 
-0. 8422+/- j2. 7618 
-0 . 9538+/- j3 .7178 


Rational  Function 
'-O.'49999+7-JOBSS 
-0.7019  +/-j 1.8028 
-0.8424  +/-J2.7621 
-0.9175  +/-J3.6829 


*  The  units  of  the  poles  are  In  ka,  where  k  Is  the  wavenumber 
and  a  Is  the  sphere  radius. 


approxlmants  are  more  accurate  than  those  extracted  via  the 
signal -flow-graph  method  but  the  asymptotic  approach  obtains  the  high 
frequency  poles  with  greater  accuracy. 

3.  Rational  Function  Approximation 

The  results  of  applying  the  rational  function  approximations  to  the 
backscatterl ng  data  of  sphere  [27]  have  been  shown  In  Figures  3-7  and 
3-8  for  a  system  order  for  (M,N)  as  (14,13).  Fourteen  data  points 
(marked  by  x's)  are  used  In  the  RFA  In  the  exact  sense.  A  comparison  of 
the  true  poles  [25]  and  those  extracted  via  the  rational  function 
approximations  and  the  signal -flow-graph  approach  are  listed  In  Table 
5-2.  The  data  used  for  the  rational  function  approxlmant  are  from 
ka  of  0.2  to  ka  of  4,  l.e.,  the  data  are  truncated  at  ka  of  4. 

Therefore,  It  Is  not  surprising  that  the  real  part  of  the  last  pole-pair 
Is  not  as  accurate  as  the  others, 

C.  DISC 

1.  Introduction 

The  diffraction  of  a  plane  electromagnetic  wave  by  a  thin, 
circular,  metallic  disc  remained  unsolved  until  the  rigorous  solution 
obtained  by  Melxner  and  Andrejewskl  In  1950*.  The  disc  Is  treated 


*  J.  Melxner  and  W.  Andrejwskl,  "A  Rigorous  Theory  on  the  Diffraction  of 
Plane  Electromagnetic  Haves  by  a  Perfectly  Conducting  Circular  Disc 
and  by  a  Circular  Aperature  In  a  Perfectly  Conducting  Plane  Sheet," 
Aner.  der  Phys.,  Vol.  7,  p.  157,  1950. 
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as  a  limiting  surface  In  an  oblate  spheroidal  coordinate  system  using  an 
eigenfunction  solution.  Numerical  results  of  fourteen  discrete 
frequencies  were  subsequently  obtained  by  Andrejewskl  In  1952  [30]. 
Flammer  [31]  derived  another  solution  using  oblate  spheroidal  vector 
wave  functions  In  1953.  Hodge  [32]  presented  Melxner's  solution  In 
In  terms  of  the  notation  Introduced  by  Flammer  [31]  and  the  far  field 
amplitude  and  phase  scattering  data  (ka  of  0.2(0.2)15.2)  for  plane  wave 
Incidence  at  broadside  were  calculated  using  his  development  of  an 
efficient  method  of  computing  the  spheroidal  eigenvalues  [32].  Using 
Hodge's  computer  programs,  J.  L.  LI  of  the  ElectroSclence  Laboratory  at 
The  Ohio  State  University  compiled  the  backscatterlng  data  (ka  of 
0.2(0. 2)  10.)  of  the  disc  at  angles  off  broadside  Incidence 
(e  of  15°(15°)90°)  for  the  Incident  electric  field  In  both  e  and  <|> 
polarizations  (see  Figure  5-1).  The  results  of  applying  the  RFA  to 
these  sets  of  data  and  UTD  approximations  to  the  resonant  modes  of  the 
disc  are  presented  In  the  following  section. 

2.  Broadside  Incidence 

Rational  function  approxlmants  have  been  developed  and  applied  to 
the  calculated  backscatterlng  data  of  a  circular  disc.  The  data  sets 
ure  In  the  range  of  ka  of  0.2  to  ka  of  15.2  for  broadside  Incidence. 
Cases  with  both  the  Incident  electric  field  and  the  Incident  magnetic 
field  parallel  to  the  plane  of  the  disc  are  Included.  These  data  sets 
are  fitted  with  the  rational  function  approxlmants  In  the  exact  sense. 


Figure  5-2.  The  RFA  for  (M,N)  as  (3,2)  fit  to  the  back scattering  data 
for  the  disc  at  broadside.  The  x's  are  the  data  points 
used.  The  solid  line  Is  the  calculated  data  and  the  dotted 
line  Is  the  RFA  (amplitude  plot).  The  3  Input  samples  are 
In  the  Rayleigh  region  and  are  marked  with  x's. 


Figures  5-3  through  5-6  show  several  examples  of  applying  the  rational 
function  approximation  (RFA)  program  to  the  real  frequency  data  of  the 
disc  for  plane  wave  Incidence  at  broadside.  Appendix  C  gives  a  detailed 
list  of  the  coefficients  of  the  RFA  using  system  order  for  (M,N)  as 
(14,13)  (Table  C-l)  and  the  corresponding  zeros,  poles  and  residues 
(Table  C-2)  by  applying  RFA  to  the  disc  backscatterlng  data  at  broadside 
for  plane  wave  excitation.  Mote  the  coefficients  of  the  RFA  at  higher 
orders  become  smaller  and  smaller.  In  order  to  estimate  the  accuracy  of 
the  rational  function  coefficients  for  the  disc  backscatterlng  data  for 
a  plane  wave  Incidence  at  broadside,  a  long  division  of  the  rational 
function  approxlmant  with  the  coefficients  listed  In  Table  C-l  was 
performed  with  the  result  In  the  form  of 

n  k  n+1 

E  (X)  -  l  cnX  ,  (5-1) 

n»l 

where  X  Is  jka  and  a  Is  the  radius  of  the  disc.  Is  shown  In  Equation 
(C-l).  It  approximates  closely  the  result  In  Equation  (C-2),  the  low 
frequency  approximation  given  by  Boersma*  (see  Appendix  C).  Figure  5-3, 
which  Is  similar  to  that  obtained  by  D.B.  Hodge,  Is  the  RFA  to  the  first 
resonant  region  [33].  The  order  (M,N)  for  the  system  as  (3,2),  l.e., 
only  three  Input  samples  (marked  with  x's  In  Figure  5-2)  are  used  In  the 
RFA.  Figure  5-3  and  5-4  are  excellent  fits  of  the  RFA  to  the  entire 
data  region  with  a  system  order  of  (M,N)  as  (14,13).  Poles,  zeros,  and 

*  Boersma,  J . ,  "Boundary  Value  Problems  In  Diffraction  Theory  and 
Lifting  Surface  Theory,"  Thesis,  Gromlnger,  Holland,  1964. 
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Figure  5-3.  The  RFA  (dashed  line)  for  (M,N)  as  (14,13)  fit  to  the 

backscatterlng  calculated  data  (solid  line)  for  the  disc  at 
broadside.  The  x's  are  the  data  points  used  In  the  RFA. 

The  last  Input  sample  is  at  ka  of  14.8  (amplitude  plot). 


M 


ft _ 9 


97 


o 


-2.5 


Figure  5-' 


a 


6 


a 


0 


-i— | - 1 - r 

-2.0 


-]  .5 


i - ; - 1 - r 

oa 
c 


-0.5 


CD 


O 

U“ 


o 

C\J 


o 

CD 


O 


c.o 


>.  A  comparison  of  the  poles  of  a  disc  extracted  using  RFA 
(c)  and  UTD  (a)  for  broadside  excitation. 


100 


real  part 

imaginary  part 

real  part 

imaginary  part 

-4.588581  X10" 

9.540524  X10n 

-5. 136516 

X10" 

-9.723603  X10" 

-4.S88586  X10° 

-9.S40520  X10® 

-5. 136S25 

X10° 

9.729801  X10° 

-1.147251  X10° 

8.001202  X10° 

6.064494 

X10-> 

4.090105  X10® 

-1 .  1*17256  X10° 

-8.001204  X10° 

6.071571 

X10'5 

'1.090107  X10# 

-2.959823  X10° 

5.158536  X10° 

-3.0066S7 

X10® 

-5.328921  X10° 

6.072063  X10-* 

-4.090106  X10° 

-3.006629 

X10° 

5 .328900  X10® 

6.071192  X10-* 

4.090109  X10° 

-1.135383 

X10° 

8.119683  X 1 0° 

-2.959807  X10° 

-5.158528  X10« 

-1.135389 

X10° 

ij.  1 19688  X10# 

-1.045099  X10° 

4.680569  X10° 

-1.043285 

X10° 

-4.888601  X10° 

-1.010567  X10° 

-1.927570  X10° 

-1.043308 

X10° 

4.888618  X10® 

-1.010567  X10° 

1.927570  X 10® 

-2.209201 

X10° 

-1.656605  X10# 

-1.045106  X10° 

-4.680579  X10® 

-2.209198 

X10° 

1.656607  X10° 

6.973723  X10*‘ 

0.000000  X10° 

-1.010578 

X10® 

-1.927627  X 103 

-1.B4061S  XI0° 

-1.046664  X10e 

-1.010577 

X10® 

1.927627  X10° 

-1340615  X10° 

1.046664  X10° 

-4.917520 

X  10~* 

-1.228591  X 1 0® 

-3^5 

>002  X10° 

0.000000  X10° 

-4.917522 

X10-' 

1.228591  X10# 

-1.44 

5699  X10-' 

0.000000  X10° 

6.973706  : 

X10-' 

0.000000  X10° 

-1.443693 

X10-* 

0.000000  X10® 

I 

9.0 

<ko) 


Figure  5-5.  The  RFA  for  (M,N)  as  (18,17)  fit  to  the  backscatterlng 
data  for  the  disc  at  broadside.  The  x's  are  the  data 
points  used.  The  solid  line  Is  the  calculated  data 
and  the  dotted  line  Is  the  RFA.  Note  the  last  Input 
sample  of  the  RFA  Is  at  ka  of  7.2  (amplitude  plot). 
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Figure  5-6.  A  comparison  of  the  poles  of  a  disc  extracted  using  RFA 
(o)  and  UTO  (a)  for  broadside  excitation. 
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residues  are  all  shown  In  the  amplitude  and  phase  plots  In  Figures  5-3 
and  5-4.  With  a  system  order  at  M  of  14,  7  pole-pairs  are  extracted. 
However,  only  5  out  of  these  7  pole-pairs  are  the  true  system  poles  (see 
the  criteria  In  Chapter  III). 

The  result  shown  In  Figures  5-5  Is  to  extract  more  poles  from  a 
shorter  frequency  range  to  see  If  secondary  dominant  poles  may  appear. 
Unfortunately,  the  poles  of  the  second  layer  are  too  weakly  excited  to 
be  extracted  via  RFA  using  these  real  frequency  data.  However,  as  was 
discussed  In  Chapter  III,  the  overdetermined  system  order  does  not 
Influence  the  extraction  of  the  real  physical  poles.  In  addition,  the 
UTD  [11]  and  an  Equivalent  Current  approach  (EC)  [29]  can  be  used  to 
predict  the  complex  natural  resonances  of  the  disc  at  broadside.  A 
detailed  description  Is  given  In  Appendix  D.  A  list  of  the  poles 
extracted  via  the  RFA  and  those  via  the  UTD  method  for  the  disc  at 
broadside  are  presented  In  Table  5-3.  Figure  5-6  Is  a  comparison  of 
poles  of  disc  for  broadside  using  RFA  and  UTD.  The  poles  of  string  I  In 
Table  5-3  from  the  RFA  are  an  edge  diffraction  mode.  The  one  pole  pair 
of  the  string  II  In  Table  5-3  Is  the  lowest  frequency  pole  of  the 
creeping  wave  mode.  This  can  be  Identified  because  the  Imaginary  part 
Is  around  ka  of  1,  l.e.,  2wa  of  X.  The  poles  of  the  disc  at  broadside 
from  these  two  methods  are  very  close  to  each  other  except  for  the  two 
lowest  frequency  pole-pairs.  It  Is  known  that  UTD  will  not  properly 
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predict  low  frequency  poles  because  It  Is  an  asymptotic  high  frequency 
approximation.  The  oscillatory  parts  of  the  poles  (string  I)  occur  In 
Increments  of  around  ka  of  3  (which  Is  very  close  to  n,  l.e.,  2a  of  X), 
because  the  diameter  (2a)  should  be  close  to  (n  +  1/2)  wavelength  at  the 
resonances  of  the  edge  diffraction  at  broadside.  A  list  of  the  first  26 
pole-pairs  for  the  disc  at  broadside  excitation  (edge  diffraction  mode) 
Is  shown  In  Appendix  0.  Note  the  first  three  low  frequency  pole-pairs 
come  from  the  rational  function  approximation  to  the  backscatterl ng 
calculated  data  for  the  disc  at  broadside.  The  higher  order  poles  for 
broadside  excitation  are  obtained  using  UTD  approximation.  The  UTD 
Equivalent  Currrent  Concept  (EC)  [29]  has  been  used  because  of  the 
caustic  field.  A  detailed  discussion  Is  presented  In  Appendix  D.  From 
the  UTD  approximation  for  broadside  excitation.  It  has  been  shown  (see 
Appendix  D)  that  the  equivalent  current  of  more  than  the  second  order 
diffraction  (Included)  forms  a  geometric  series.  Obviously,  this  Is  the 
edge  diffraction  mode.  However,  the  first  order  diffraction  which  Is 
dominant  and  Is  not  Included  In  the  geometric  series  forms  the  lowest 
resonant  mode  of  the  disc.  The  first  order  equivalent  current  on  the 
rim  Is  shown  In  Equations  (D-2)  and  (D-3)  of  Appendix  D.  Apparently, 
both  the  equivalent  electric  and  magnetic  currents  have  only  one  period 
of  variation  on  the  rim,  l.e.,  2ira  of  x.  Therefore,  the  imaginary  part 
of  the  poles  (-0.49+/-J1.21)  Is  around  one  unit  of  ka,  where  a  Is  the 
radius  of  the  disc  and  k  Is  the  wavenumber. 
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TABLE  5-3 


A  LIST  OF  POLE-PAIRS*  FOR  THE  DISC  BACKSCATTERING  AT  BROADSIDE 
USING  THE  RFA  AND  UTD  APPROXIMATIONS 


RFA 

String  I _ String  II 

1.  -2.22+/-J1.656  -0.49+/-J1.23 

2.  -1.04+/-J4.87 

3.  -1.16+/-J8.10 

4.  -1.24+/-J11.28 

5.  — 1 . 27+/- j 14-4 


UTD  method 


-1.04  +/-J5.0542 
-1. 1616+/- j8. 2116 
-1. 2418+/- jl 1.3611 
-1 . 3024+/- jl4 . 5075 


*  All  poles  are  In  units  of  ka,  where  a  Is  the  radius  of 
disc  and  k  Is  the  wavenumber. 


3.  Poles  of  the  Disc  at  Oblique  Incidences 


Here  the  RFA  Is  applied  to  the  calculated  backscatterl ng  data  of 

A 

the  disc  at  0  of  (15°  (15°)  90°)  with  ♦  polarization  (electric  field 
parallel  to  the  disc).  The  angle  0  Is  measured  from  broadside  Incidence. 
The  first  dominant  pole-pair  (see  Table  5-4)  Is  exactly  the  same 
as  that  obtained  at  broadside, l.e.,  0. 49+/- j 1.23.  However,  there 
are  more  poles  extracted  In  this  string  at  oblique  Incidence  (In¬ 
stead  of  only  0.49+/-J1.23  for  the  broadside  excitation).  The  oscilla¬ 
tory  parts  of  this  string  have  an  Increment  of  approximately  one  ka 
unit.  Obviously,  the  dominant  string  of  poles  (edge  diffraction  mode) 
at  broadside  Is  not  strongly  excited  for  off-broadside  Incidences. 

A  calculation  of  the  caustic  distances  for  waves  multiply  diffracted 
at  Ql»  02  (see  Figure  D-4)  Is  presented  In  Appendix  D.  The  caustic 
distance  for  the  Nth  order  diffraction  Is 


o..  s  -a(l+2Ns1n0) 
N  nTTTFTTsTnS 


(5-2) 


Obviously,  the  caustic  distances  of  each  multiple  diffraction  Is  not  the 
same  at  oblique  Incidences.  There  Is  no  way  to  form  a  geometrical 
series  at  angles  off-broadside  Incidence.  However,  the  caustic 
distances  are  the  same  for  all  diffractions  at  broadside  (theta  of  0 
degrees).  Actually  the  caustic  distances  are  -a  for  any  order  of 
diffraction  at  broadside.  Therefore,  the  resonant  modes  are  not  the 
same  for  broadside  and  angles  off-broadside  Incidence. 


TABLE  5-4 


A  LIST  OF  POLES*  FOR  THE  DISC  AT  ANGLES  OFF  BROADSIDE 
INCIDENCE  (BACKSCATTER)  USING  THE  RFA  AND  SENIOR'S 
CREEPING  WAVE  AT  EDGE-ON  INCIDENCE 


Rational  Function  Approximation 
8*15°  8=30° 


8=90° 

-0.49  +/-J1.23 
-0.614+/-J2.31 
-0. 688+/- j 3. 36 
-0.763+/-J4.42 
-0. 848+/- j 5. 40 
-0.925+/-J6.42 
-1.08  +/-J7.43 
-1.27  +/-J8.53 


Senior's  creeping  wave  node 
e»90° 

-0. 506+/- j 1.23 
-0.603+/-J2.28 
-0.675+/-J3.33 
-0.733+/-J4.37 
-0,78  +/-J5.40 
-0.83  +/-J6.43 
-0.87  +/-J7.46 
-0.90  +/-J8.48 
-0.94  +/-J9.50 
-0.97  +/-J10.52 


-0. 49+/- j 1.23 
-0.64+/-J2.33 
-1.38+/-J3.3 
-0.98+/-J4.82 

0*45° 

-0. 492+/- j 1.23 
-0.62  +/-J2.32 
-0.66+/-J3.20 
-0.67+/-J4.22 


-0. 49+/- j 1.23 
-0.61+/-J2.32 
-1.39+/-J3.71 
-0.95+/-J4.52 

0=60° 

-0.49  +/-J1.23 

-0.687+/-J2.31 

-0.71+/-J3.39 

-0.74+/-J4.40 

-0.83+/-J5.43 

-1.04+/-J6.43 

-0.92+/-J7.43 

-0.99+/-J8.50 


*  The  units  of  poles  are  In  ka,  where  k  Is  the  wavenumber  and 
a  Is  the  radius  of  the  disc. 


The  creeping  wave  around  the  rim  of  the  disc  Is  not  an  easy 
problem.  Two  papers  [34,35]  which  specifically  deal  with  the  creeping 
wave  on  the  rim  of  the  disc  can  be  found  In  the  literature.  In  this 
study,  the  creeping  wave  formula  developed  using  both  the  Geometrical 
Theory  of  Diffraction  (GTD)  and  measured  technique  In  [34]  Is  used. 

After  some  manipulation,  the  normalized  electric  field  (l.e.,  the 
squared  root  of  the  normalized  Radar  Cross  Section  (RCS),  the  RCS 
normalized  to  the  area  of  the  disc)  Is 

5*/6 

y  -1/2  -1/3  j  „  1/3 

En  *  2(1  •  •)  u  e  exp[-iru  +  -  u  ]  ,  (5-3) 

cw  n  2 

where  y  Is  0.8*  (actually  a  trial  value),  u  Is  jka,  and  a  Is  the  radius 
of  the  disc  (see  Appendix  E).  The  complex  natural  resonances  are 
obtained  by  Imposing  a  self con si  stent  condition  on  the  creeping  waves 
which  revolve  around  the  rim  of  the  disc,  l.e.,  after  one  revolution 
these  fields  should  exhibit  phase  coherence  with  the  fields  at  the 
starting  points.  Consequently,  the  zeros  of  the  following  equation  (see 
Appendix  E)  are  the  poles  of  the  creeping  mode  of  the  disc. 
i/3  „ 

e„  U  -2iru  s  i  ,  (5-4) 

where  u  Is  jka  and  a  Is  the  radius  of  the  disc.  A  list  of  the  dominant 
string  of  poles  at  angles  off  broadside  Incidence  extracted  from  RFA 
and  those  extracted  via  Senior's  creeping  wave  mode  [34]  are  listed  In 
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Table  5-4.  The  details  of  the  extraction  of  poles  using  Senior's 
creeping  wave  mode  are  presented  In  Appendix  E.  A  comparison  of  the 
scattering  by  the  disc  at  edge-on  Incidence  between  Senior's  analytical 
method  and  Hodge's  calculated  data  Is  presented  In  Appendix  F. 

Table  5-5  Is  a  list  of  poles  using  the  RFA  to  the  calculated  data 
for  the  disc  at  several  aspect  angles  (0°(15°)90°)  and  two  polarizations 

A  A 

(e  and  4),  the  natural  resonances  are  Indeed  Independent  of  the 
Incidence  and  the  field  polarization.  For  angles  off-broadside 
Incidence,  It  Is  the  creeping  wave  mode  which  Is  dominant  for  resonance. 
A  detailed  list  of  the  first  30  pole-pairs  of  the  disc  using  Senior's 
creeping  wave  formula  for  edge-on  excitation  Is  shown  In  Table  5-6. 
Figure  5-7  Is  a  plot  of  the  poles  of  the  disc  for  broadside  using  RFA 
to  the  calculated  data,  and  the  poles  of  creeping  wave  mode  using 
Senior's  analytical  method.  This  set  of  poles  (Figure  5-7)  may  be 
used  to  generate  a  K-pulse  waveform  [6]  for  the  backscatterlng  of  the 
disc. 


4.  Andre jewskl's  Magnetic  Near  Field  Data  of  the  Disc 

Backscatterlng  at  Broadside 

In  Andrejewskl 's  dissertation  [30],  fourteen  normalized  scattered 
magnetic  field  data  points  for  the  disc  backscatterlng  at  broadside  were 
plotted  In  the  complex  plane.  Here,  the  normalized  magnetic  scattering 
data  at  these  14  points  are  reproduced  at  ka  of  0.5(0. 5), 4.  and 
ka  of  4.(1.), 10.  Details  of  the  calculation  are  shown  in  Appendix  G. 
These  data  are  shown  fitted  with  a  RFA  in  Figure  5-8.  The  extracted 
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TABLE  5-5 


A  LIST  OF  POLES*  FOR  THE  DISC  USING  THE  RFA  TO  THE  CALCULATED 
DATA  FOR  DIFFERENT  ASPECT  ANGLES  AND  POLARIZATIONS.  THE  RFA 
ARE  FITTED  TO  ka  of  0.2  to  6  (30  SAMPLES) 


Rational  Function  Approximation 

_ poles _ 

A 

♦  Polar-  PI  P2  P3  P4 

Izatlon 

0.  -1.04+/-J4.87  -1.16+/-08.10  -1.24+/-J11.28 

-0.49+/-J1.23 

15.  -0.49+/-J1.23  -0.64+/-J2.33  -1.38+/-J3.30  -0.98+/-04.82 

30.  -0.49+/-J1.23  -0.61+/-J2.32  -1.39+/-J3.71  -0.95+/-04.52 

45.  -0.49+/-J1.23  -0.62+/-J2.32  -0.67+/-J3.40  -0.75+/-J4.30 

60.  -0.49+/-J1.23  -0.63+/-02.30  -0.97+/-J3.80 

75.  -0.49+/-J1.25  -0.61+/-2.33  -0.73+/-J3.35  -0.67+/-J4.50 

90.  -0.49+/-J1.23  -0.60+/-J2.31  -0.74+/-J3.42  -0.65+/-J4.28 


e  Polar¬ 
ization 

15.  -0.49+/-J1.22  -0.64+/-J2.27  -1.08+/-J3.30  -0.82+/-J5.02 
30.  -0.49+/-J1.23  -0.62+/-J2.33  -1.23+/-J3.52  -0.94+/-J4.57 
45.  -0.49+/-J1.23  -0.64+/-J2.34  -0.87+/-J2.90  -0.86+/-J4.80 
60.  -0.49+/-J1.23  -0.62+/-J2.32  -0.65+/-J3.21  -1.52+/-J5.40 
75.  -0.49+/-J1.23  -0.63+/-U2.32  -0.64+/-J3.24  -1.26+/-J5.80 


*  The  units  of  poles  are  In  ka,  where  a  Is  the  radius  of  the  disc 
and  k  Is  the  wavenumber. 
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POLES*  OF  DISC 

n 

TABLE  5-6 

USING  SENIOR'S  CREEPING  WAVE  MODE 

Poles 

< 

< 

- 

1. 

(-0.5057350, ±1.208781) 

’  j 

2. 

(-0.6030402, ±2. 282039) 

3. 

( -0 . 6746498 ,±3.331739) 

4. 

(-0.7328901, ±4. 370643) 

5. 

(-0.7826892, ±5. 403155) 

6. 

(-0.8265909, ±6. 431364) 

r: 

7. 

(-0.86607 54, ±7. 456467) 

*  < 

8. 

(-0.9021299, ±8. 479179) 

9. 

(-0.9354153, ±9. 499999) 

10. 

(-0.9664097, ±10. 51928) 

11. 

(-0.9954736, ±11. 53727) 

12. 

(-1.022884, ±12. 55416) 

i 

13. 

(-1.048859, ±13. 57012) 

i 

14. 

(-1.073574, ±14. 58525) 

15. 

(-1.097172, ±15. 59966) 

16. 

(-1.119781, ±16. 61343) 

1 

17. 

(-1.141471, ±17. 62662) 

18. 

(-1.162356, ±18. 63929) 

■ 

19. 

(-1.182501, ±19. 65149) 

20. 

(-1.201968, ±20. 66326) 

21. 

(-1.220811, ±21. 67463) 

22. 

(-1.239077, ±22. 68564) 

23. 

(-1.256811, ±23. 69632) 

24. 

(-1.274049, ±24. 70668) 

25. 

(-1.290820, ±25. 71675) 

a 

26. 

(-1.307164, ±26. 72656) 

27. 

(-1.323088, ±27. 73611) 

28. 

(-1.338675, ±28. 74543) 

29. 

(-1.353844, ±29. 75452) 

30. 

(-1.368698, ±30. 76340) 

4 

*  Units  of  ka,  where  a  is 

the  radius  of  the  disc. 
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rail  part 

laginary  part 

raal  part 

imaginary  part 
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-4.460859  XIO'1 
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o 
o 
•  — "1 

0.000000  X10° 
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1.023446  XIO*' 
-1.023446  XIO’1 

«'  A 
-I  \ 


M=  14 
N=  13 


f  \  /^v 

/  V/  N- 


Figure  5-8.  The  RFA  (dashed  line)  fit  to  Andrejewskl  s  data  (magnetic 
near  field  denoted  by  the  x's)  for  the  disc  at  broadside. 
The  x's  are  data  points  used,  and  the  dashed  line  Is  the 
RFA  (amplitude  plot). 


poles  are  listed  In  Table  5-7.  The  resultant  poles  are  very  close  to 
those  obtained  from  the  RFA  to  Hodge's  data  for  the  disc  (electric 
field),  also  given  In  Table  5-7  and  plotted  In  Figure  5-9.  The  great 
discrepancy  of  the  third  pole-pairs  Is  due  to  the  truncation  In 
Andrejewskl 's  data  at  ka  of  10.  The  extracted  complex  natural 
resonances,  are  accurate  using  the  RFA  both  for  electric  far  field  or 
magnetic  near  field. 

5.  Summary 

Two  strings  of  poles  for  back scattering  from  discs  have  been  found. 
One  Is  the  resonant  mode  due  to  the  multiple  diffraction  at  the  edges 
(rim)  of  the  disc,  the  other  Is  the  creeping  wave  mode  formed  by  the 
wave  creeping  around  the  circumference  of  the  disc.  For  broadside 
Incidence,  the  dominant  string  of  poles  Is  the  edge  diffraction  mode. 
Simultaneously,  the  creeping  wave  mode  also  exists.  However,  only  the 
first  pole-pair  of  the  creeping  wave  mode  can  be  extracted  by  RFA  at 
broadside  Incidence.  A  detailed  discussion  Is  presented  In  Appendix  D. 
The  extracted  poles  are  Independent  of  the  angle  of  Incidence  and 
polarization  (see  Table  5-5).  This  conclusion  comes  from  applying  the 
RFA  to  the  calculated  data  for  the  disc  at  angles  off  broadside 
excitation. 


TABLE  5-7 


A  COMPARISON  OF  THE  COMPLEX  NATURAL  RESONANCES*  FOR  A  DISC 
USING  RFA  TO  HODGE'S  DATA  AND  ANDRE J EHSKI *S  DATA 
(BACKSCATTERING  FOR  BROADSIDE) 


RFA  to  Hodge's  Data  RFA  to  Andrejewskl 's  Data 

(Electric  Far  Field) _ (Magnetic  Near  Field) 


-0.49+/-J1.23 

-1.04+/-J4.87 

-1.16+/-J8.10 

-1.24+/-J11.28 

-1.27+/-J14.4 

data  up  to  ka*15.2 


-0. 547+/- j 1.19 
-0.998+/-J4.62 
-1.38  +/-J8.14 
data  up  to  kaa10. 


The  units  of  the  poles  are  In  ka,  where  k  Is  the  wavenumber  and 
a  Is  the  disc  radius. 


O 


Figure  5-9.  A  comparison  of  the  first  3  pairs  of  poles  of  the  disc  — 

for  broadside  using  RFA  to  H  Ige's  calculated  data  IQ)  ■ — - 

and  Andrejewskl 's  data  (a).  Note  third  pole-pairs  using 

Andre jewskl 's  data  Is  not  very  close  to  those  using 

Hodge's  data,  because  Andrejewskl 's  data  Is  truncated 

at  ka  of  10. 
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D.  WIRE  SCATTERERS 

1.  Introduction 

The  electromagnetic  resonances  of  thin  wire  scatterers  have  been 
studied  for  a  long  time  [15].  Kennaugh  solved  the  problem  using  a 
signal -flow-graph  method  [10].  The  results  were  compared  with  those  from 
an  Impedance  method  by  Richmond*  with  excellent  agreement.  Since  the 
poles  of  the  thin  wire  scatterers  are  well  known  and  they  may  be  the 
basic  building  elements  for  more  complex  targets.  It  Is  wise  to  apply 
the  rational  function  approximations  to  the  thin  wire  scatterers  and 
compare  the  results  with  the  known  poles. 

2.  Data  Sets  For  Thin  Wire  and  Thick  Wire 

Computed  values  of  backscatterl ng  (and  bl static  scattering  as  well 
as  Induced  current  distributions)  for  thin  wire  and  thick  wire**  were 
made  In  1964,  using  a  program  developed  by  J.  Richmond,  for  L/X  from  0.1 
to  1.4  at  the  Increment  of  0.1  (14  data  points)  for  plane  wave 
Incidences  at  broadside  and  30°  from  end-on.  (L  Is  the  wire  length  and 
D  Is  the  wire  diameter.)  These  sets  of  phasor  signals  were  fitted  with 
RFA  using  system  orders  for  (M,M)  as  (12,11)  ,  l.e.,  12  out  of  the  14 
data  points  were  used  In  the  approximations.  The  selection  of  the  data 


*  J.H.  Richmond, "Digital  computer  solutions  of  the  rlglous  equations 
for  scattering  problems",  IEEE,  v.  57,  pp.  796-804,  August  1965. 

**  The  wire  length  (t)  to  wire  diameter  (D)  ratio  are  20  and  2000  for 
thick  wire  and  thin  wire  respectively. 
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points  Is  very  arbitrary;  however,  the  set  of  data  points  which  can 
yield  a  smooth  Interpolation  Is  used  because  there  are  some  errors  In 
the  data  sets.  These  errors  may  be  due  to  numerical  calculation  or 
Imperfect  modelling.  The  amplitude  plots  of  the  fits  to  the  thin  wire 
scatterers  are  shown  In  Figures  5-10  to  5-13  for  both  broadside  and  30° 
from  end-on  Incidence.  Fifty  data  points  are  reproduced  In  Figure  5-10 
to  5-13  using  the  coefficients  of  the  RFA  fitted  to  the  14  samples  of 
the  wire  scatterers  In  each  case.  These  are  good  examples  of  using  the 
RFA  for  complex  data  Interpolation  as  well  as  pole  extraction.  The 
pole-pairs  extracted  via  the  RFA  and  those  extracted  via  other  methods 
are  listed  In  Table  5-8.  As  Is  well  known,  at  broadside  Incidence  only 
even  modes  are  excited.  However,  for  30°  off  broadside,  both  even  and 
odd  modes  are  excited. 

E.  CIRCULAR  LOOP 

1.  Introduction 

Complex  natural  resonances  of  a  circular  metallic  loop  are 
extracted  using  a  space-frequency  Integral  equation  and  numerical 
computation.  Then,  the  backscatterlng  data  for  the  loop  were  generated 
for  broadside  Incidence  (ka  of  0.2(0.2)10.)  and  for  edge-on  Incidence 
(ka  of  0.2(0.04)5.2).  Using  these  sets  of  data,  poles  were  extracted  by 
applying  the  RFA  to  the  calculated  data  and  the  results  are  compared 
with  those  obtained  from  the  Integral  equation  formulation. 


TABLE  5-8 


A  LIST  OF  POLES  FOR  A  THIN-WIRE  AND  A  THICK-WIRE 
EXTRACTED  USING  RFA,  INTEGRAL  EQUATION  AND  NUMERICAL  SEARCH 
AND  SIGNAL -FLOW  GRAPHS  FOR  BROADSIDE  AND  30°  OFF  END-ON 
EXCITATION  (BACKSCATTER) 


Poles* 

Richmond  Rational  Function  Approxlmants 

L/D=2000  L/D=2000 


Broadside  Theta=30° 


-0.0128+/-J0.2398 

-0.0128+/-J0.239 

-0 . 0128+/ -JO . 2399 

-0.0175+/-J0.4881 

-0 . 01 74+/- JO . 4884 

-0.0208+/-J0.7366 

-0.0234+/-J0.9854 

-0.018+/-J0.743 

-0 . 01998+/-J0 . 735 

-0.0255+/-J 1.2344 
-0.0273+/-J1.4835 

-0 . 48588+/ -JO . 2825 

Signal  Flow  Graphs 
L/D=2000 

(Bagby  and  Kennaugh) 

-0. 0123+/- JO. 2389 
-0. 0172+/- jO. 4872 
-0. 0204+/- jO. 7360 
-0. 0229+/- jO. 9850 
-0. 0249+/- j 1.2341 
-0.0266+/-J1 .4834 


Rational  Function  Approxlmants 
L/D=20 

Broadside  Theta*30°. 


-0.03086+/-J0.2149 

-0.052128+/-J0.677 


-0.03015+/-J0. 21495 
-0.04958+/-J0.4469 
-0 . 05349+/ -J 0 . 6887 


*  L  Is  the  wire  length,  D  Is  the  wire  diameter,  and  a  Is  the 
units  of  the  poles  are  In  ka,  where  k  Is  the  wavenumber. 
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5-10.  The  RFA  for  (M,N)  as  (12,11)  fit  to  the  back scattering  data 
for  the  thln-wlre  with  a  wire  length-to-wlre  diameter  ratio 
of  2000  at  broadside  excitation.  The  x’s  are  the  data 
points  used  In  RFA.  The  dashed  line  Is  the  rational 
function  approximation  (amplitude  plot). 
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Figure  5-11.  The  RFA  for  (M,N)  as  (12,11)  applied  to  the 

backscatterlng  data  for  the  thln-wlre  with  a  wire 
length-to-wlre  diameter  ratio  of  2000  at  30°  from  end-on 
Incidence.  The  x's  are  the  data  points  used  In  RFA.  The 
dashed  line  Is  the  RFA  (amplitude  plot). 
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Figure  5-12. 


The  RFA  for  (M,N)  as  (12,11)  applied  to  the  backscatterlng 
data  for  the  thick-wire  with  a  wire  length  (L)-to-wlre 
diameter  (D)  ratio  of  20  for  broadside  excitation.  The 
x's  are  the  data  points  used  In  RFA.  The  dashed  line  Is 
the  RFA  (amplitude  plot).  The  peak  In  the  amplitude  of 
the  first  resonance  can  be  seen  In  the  RFA,  because  more 
data  points  (50)  are  reproduced  using  the  coefficients  of 
the  RFA. 
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-3.014495  XIO** 
-3.014495  XIO** 


6.83SSS4  XIO*' 
-6.886S54  XIO'1 
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4.469083  XIO*' 
-1.429911  XIO’1 
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Figure  5-13.  The  RFA  for  (M,N)  as  (12,11)  applied  to  the 
backscatterlng  data  for  the  thick -wire  for  a 
wife-length  (L)  to  wire  diameter  (D)  ratio  of  20  at 
30  from  end-on  Incidence.  The  x's  are  the  data  points 
used  In  RFA.  The  dashed  line  Is  the  RFA  (amplitude 
plot). 


121 


2.  Integral  Equation  and  Numerical  Search 


Because  of  the  rotational  symmetry  of  the  circular  loop  (Figure 
5-14),  the  Induced  currents  1(4)  can  be  expressed  In  terms  of  In  e^n4, 
where  In  are  constants.  Using  Harrington's*  Impedance  matrix  formulation 
for  the  circular  loop 

Zjj|n  =  0  (m*n)  ,  (5-5) 

and 


Z™  -  i”*b  IJ  W  1  Kn*l  -  tfe) 2  U 


where 


2ir  -jkb/.  .  2  4  +  (a)z  . 
K  =  1  f  e  4s1n  j  V 

n  27  °  - ; - 7T 

[W|*  («)2F 


d<t> 


(5-6) 


(5-7) 


A  program  was  written  to  find  the  zeros  of  Znn  In  Equation  (5-6) 
using  Muller's  method  of  Iteration  and  Simpson's  rule  of  Integration. 
All  poles  found  with  this  program  are  consistent  with  those  found  In 
[22].  Table  5-9  Is  a  list  of  poles  of  the  first  layer**  at  different 
wire  to  loop  ratios.  A  plot  of  those  poles  (Table  5-9)  Is  shown  In 


*  R.F.  Harrington,  Field  Computation  by  Moment  Method,  New  York, 
Macmillan,  1968. 

**  The  classification  of  Layers  are  shown  In  Figure  5-17. 
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TABLE  5-9 


A  LIST  OF  POLES*  (LAYER  I)  OF  LOOPS  AT  THREE  DIFFERENT 
WIRE  TO  LOOP  RADIUS  RATIOS.  THE  INTEGRAL  EQUATION  AND 
NUMERICAL  SEARCH  IS  USED  TO  FIND  THE  POLES. 


Poles 


a/b  ( .5E-3)  x  it 


(2.E-3)  x  it  (l.E-3)  x  it 


-.109  +/-J 1.0485 
-.164  +/-J2.066 
-.2085  +/-J3.079 
-.248  +/-J4.085 
-.2845  +/-J5.095 
-.3185  +/-J6.11 


-0.086  +/-J1.04 
-0.125  +/-J2.058 
-0.1562+/-J3.068 
-0 . 1825+/- J  4.078 
-0. 2068+/-J5. 06 
-0.2288+/-J6.09 


-0.073317+/-J1. 03822 
-0.10591  +/-J2. 05176 
-0.130826+/-J3. 06151 
-0.151932+/-J4. 06949 
-0.17066  +/-J5.0764 
-0. 187729+ /-J6. 0826 


*  The  units  of  the  poles  are  In  kb,  where  k  is  the  wavenumber  and 
b  is  the  loop  radius.  (  it  is  3.14159  and  a  is  the  wire  radius.) 
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Figure  5-14.  The  coordinate  system  for  a  circular  metallic  loop 


Figure  5-15.  Figure  5-16  Is  a  plot  of  poles  for  wire  radius  to  the  loop 
radius  ratio  of  3.14159E-3. 

The  numbers  shown  with  the  locations  of  poles  In  Figure  5-17  are 
n's  of  Znn.  The  poles  In  Figure  5-17  are  divided  into  layers.  Table 
5-10  Is  a  list  of  poles  for  the  circular  loop  via  Integral  equation  and 
numerical  search  at  wire  radius  to  the  loop  radius  ratio  of  3.14159E-3. 
Figure  5-18  Is  a  plot  of  all  poles  showing  constant  n  results  In  Znn. 
From  the  result  of  applying  RFA  to  the  backscatterlng  data  for  the  loop. 
It  will  be  seen  that  layer  I  poles  are  found  at  edge-on  incidence  and 
layer  II  poles  are  found  at  broadside  (also  the  first  one  of  layer  I). 
Layer  III  poles  are  not  found  using  RFA  approach  either  at  broadside  or 
at  edge-on  Incidence.  However,  this  does  not  mean  that  they  are  not 
there. 


3.  Rational  Function  Approximation 

The  backscatter  of  the  loop  for  broadside  and  edge-on  (TE  case) 
Incidence  (Figures  5-19  and  5-20)  were  calculated  using  Kouyoumjlan's 
equations  In  his  report*.  The  backscatter  of  both  broadside  and  edge-on 
Incidence  are  shown  In  Figures  5-19  and  5-20  Individually. 

The  variational  method  has  been  used  to  determine  the  radar  cross 
section  (RCS)  of  a  perfectly-conducting,  thin,  circular  loop  by  R.G. 


R.G.  Kouyoumjlan,  "The  Backscatterlng  from  a  circular  loop  (antenna)," 
Bulletin  Ho.  162,  The  Ohio  State  University,  The  Engineering 
Experiment  Station,  November  1956. 
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TABLE  5-10 


A  LIST  OF  POLES*  OF  A  LOOP  (a/b**  of  3.14159E-3)  EXTRACTED 
USING  INTEGRAL  EQUATION  AND  NUMERICAL  SEARCH 


layer  1 

-0.  +/-J0. 

-0.073317+/-J1. 03822 
-0.105910+/-J2. 05177 
-0.130827+/-J3. 06151 
-0.15193  +/-J4. 06950 
-0.170661+/-J5.0769 
-0.187729+/-J6.0826 
-0.203555+/-J7.0883 


_ Poles _ 

layer  3 

-1.48611  +/-J0. 
-1.47916  +/-J1. 54327 
-1.6727  +/-J2. 76692 
-1.82195  +/-J3. 90846 
-1.9425  +/-J5. 01545 
-2.048715+/-J6. 10285 

layer  4 

-2.15481+/-J0. 

-2.70814+/-J1.0550 

-3.12022+/-J2.05139 

-3 . 45635+/-J  3 . 02946 

-3 . 74378+/-J3 . 99996 

-3.99684+/-J4.96701 


layer  5 
-3.50331+/-J0. 
-4.10168+/-J0.97958 
-4 . 59325+/-J 1 .93514 
-5.0157  +/-J2. 88079 


layer  6 
-4.84712+/-J0. 
-5.4633  +/-J0. 94868 
-5.9950  +/-J1. 88392 
-6.4665  +/-J2. 81279 
-6.8921  +/-J3. 73865 


layer  2 

-0.073317+/-J1. 03822 
-1.47916  +/-J1. 54327 
-1.25879  +/-J5.0763 
-1.31026  +/-J8.3208 
-1.34771  +/— J 11.5139 


layer  7 
-6.1859+/-J0. 

-6.817  +/-J0. 93148 
-7. 3681+/-J 1.85427 

*  The  units  of  the  poles  are  In  kb,  where  k  Is  the  wavenumber  and  b  is 
the  loop  radius. 


a  Is  the  wire  radius 


Figure  5-15.  A  plot  of  the  complex  natural  resonances  (layer  I)  for 
the  loops  at  different  wire  to  loop  radius  ratios  using 
the  Integral  equation  and  numerical  search  to  find 
the  poles. 


Figure  5-16.  A  plot  of  the  poles  for  a  circular  loop  at  the  wire  radius 
to  the  loop  radius  ratio  of  3.14159E-3.  The  Integral 
equation  and  numerical  search  Is  used  to  find  the  poles. 
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Figure  5-17. 


A  plot  of  the  poles  for  a  circular  loop  {a/b  of 
3.14159E-3).  The  poles  are  extracted  via  the  Integral 
equation  and  numerical  search.  The  number  around  the 
location  of  each  pole  Is  the  n  of  Znn.  Three  types  of 
poles  are  defined. 
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Figure  5-18.  A  plot  of  the  poles  for  a  circular  loop  (a/b  of 

3.14159E-3  extracted  via  the  Integral  equation  and 
numerical  search.  The  number  around  the  location  of  each 
pole  Is  the  n  of  Znn.  The  poles  which  are  extracted  from 
the  same  Znn  are  Identified  by  a  solid  line. 


Figure  5-19.  The  application  of  RFA  for  (M,N)  as  (14,13)  to  the 

backscatterlng  data  (solid  line)  for  a  circular  loop  at 
broadside.  The  wire-radius  to  the  loop  radius  ratio  is 
3.14159E-3.  The  x's  are  the  data  points  used  In  the  RFA. 
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Figure  5-20.  The  RFA  (dotted  line)  fit  to  the  back scattering  data 

(solid  line)  for  the  loop  at  edge-on  Incidence,  a/b  Is 
3.14159E-3,  where  a  Is  the  wire  radius  and  b  Is  the  loop 
radius.  The  x's  are  the  data  points  used  In  RFA. 


Kouyoumjlan  at  the  Ohio  State  University.  The  equivalent  source  of  the 
scattered  field  Is  approximated  by  the  current 


I($)  *  1  (cncosn$  +  dns1nn$)  ,  (5-9] 

n=0 

where  I ( 4> )  Is  assumed  to  flow  parallel  with  the  axis  of  the  wire.  Good 
agreement  between  the  calculated  (broadside  and  edge-on)  and  measured 
values  were  obtained.  More  data  points  are  necessary  to  test  the  RFA. 
Therefore,  the  following  formula  were  used  to  calculate  additional 
backscatterlng  data  (TE  case)  In  this  study. 

En(kb)  =•  2  2  ir(kb)2  \  Jjj(*b  S<nfl) 

I  kzbzK. (kb  «) 

A  D 

•  2 
-  I  [Jn-l(kbslne)  -  Jn+l(kbs1n0)]  I 

+  l  (-1)"  n*l _ _ _ I 

^[K^jtkb,  |)  +  K^kb,  «))  -  2n*Kn(kb,  |)"l  (5-10 


and  at  broadside 


En( kb)  *  2  2w(kb)z _ 

|kzbz[KQ(kb,A)  +*2  (kb, a]  -  ^l(kb,A|z 


(5-11 


where  A»a/b,  and 


7  o,l/2 

*  cosn$  exp{-1kb[4s1n4*  +  (*rJ  } 

Kn(kb  «)  *  J  _ I_JL__d*  (5-12) 

^  0  o  ,  o,l/2 

[2sinZ*  *  1  («)Z3 
7  7  F 

En  Is  the  normalized  backscattered  amplitude  and  phase,  i.e.,  the 
squared  root  of  the  normalized  (wa2)  RCS,  and  Jn  Is  the  Bessel 
function  of  the  first  kind  and  order  n. 

The  backscatterlng  data  for  the  loop  at  broadside  and 
edge-on  (TE  case)  Incidences  were  calculated  using  Equations  (5-10)  to 
(5-12).  It  Is  noted  that  the  Kn  of  Equation  (5-12)  were  Integrated  by 
Simpson's  rule  for  5000  points  In  [0,*].  Also,  three  unequal 
subsections  were  divided  In  [0,*],  I.e.,  more  data  points  were  put  In 
the  steepest  region  of  Kn.  Thus,  more  accuracy  of  Kn  can  be  reached 
without  much  more  effort  In  the  Integration.  Fifty  data  points  (kb  of 
0.2(0.2)50)  were  calculated  using  the  Equation  (5-11)  for  the  broadside 
excitation. 

More  effort  Is  Involved  In  the  calculation  of  backscatterlng  data 
for  the  loop  at  edge-on  Incidence.  There  Is  a  summation  of  an  Infinite 
series  In  Equation  (5-10).  The  Infinite  series  was  truncated  at  the 
point  where  five  more  successive  terms  made  little  contribution 
to  the  total  sum  at  that  point.  Thus,  the  backscatterlng  data  of 
the  loop  at  edge-on  Incidence  were  calculated  at  ka  of  0.2(0.1)10. 
Figures  5-19  and  5-20  show  the  application  of  the  RFA  to  the 
backscatterlng  data  for  the  loop  at  broadside  and  edge-on  respectively 
(amplitude  plots). 


The  RFA  is  not  only  an  excellent  approximation  In  finding  the 
complex  natural  resonances  but  Is  also  a  good  approach  for  complex  data 
Interpolation.  Here  Is  a  good  example  of  using  the  RFA  for  the  data 
Interpolation  because  the  calculated  data  are  very  expensive  and 
time-consuming  to  obtain.  For  broadside  excitation  the  dominant  string 
of  poles  was  extracted  using  RFA  (Figure  5-19).  Table  5-11  is  a  list  of 
poles  (layer  I)  at  broadside  via  the  Integral  equation  (using  numerical 
search)  and  RFA.  It  Is  Interesting  to  note  that  all  of  these  poles  are 
the  same  as  those  searched  at  Zn  via  the  Integral  equation  and 
numerical  search.  From  the  RFA  viewpoint,  only  one  pair  of  poles  (the 
lowest  frequency  ones)  belong  to  layer  I.  The  rest  belong  to  layer  II. 

The  pole-pairs  of  a  loop  at  broadside  are  similar  to  the  pole-pairs 
of  a  disc  at  broadside.  For  edge-on  Incidence  (Figure  5-20),  the 
dominant  string  of  poles  Is  shown  on  Table  5-12.  They  belong  to  the 
family  Zll,  Z22,  Z33,  ...,  etc.  respectively.  This  string  of  poles  has 
the  same  spacing  In  the  Imaginary  (oscillatory)  part  as  that  of  creeping 
wave  mode  of  disc. 

Some  of  the  pole-loci  of  a  loop  are  similar  to  those  of  a  disc, 
because  some  of  the  scattering  mechanisms  are  the  same.  One  Is  due  to 
the  creeping  wave  around  the  cl rcumference  of  the  scatterer,  the  other 
Is  due  to  the  multiple  diffractions  of  the  thin  wire  (diffractions 
across  the  center  of  the  loop).  A  Table  of  25  pole-pairs  of  the  first 
layer  via  the  Integral  equation  and  numerical  search  are  also  shown  In 
Table  5-13.  Figure  5-21  Is  a  comparison  of  the  poles  extracted  using 
RFA  and  Integral  equation  with  numerical  search  for  the  loop. 


TABLE  5-11 


A  LIST  OF  POLES  OF  A  LOOP  (a/b*  of  3.14159E-3)  EXTRACTED  USING 
THE  RFA  AND  INTEGRAL  EQUATION  AND  NUMERICAL  SEARCH  FOR 
BROADSIDE  EXCITATION. 


Poles 

Integral  Equation 

and  Numerical  Search  Rational  Function  Approximation 


Layer  I 


-0 . 073317+/-J 1 . 03822  -0 . 073548+/-J 1 . 03834 


Layer  II 


-1.47916+/-J 1.54327 
-1.25879+/-J5.0763 
-1 . 31026+/-J8. 3208 
-1 . 34771+/— 11 . 5139 


-1.475459+/-J 1.540965 
-1.261053+/-J5. 07532 
-1.305286+/-J8. 33954 


*  The  units  of  the  poles  are  kb,  where  k  Is  the  wavenumber, 
b  Is  the  loop  radius,  and  a  Is  the  wire  radius. 


TABLE  5-12 


POLES  OF  A  CIRCULAR  LOOP  AT  EDGE-ON  INCIDENCE  WITH 
a/b*  of  3.14159E-3 


Integral  Equation 
and  Numerical  Search 


Rational  Function  Approxlmant 


Layer  I 


-0.73317E-1+/-J 1.03822  -0.7354E-1+/-J1. 03834 
-1.05310E-1+/-J2. 05176  -1.0628E-1+/-J2.0515 
-1.30826E-1+/-J3. 06151  -1.308  E-1+/-J3.06 


Layer  II 

-1.47546+/-J1. 540965  -1.62  +/-J1.69 


*  The  units  of  the  poles  are  kb,  where  k  Is  the  wavenumber, 
b  Is  the  loop  radius,  and  a  Is  the  wire  radius. 
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TABLE  5-13 


A  LIST  OF  THE  POLES  AND  THE  CORRESPONDING  RELATIVE  IMPEDANCE 
(Znn)  FOR  THE  LOOP  (a/b*  of  3.14159E-3).  THE  IMPEDANCE  {Znn) 
MAY  BE  AS  HIGH  AS  1.0  E3  FOR  THE  LOCATION  SLIGHTLY  AWAY  FROM  THE 

POLE  LOCATIONS. 


N 

P0LES(2kb) 

Z  (NN) 

0 

(0.,0) 

1.78470E-5 

1 

(-1.46634E-1,±2.07644) 

1.13116E-4 

2 

(-2.1 1820E-1 , ±4 . 10353 ) 

3.79588E-4 

3 

(-2.61653E-1,±6.12302) 

3.21737E-4 

4 

(-3. 03864E-1, ±8.13899) 

6.98320E-4 

5 

(-3.41321E-1,±10.1528) 

5.67344E-4 

6 

( -3 . 75458E-1 ,±12.1652) 

2.01817E-5 

7 

(-4.07110E-1,±14.1765) 

3.44889E-4 

8 

( -4 . 36814E-1 ,  ±16 . 1869 ) 

7.87555E-4 

9 

( -4 . 64939E-1 ,  ±18 . 1967 ) 

2.44959E-5 

10 

(-4.91752E-1,±20.2058) 

7.70241E-4 

11 

{ -5 . 17453E-1 , ±22 . 2145 ) 

7.71799E-4 

12 

( -5 . 42195E-1 , ±24 . 2228 ) 

6.32431E-4 

13 

(-5.66101E-1,±26.2307) 

7.59797E-4 

14 

(-5. 89269E-1, ±28.2382) 

9.36263E-4 

15 

(-6.11779E-1,±30.2455) 

3.61295E-5 

16 

( -6 . 33699E-1 , ±32 . 2525 ) 

1.81433E-3 

17 

( -6 . 55086E-1 , ±34 . 2593 ) 

8.89091E-4 

18 

( -6 . 75987E-1 ,±36.2658) 

5.91972E-4 

19 

(-6. 96445E-1 ,±38.2722) 

5.83594E-4 

20 

(-7.16494E-1,±40.2783) 

8.52026E-5 

The  units  of  the  poles  are  2kb,  where  k  Is  the  wvaenumber 
and  b  Is  the  loop  radius.  {  a  Is  the  wire  radius.) 


Figure  5-21.  A  comparison  of  the  poles  extracted  using  RFA  and  Integral 
equation  and  numerical  search  ("a  "  Impedance  method, 

”o"  RFA  fit)  for  the  loop. 


From  the  backscatterlng  data  of  a  loop  for  broadside  incidence,  an 
Impulse  response  in  the  time  domain  {Figure  5-22)  can  be  found  using  an 
Inverse  Fourier  transform  of  the  data.  Figure  5-23  Is  a  plot  of  the 
Impulse  response  of  the  loop  for  broadside  where  the  specular  return  of 
the  waveform  has  been  removed.  Figure  5-24  Is  a  comparison  of  the 
Impulse  response  waveforms  for  a  loop  and  a  disc.  Obviously,  the  loop 
Is  a  high  Q  scatter  compared  to  the  disc.  Figure  5-25  Is  a  ramp  response 
of  the  loop  using  the  first  10  harmonics  for  Fourier  synthesis.  The 
duration  of  the  response  Is  longer  than  15  TAU  {2b/c,  I.e.,  the  loop 
diameter  transit  time  )  although  the  plot  only  shows  up  to  15  TAU. 

Landt  and  Miller  have  calculated  the  short  pulse  response  of  a  loop  (the 
derivative  of  a  Gaussian  pulse)  using  a  time  domain  approach  (space-time 
Integral  equation).  The  convolution  of  the  impulse  response  obtained 
here  with  a  short  pulse  yields  a  result  very  similar  to  that  of  Landt 
and  Miller's  short  pulse  response  [36].  A  comparison  Is  shown  In  Figure 
5-26. 

A  plot  of  the  Individual  contribution  of  each  pole-pair  and  the 
corresponding  residues  of  the  Impulse  response  for  a  loop  (backscatter) 
at  broadside  Is  shown  In  Figure  5-27.  More  transient  response  waveforms 
for  the  loop  at  broadside  and  edge-on  Incidence  are  shown  In  Appendix  H. 
The  dominant  string  of  poles  of  a  loop  at  edge  on  Incidence  belong  to 
layer  I  (see  Figure  5-17).  Only  one  pole  belonging  to  layer  II  was 
found  using  a  RFA.  Other  poles  were  too  weakly  excited  for  extraction 
using  RFA. 
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Figure  5-22.  The  Impulse  response  (with  the  specular  term  to  be 
removed) for  a  loop  (backscatter)  at  broadside.  The 
specular  return  Is  seen  around  the  time  origin, 
a/b  Is  3.14159E-3,  where  a  Is  the  wire  radius  and  b 
Is  the  loop  radius. 


-o.60  -o.m 


TAU=2b/c 

b  is  the  loop  radius 


Figure  5-23.  The  Impulse  response  (with  the  specular  return  removed) 
for  a  loop  at  broadside  (back scatter),  a/b  Is  3.14159E-3 
where  a  Is  the  wire  radius  and  b  Is  the  loop  radius. 


Figure  5-24.  A  comparison  of  the  Impulse  response  waveforms  for  a 
circular  loop  (a/b  Is  3.14159E-3)  and  a  circular 
metallic  disc  (back scatters). 


Figure  5-26.  The  derivative  of  a  short  Gaussian  pulse  response  for  a 


loop,  (a)  Is  obtained  using  the  convolution  of  the 


Impulse  response  and  the  Input  short  pulse,  (b)  Is 
calculated  using  the  time  domain  moment  method  [36]. 


Note  that  wire  to  the  loop  radius  ratio  Is  0.00314  In  (a) 


TAU“2b/c 

b  is  the  loop  radius 


(c)  third  pair  of  poles 
<-1.259+/-j5.076) 


(d)  fourth  pair  of  poles 
(-l.310V-j8.321) 


Figure  5-27.  The  Individual  contribution  of  each  pole-pair  and  the 
corresponding  residues  to  the  Impulse  response  for  a 
loop  backscatter  for  broadside. 


F.  FINITE  AND  SEMI-INFINITE  OPEN  CIRCULAR  WAVEGUIDES 


The  frequency  domain  axial  backscatterlng  for  the  semi -Infinite* 
circular  waveguide  has  been  studied  and  computed.  In  the  resonance 
region,  using  exact  Wlener-Hopf  solution  by  Johnson  and  Moffatt  [37]. 
Later,  a  UTD  approximation  was  done  by  Huang  [38].  Axial  backscatter 
from  finite  circular  waveguides**  with  the  far  end  shorted  or  opened  have 
been  computed  using  the  UTD  approximation,  and  a  comparison  has  been  made 
to  a  moment  method  solution  [38,  39].  Geometrical  descriptions  of  the 
finite  and  semi-infinite  waveguides  are  shown  In  Figure  5-28. 

Figures  5-29  to  5-34  show  plots  of  RFA  fits  to  portions  of  the 
scattering  data  In  the  frequency  domain  obtained  using  a  moment  method 
solution.  Table  5-14  Is  a  list  of  poles  extracted  from  RFA  and  the  window 
technique  developed  In  Chapter  III  for  these  three  objects.  Figure  5-35 
Is  a  plot  of  the  extracted  poles  for  the  three  geometries.  Note  that  the 
poles  of  the  finite  open  waveguides  are  very  close  for  both  the  far  end 
shorted  or  opened,  and  definitely  more  pole-strings  exist  for  the  the 
finite  waveguides  than  for  the  the  semi -Infinite  guide  because  of 


*  For  the  semi-infinite  circular  waveguide  the  data  are  D/x  of  01 ( .01) .7 
from  Johnson  and  Moffatt' s  exact  solution.  The  data  D/x  of  l.( .01)2.6 
are  the  GTD  high  frequency  approximation.  The  rest  of  the  data  are 
fitted  and  predicted  by  the  rational  function  approximation.  D  Is  the 
diameter  of  the  waveguide. 

**  For  the  finite  circular  waveguide  with  the  far  end  shorted  or  opened 
the  data  D/x  of  ,01(. 01)1. 20  are  from  Huang's  UTD  solution.  D  Is  the 
the  diameter  of  the  waveguide. 


TABLE  5-14 


A  LIST  OF  POLES*  OF  FINITE  AND  SEMI-INFINITE  WAVEGUIDE  OBTAINED 
USING  THE  RFA  TO  THE  BACKSCATTERING  CALCULATED  DATA. 


Finite  Waveguide 

Finite  Waveguide 

Semi -infinite 

The  Rear  End  Open 

Both  Ends  Open 

Waveguide 

-0.22+/-J0.54 

-0.22+/-J0.55 

-0.20+/-J0.82 

-0.21+/-J0.82 

-0.21+/-J1.06 

-0.19+/-J1.04 

-0.19+/-J1.29 

-0.16+/-J1.28 

-0.95+/-J1.36 

-0.15+/-J1.51 

-0.13+/-J1.53 

-0.19+/-J1.63 

-0.17+/-J1.62 

-0.33+/-J1.79 

-0.10+/-J1.81 

-0.11+/-J1.90 

-0.02+/-J2.08 

-0.01+/-J2.05 

-0.05+/-J2.22 

-0.10+/-J2.21 

-0.08+/-J2.42 

-0.17+/-J2.45 

-0.30+/-j2.57 

-0.47+/-J2.57 

-0.10+/-J2.64 

-0.22+/-J2.69 

-0.13+/-J2.88 

-0.31+/-J2.82 

-0.13+/-J3.04 

-0.21+/-J3.04 

-0.11+/-J3.14 

-0.20+/-J3.14 

-0.13+/-J3.26 

-0.19+/-J3.24 

-0.40+/-J3.31 

-0.86+/-J4.95 

-0.80+/-J6.14 

-0.53+/-J7.89 

*  The  poles  are  In  the  units  of  ka,  where  a  Is  the  radius  of  the 
circular  waveguide  and  k  Is  the  wavenumber. 
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OPEN  CIRCULAR  WAVEGUIDE 


SEMI -INFINITE  CIRCULAR  WAVEGUIDE 


Figure  5-28. 


The  geometries  of  finite  waveguides  with  both  ends  open 
(a),  the  rear  end  shorted  (b),  and  semi-infinite 
waveguide  (c)« 
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Figure  5-29.  The  RFA  for  (M,N)  as  (14,13)  (dashed  line)  fit  to  the 
axial  backscatterlng  data  (solid  line)  for  the  finite 
waveguide  (both  ends  open)  within  the  window  of  Interest. 
The  x's  are  data  points  used  In  RFA  (amplitude  plot). 
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Figure  5-30.  The  RFA  for  (M,N)  as  (14,13)  (dashed  line)  fit  to  the 
axial  backscatterlng  data  (solid  line)  for  the  finite 
waveguide  (both  ends  open)  within  the  window  of  Interest. 
The  x's  are  the  data  points  used  In  RFA  (phase  plot). 
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Figure  5-31.  The  RFA  for  (M,N)  as  (10,9)  (dashed  line)  fit  to  the  axial 
back scattering  data  (solid  line)  for  the  finite  circular 
waveguide  (the  rear  end  shorted)  within  the  window  of 
Interest.  The  x's  are  the  data  points  used  In  RFA 
(amplitude  plot). 
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Figure  5-32.  The  RFA  for  (M,N)  as  (10,9)  (dashed  line)  fit  to  the  axial 
backscatterlng  data  (solid  line)  for  the  finite  circular 
waveguide  (the  rear  end  shorted)  within  the  window  of 
Interest.  The  x's  are  the  data  points  used  In  RFA  (phase 
plot). 
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Figure  5-33.  The  RFA  for  (M,N)  as  (16,15)  (dashed  line)  fit  to  the 
axial  back-scattering  data  (solid  line)  for  a 
semi-infinite  circular  waveguide.  The  x's  are  the  data 
points  used  In  the  RFA.  D  Is  the  diameter  of  the  clrcul 
waveguide  (amplitude  plot). 
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Figure  5-34.  The  RFA  for  (M,N)  as  (16,15)  (dashed  line)  fit  to  the 
axial  backscatterlng  data  (solid  line)  for  a 
semi-infinite  circular  waveguide.  The  x's  are  the  data 
points  used  In  the  RFA.  D  Is  the  diameter  of  the 
circular  waveguide  (phase  plot). 
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Figure  5-35.  The  complex  natural  resonances  for  the  three  cavities 
a:  finite  waveguide  (both  ends  opened) 
finite  waveguide  (the  rear  end  opened) 
t:  semi -Infinite  waveguide. 


basically  dipole-type  inodes.  This  can  be  seen  either  In  the  plot  of 
poles  In  Figure  5-35  or  In  Table  5-14. 

The  complex  natural  resonances  of  finite  waveguides  and 
semi -Infinite  waveguide  are  compared  Figure  5-35.  The  complex  natural 
resonances  are  very  close  for  the  finite  waveguides  with  the  far  end 
opened  or  closed  (as  least  In  the  Interesting  spectrum  used  for  the 
RFA).  The  major  resonant  modes  (dipole  modes)  of  the  finite  guides  are 
nearly  the  same  as  can  be  seen  In  either  Figure  5-35  or  Table  5-14. 

The  edge  diffraction  and  rim  creeping  modes  are  dominant  for  the 
semi -Infinite  waveguide  and  the  dipole  modes  are  dominant  for  the  finite 
waveguide.  However,  the  scattering  mechanisms  due  to  the  front  rim  for 
any  of  the  open  waveguides  are  the  same  for  all  three  cavities. 

The  window  technique  was  used  with  the  RFA  to  the  calculated  data 
for  the  finite  waveguides,  therefore  the  poles  due  to  the  resonant  modes 
from  the  front  rim  were  missed  for  both  finite  waveguides.  If  the  rim 
resonances  were  to  be  extracted,  a  very  wide  window  could  be  needed. 

G.  AIRCRAFT  F104 

Table  5-15  shows  calculated  data  for  the  F104  aircraft.  (A  stick 
model  obtained  at  the  ElectroSclence  Laboratory  by  E.  Lin  [40].)  The 
simplest  wire-grid  airplane  model  utilizing  8  wire  segments  was  used  In 
the  claculatlon.  The  stick  model  shown  In  Figure  5-36  Is  not  the  wire 
grid  model  used  by  E.  Lin;  It  only  shows  the  orientation.  This  set  of 
low  frequency  backscattered  data  at  e  of  0°  and  $  of  0*  for  horizontal 
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TABLE  5-15 

THE  CALCULATED  DATA  OF  F104  AIRCRAFT  MODEL* 
AT  1  to  12  MHz 


"Aspect"  9=0°  ,  <j>=0°,  polarization 


Frequency (MHZ) 
1 
2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


Normalized  Echo  Signal 
. 010591-jO. 001060 
. 044214- jO. 009941 
.107527-jO. 025901 
. 186288- j. 061203 
.281474-j. 119757 
. 390186- j. 208745 
. 508891- j. 337996 
. 632569- j. 522298 
. 751477- j. 783962 
.843199- Jl. 163223 
.842512- jl. 698540 
. 612234- j2. 427848 


The  fuselage  length  Is  16.69  m  and  the  wing  length  Is  6.68 


U)  polarization  as  fitted  by  the  rational  functional  approxlmant  Is 
shown  In  Figures  5-37  and  5-38. 

An  extracted  pole  Is  listed  In  Table  5-16  and  Is  compared  with  that 
found  by  Moffatt  and  Chuang  [41]  using  Prony's  method  for  the  same  data. 
The  reason  that  only  one  pole-pair  was  extracted  is  that  the  data  points 
available  were  limited  to  the  Rayleigh  region.  Basically,  the  wire 
model  used  for  the  aircraft  Is  really  a  very  simple  one  In  the 
calculations  for  these  data.  This  example  shows  that  complex  natural 
resonances  can  sometimes  still  be  extracted  even  If  the  imaginary  parts 
of  the  poles  are  outside  the  spectrum  of  the  available  data. 
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TABLE  5-16 


THE  DOMINANT  POLES*  OF  F104  EXTRACTED  VIA 
THE  RATIONAL  FUNCTION  APPROXIMANT  AND  PRONY'S  METHOD 


Rational  Function  Prony's  Method 

(Moffatt  A  Chuang) 

-0.475+/-J2.526  -0.545+/-J2.62 


*  The  units  of  the  poles  are  ka,  where  k  Is  the  wavenumber  and  a  Is 
one  half  length  of  the  fuselage. 


AMPLITUDE 


FREQUENCY  (  ko  ) 


Figure  5-37.  The  RFA  for  (M,N)  as  {12, II)  fit  to  the  calculated 

back scattering  data  (solid  line)  for  the  F104  aircraft 
at  9  of  0°,  of  0°  Incidence  and  polarization.  The 
x*s  are  the  data  points  used  In  the  RFA  (amplitude  plot) 
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CHAPTER  VI 


EXTRACTING  COMPLEX  NATURAL  RESONANCES 
FROM  MEASURED  DATA 


A.  INTRODUCTION 

The  primary  goal  In  tills  chapter  Is  to  describe  In  some  detail  a 
relatively  general  method  for  extracting  complex  natural  resonances  by 
applying  the  RFA  to  measured  scattering  data*.  Because  of  the 
complexity  of  the  target  geometry  for  objects  such  as  aircraft  or  ships. 
It  Is  difficult  (If  not  Impossible)  to  formulate  and  solve  the 
scattering  problem  even  approximately  using  analytical  methods.  Thus 
measured  data  on  some  complex  targets  are  necessary,  particularly  If 
broadband  results  are  needed.  In  this  chapter,  measured  data  of  two 
commercial  passenger  aircraft  scale  models  are  used.  Sufficient  detail 


*  The  scattering  data  on  complex  targets  were  obtained  from  measurements 
being  made  on  a  noncooperative  target  recognition  (NCTR)  program  In 
progress  at  the  ElectroSclence  Laboratory.  The  targets  are  primarily 
fighter  or  other  aircraft.  Since  no  existing  complex  natural 
resonances  are  available  for  comparison  for  these  targets,  we  shall 
only  Identify  them  by  type,  e.g.,  aircraft  A,  B  etc.  The  analysis 
methods  are  adequately  explained  In  applications  to  simple  targets. 
This  approach  will  circumvent  any  possible  future  classification 
problems  while  In  no  way  diminishing  the  application  results. 


of  RFA  applications  have  been  given  on  the  simple  targets,  but  It  Is 
Important  that  the  techniques  developed  herein  be  applied  to  measured 
data  on  real  complex  (geometry)  targets.  At  the  same  time, 
classification  problems  are  avoided  by  not  specifically  Identifying  the 
targets. 

B.  AMPLITUDE  AND  PHASE  PLOTS 

An  advanced  state-of-the-art  broadband  radar  reflectivity 
measurement  facility  Is  under  development  at  the  ElectroSclence 
Laboratory  of  the  Ohio  State  University  [42].  Detailed  procedures  for 
the  measurement  and  data  processing  are  not  discussed  here  but  It  Is 
helpful  to  present  some  of  the  scattering  data  as  a  function  of  frequency 
for  the  scale  model  targets.  The  coordinate  system  and  the  target 
orientation  (aircraft  shown  In  stick  form  only)  which  were  used  In  the 
measurements  are  shown  In  Figure  6-1.  The  aircraft  Is  aligned  with  wings 
and  fuselage  In  the  xy  plane  (nose  In  the  positive  x  with  direction  and 
vertical  stabilizer  In  the  z  direction).  Vertical  polarization  Is 
Identified  as  measurements  with  the  electric  field  polarized  in  the  z 
direction  and  horizontal  polarization  as  measurements  with  the  electric 
field  polarized  In  the  xy  plane. 

For  horizontal  polarization,  the  complex  data  were  taken  In  the 
frequency  band  of  1.5  GHZ  to  6.5  GHZ  In  500  steps  for  aircraft  A.  Also 
measured  data  In  1-2  GHz,  2-4  GHz,  and  4-6  GHz  bands  (200  steps  In  each 
band)  were  taken  for  vertical  polarization  for  aircraft  A  and  1.0  to  6.5 


GHz  (550  steps)  for  aircraft  B.  Figures  6-2  and  6-3  (amplitude  plots) 
and  Figures  6-4  and  6-5  (phase  plots)  are  the  spectrum  at  aspect  angles 
of  0°,  20°,  30°,  45°,  50°  off  the  nose  for  horizontal  polarization  for 
aircraft  A.  Figures  6-6  and  6-7  are  amplitude  and  phase  plots  for 
vertical  polarization  at  aspect  angles  of  0°,  15°,  and  30°  off  nose-on 
for  aircraft  A,  and  Figures  6-8  to  Figure  6-11  are  plots  (Figure  6-8  and 
6-9  amplitude,  Figure  6-10  and  6-11  phase)  for  vertical  polarization  at 
0®,  30°,  60°,  90°,  135°,  180°  off  nose-on  for  aircraft  B.  The  amplitude 
is  /o' and  is  in  the  units  of  cm,  where  a  Is  the  radar  cross  section  (RCS) 
on  a  linear  scale.  (The  RCS  is  In  the  units  of  cm2.) 

The  plots  of  the  spectra  are  shown  here  because  the  number  of  peaks 
In  the  amplitude  plots  are  clues  for  determining  the  system  order. 
Furthermore,  as  will  be  seen  later,  the  Imaginary  part  of  the  extracted 
poles  are  near  most  of  those  peaks.  A  comparison  of  the  amplitude  plots 
of  the  measured  data  for  aircraft  A  at  horizontal  and  vertical 
polarizations  from  nose-on  Is  shown  In  Figure  6-12.  All  of  these  plots 
shown  In  this  Chapter  have  been  smoothed  by  the  zero-phase-shift  digital 
filter  given  In  Appendix  B. 

It  Is  difficult  to  smooth  noisy  data  precisely  (without  any 
distortion  of  the  signal)  because  the  transient  responses  vary  in  length 
for  various  targets,  aspects,  and  polarizations.  However,  strongly 
excited  poles  will  not  be  Influenced  much  as  long  as  the 
cut-off  frequency  of  the  zero-phase  low-pass  digital  filter  Is  properly 
designed. 
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Figure  6-2.  The  backscattered  amplitude  for  commercial  transportation 
aircraft  A  at  0°,  10°,  and  20°  from  nose-on  for  horizontal 
polarization. 
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Figure  6-3.  The  backscattered  amplitude  for  commercial  transportation 
aircraft  A  at  30  ,  45°,  and  50°  from  nose-on  for  horizontal 
polarization. 
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Figure  6-4.  The  back scattered  phase  for  commercial  transportation 

aircraft  A  at  0°,  10°,  and  20°  from  nose-on  for  horizontal 
polarization. 
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Figure  6-5.  The  backscattered  phase  for  commercial  transportation 

aircraft  A  at  30°,  45°,  and  50°  from  nose-on  for  horizontal 
polarization. 
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Figure  6-8.  The  backscattered  amplitude  for  the  commercial 

transportation  aircraft  B  at  0°,  30°,  and  60°  from  nose-on 
for  vertical  polarization. 
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Normally  4  TAU  width  for  the  time-domain  window  (TAU  Is  2a/c,  where 
a  Is  a  characteristic  target  length)  Is  a  good  choice  for  most  targets. 
Some  targets  however  should  have  windows  (the  3  dB  cut-off  frequency  of 
the  digital  low-pass  filter)  wider  than  4  TAU.  Such  targets  have 
resonances  with  small  damping  such  as,  e.g.,  the  circular  loop.  The 
cut-off  frequency  of  the  filter  In  the  plots  of  Figures  6-2  to  6-8  Is 
set  to  23  Hz;  and  the  sampling  rate  of  the  digital  filter  Is  assumed  to 
be  500  Hz.  The  folding  frequency  Is  250  Hz.  The  corresponding 
time-domain  window  width  Is  therefore  roughly  4  TAU  for  this  cutoff 
frequency  of  23  Hz.  The  sampling  Intervals  In  the  spectra  are  10  MHz. 
Typical  examples  of  oversmoothing  using  zero-phase-shift  digital  filters 
are  shown  In  Figures  6-13  to  6-16.  The  3dB  cut-off  frequency  Is  set  to 
fc  of  11.5  Hz,  l.e.,  2  TAU  lengths  In  the  time  domain.  It  has  been 
found  that  some  of  the  true  poles  were  missed  If  oversmoothed  data  are 
used  In  the  rational  function  approximation.  Typical  plots  of  the 
oversmoothed  data  for  aircraft  A  are  shown  In  Figures  6-13  and  6-14 
(horizontal  polarization),  and  In  Figures  6-15  and  6-16  (vertical 
polarization)  at  nose  on  Incidence.  Obviously,  these  plots  show  that 
higher  order  scattering  was  filtered  out.  Oversmoothing  results  In 
loss  of  some  weak  resonances.  To  filter  out  noise  and  clutter,  plots  of 
the  Impulse  response  are  necessary.  The  Impulse  response  helps  select 
fc,  the  3  dB  cutoff  frequency,  of  the  digital  filter  as  will  be  shown  In 
section  D  of  this  chapter. 
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E  Figure  6-13.  An  amplitude  plot  of  the  backscatterlng  data  for 

commercial  transportation  aircraft  A  at  nose-on. 
The  measured  data  are  oversmoothed  by  a 
t’  zero-phase-shift  digital  filter. 
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Figure  6-14.  A  phase  plot  of  the  hack scattering  data  for  commercial 
transportation  aircraft  A  at  nose-on.  The  measured 
data  are  oversmoothed  by  a  zero-phase-shift  digital 
filter. 
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i— 15*  An  Amplitude  plot  of  the  back scattering  data  for 

commercial  transportation  aircraft  A  at  nose-on  for 
vertical  polarization.  The  data  are  oversmoothed  by  the 
zero-phase-shift  digital  filter. 


Figure  6-16.  A  phase  plot  of  the  backscatterlng  data  for  commercial 
transportation  aircraft  A  at  nose-on  for  vertical 
polarization.  The  data  are  oversmoothed  by  the 
zero-phase-shift  digital  filter. 


C.  EXTRACTING  THE  COMPLEX  NATURAL  RESONANCE  FROM  THE  MEASURED  DATA 


After  RFA  have  been  developed,  It  Is  relatively  routine  to  apply 
the  programs  to  measured  data.  The  following  features  should 
be  noted: 

1.  The  measured  data  sets  are  smoothed  by  a  low  pass  10th  order 
zero  phase  digital  filter.  The  3dB  cutoff  frequency  of  the  low 
pass  digital  filter  Is  movable. 

2.  The  window  technique  Is  used  with  the  RFA.  This  window  Is 
generated  simply  by  taking  sample  points  over  a  band  limited 
region  of  the  spectrum.  The  window  used  for  the  RFA  Is  a  very 
powerful  tool  for  extracting  the  poles,  because  there  are 
usually  too  many  resonances  within  the  spectrum  to  be  covered 
by  one  RFA. 

3.  The  system  order  Is  set  at  around  twice  the  number  of 
amplitude  peaks  In  the  spectrum  covered  by  the  window. 

4.  Several  different  sets  of  data  points  and  system  orders  are 
tested  at  each  window,  and  only  the  poles  which  are  unchanged 
for  more  than  two  selected  data  sets  are  considered  as 
candidate  poles. 

5.  The  window  Is  chosen  such  that  the  magnitude  of  the  fitted 
data  lying  outside  the  window  Is  as  small  as  possible.  Thus, 
both  end  points  of  the  window  should  be  at  relatively  minimum 
amplitudes  or  at  least  on  a  portion  of  the  negative  slope. 


1.  Aircraft  A 


The  complex  natural  resonances  of  the  transportation  aircraft  A  at 
0°,  30* ,  45®,  60®,  90®,  120®,  150“ ,  165®,  and  180®  aspect  angles  as 
extracted  from  the  backscatterlng  data  In  the  1-2  GHz  band  for 
horizontal  polarization  are  given  In  Table  6-1.  The  1-2  GHz  and  2-4  GHz 
bands  are  treated  Individually  because  the  measured  data  sets  for  these 
bands  Mere  taken  separately.  Table  6-2  lists  the  extracted  poles  from 
the  backscattered  data  at  0®,  15®,  30®,  45",  60®,  75®,  105®,  165®,  and 
180®  aspect  angles  for  the  aircraft  A  and  for  the  horizontal 
polarization  In  the  2-4  GHz  band.  HoMever,  the  rational  function 
approximation  windows  are  set  at  2-3  GHz  only.  The  pole  locations  found 
are  relatively  aspect  Independent. 

Table  6-3  shows  the  extracted  poles  for  horizontal  polarization  at 
aspect  angles  of  0®,  10®,  20®,  and  30*  as  found  using  1.5  to  6  GHz  band 
data.  In  order  to  compare  with  the  other  polarization  (vertical 
polarization)  only  data  In  the  1.5  GHz  to  3.0  GHz  band  were  used.  The 
extracted  poles  are  relatively  aspect  Independent.  Note  that  more  poles 
are  excited  for  horizontal  polarization  than  for  vertical  polarization. 

2.  Aircraft  B 

The  complex  natural  resonances  of  aircraft  B  were  extracted  using 
the  RFA  and  window  techniques  for  data  at  0®,  90®,  and  135®  from  nose-on 
In  tiie  1  to  2  GHz  band.  The  extracted  poles  are  listed  In  Table  6-4. 

The  poles  are  relatively  aspect  Independent  and  the  pole  sets  for  the 


TABLE  6-1 


THE  COMPLEX  NATURAL  RESONANCES*  FOR  COMMERCIAL 
AIRCRAFT  A  OBTAINED  FROM  THE  RATIONAL  FUNCTION 
APPROXIMATION  TO  THE  MEASURED  DATA  (BACKSCATTER, 
VERTICAL  POLARIZATION,  1-2  GHz). 


1-2  GHz 


O*. 

-0.082  +/-J1.082 
not  located 
-0.092  +/-J1.318 
-0.783  +/-J1.588 
-0.645  +/-J1.765 
-0.0755+/-J1 .915 

60° 

-0.068+/-J1.077 
-0. 082+/- j 1.187 
-0. 102+/- j 1.34 
-0. 082+ /-j 1.53 
-0.056+/-J1.71 
-0. 09 1+/-J 1.956 

150° 

-0.051  +/-J1.053 
-0.052  +/-J1.18 
-0.092  +/-j 1.351 
-0.09  +/-J1.559 
-0.0936+/-J1.72 
-0.093  +/-j 1.959 


30° 

-0.104+/-J1.034 
-0. 093+/- j 1.16 
-0.116+/-J1.338 
-0.085+/-J1.554 
-0.045+/-J1.73 
-0.05  +/-J2.06 

90° 

-0. 035+/ -J 1.057 
—0.031+/— jl . 171 
-0.11  +/-J1.39 
-0.15  +/-J1.574 
-0.063+/-J1.775 
-0.11  +/-J1.95 

165° 

-0.53  +/-J1.036 
-0. 057+/- J 1.141 
-0. 056+/- j 1.383 
-0.057+/-J1.572 
-0.055+/-J 1.726 
-0.11  +/-jl.94 


45° 

-0.05+/-J1.07 
-0.06+/-J1.18 
-0. 10+/- j 1.33 
-0.07+/-J1.52 
-0.09+/-J1.71 
-0.14+/-J1.95 

120° 

-0. 039+/- j 1.03 
-0.06  +/-J1.11 
-0.10  +/-J1.33 
-0.07  +/-J1.54 
-0.08  +/-J1.77 
-0. 073+/- j 1.934 

180° 

-0. 058+/- j 1.05 
-0.062+/-J1.161 
-0. 084+/- j 1.324 
-0.057+/-J1.512 
-0. 040+/- j 1.749 
-0. 029+/- j 1.905 


*  The  units  of  the  poles  are  In  GHz. 


186 


TABLE  6-2 


THE  COMPLEX  NATURAL  RESONANCES*  FOR  COMMERCIAL 
AIRCRAFT  A  OBTAINED  FROM  THE  RATIONAL  FUNCTION 
APPROXIMATION  TO  THE  MEASURED  DATA  {BACKSCATTER, 
VERTICAL  POLARIZATION,  2-3  GHz). 


2-3  GHz 

Aspect  Angle  0° 

Aspect  Angle  15* 

Aspect  Angle  30' 

-0.02  +/-J1.996 
-0.108+/-J2.328 
-0.07  +/-J2.681 
-0.068+/-J2.974 

-0.146+/-J2.291 

-0.107+/-J2.665 

-0.100+/-J2.967 

-0.008+/-J2.019 
-0.068+/-J2.20 
-0.10  +/-J2.51 
-0.08  +/-J2.72 
-0.192+/-J3.00 

45° 

-0.05  +/-J2.13 
-0.09  +/-J2.37 
-0.027+/-J2.65 
-0.047+/-J2.94 

60° 

-0. 047+/- j 1.94 
-0.091+/-J2.22 
-0.09  +/-J2.49 
-0.12  +/-J2.79 
-0.10  +/-j 2.93 

75° 

-0.053+/-J2.109 
-0.127+/-J2.317 
-0.08  +/-J2.593 
-0.097+/-J2.922 

105° 

-0.051+/-J2.019 
-0.079+/-J2.36 
-0.07  +/-J2.61 
-0.047+/-J2.86 
-0.06  +/-J2.97 

165° 

-0.099+/-J2.089 

-0.128+/-J2.21 

-0.016+/-J2.362 

-0.037+/-J2.62 

-0.032+/-J2.97 

180° 

-0.117+/-J2.16 
-0.105+/-J2.29 
-0.09  +/-J2.47 
-0.09  +/-J2.88 
-0.16  +/-J3.00 

The  units  of  the  poles  are  In  GHz 


TABLE  6-3 


THE  COMPLEX  NATURAL  RESONANCES*  FOR  COMMERCIAL 
AIRCRAFT  A  OBTAINED  FROM  THE  RATIONAL  FUNCTION 
APPROXIMATIONS  TO  THE  MEASURED  DATA  (BACKSCATTER 
HORIZONTAL  POLARIZATION,  1.5-3  GHz). 


0* 

10° 

20° 

30“ 

(-0.12,1.55) 

(-0.13,1.58) 

(-0.057,1.513) 

(-0 

(-0.14,1.69) 

(-0.12,1.69) 

(-0.12,1.67) 

(-0 

(-0.11,1.75) 

(-0.12,1.73) 

(-0.11,1.79) 

(-0 

(-0.10,1.81) 

(-0.13,1.801) 

(-0.14,1.85) 

(-0 

(-0.08,1.93) 

(-0.12,1.95) 

(-0.10,1.93) 

(-0 

(-0.24,2.00) 

(-0.14,2.02) 

(-0.10,2.06) 

(-0 

(-0.08,2.16) 

(-0.09,2.18) 

(-0.07,2.17) 

(-0 

(-0.15,2.22) 

(-0.10,2.27) 

(-0.57,2.21) 

(-0 

(-0.10,2.34) 

(-0.13,2.34) 

(-0.11,2.36) 

(-0 

(-0.07,2.57) 

(-0.16,2.52) 

(-0.05,2.54) 

(-0 

(-0.09,2.63) 

(-0.23,2.62) 

(-0.19,2,67) 

(-0 

(-0.10,2.76) 

(-0.104,2.80) 

(-0.57,2.86) 

(-0 

(-0.14,2.99) 

(-0.13,2.98) 

(-0.14,2.94) 

(-0 

08,1.52) 

056,1.65) 

086,1.78) 

095,1.88) 

052,1.93) 

046,2.03) 

098,2.15) 

063,2.23) 

16,2.37) 

123,2.55) 

03,2.66) 

06,2.79) 

03,3.04) 


The  poles  are  In  the  units  of  GHz 


TABLE  6-4 


THE  COMPLEX  NATURAL  RESONANCES*  OF  AIRCRAFT  B  AT 
0“,  90° ,  and  135°  FROM  NOSE -ON  INCIDENCE 
FOR  VERTICAL  POLARIZATION 


0° 

-0. 073+/- j 1.095 
-0. 19+/- j 1.232 
-0. 24+/- j 1.450 
-0.48+/- jl. 628 
-0.21+/-J1.956 


90° 

— 0.21+/— Jl .097 
—0 . 26+/- Jl . 25 
-0.19+/-J1.47 
-0.19+/-J1.66 


135° 

-0.165+/-J1.023 

-0.183+/-J1.214 

-0.129+/-J1.437 

-0.51+/-J1.65 

-0.07+/-J1.88 

-0.18+/-J1.953 


*  The  poles  are  In  the  units  of  GHz.  The  RFA  fit  In  the  1-2  GHz 


two  aircraft  are  not  the  same.  Therefore,  It  would  appear  that  the  pole 
sets  could  be  used  for  the  purpose  of  classification. 

One  of  the  characteristics  of  the  complex  natural  resonances, 
polarization  Invariance,  can  be  tested  by  comparing  Table  6-1  to  6-3. 

It  Is  well-known  that  the  complete  set  of  complex  natural  resonances  of 
an  object  may  not  be  optimally  excited  for  one  polarization  or  for  one 
aspect  angle.  For  example,  a  wire  for  broadside  Incidence  has  only  even 
modes  excited;  at  other  aspects  both  even  and  odd  modes  are  excited. 

But  the  total  resonances  of  the  wire  Include  both  even  and  odd 
resonances.  Similarly,  If  the  target  has  poles  which  are  more  strongly 
excited  with  horizontal  polarization  than  with  vertical  polarization 
then  more  poles  are  extracted  from  the  measured  aircraft  data  for 
horizontal  polarization.  Portions  of  the  extracted  poles  are  nearly  the 
same  for  both  polarizations. 

0.  IMPULSE  RESPONSES  OF  THE  AIRCRAFT 

Our  purpose  here  Is  to  plot  some  of  the  band  limited  Impulse 
response  waveforms  of  the  aircraft  using  Fourier  synthesis.  We  use  an 
Inverse  Fast  Fourier  Transform  (IFFT)  to  Implement  the  Fourier  synthesis 
and  a  smoothing  window  is  used  to  eliminate  the  Gibbs'  phenomenon.  The 
window  (W(n))  used  here  Is 


pi* 

J 


190 


«1n2j!>  _  (»)! 


— 2 - for  n  not  equal  N/2 

r  (“)  -  n2 

for  n  of  N/2 


,  (6-1) 


where  N  Is  the  total  number  of  samples.  Study  of  the  transient  response 
waveforms  for  the  aircraft  Is  not  for  the  purpose  of  target 
Identification  but  rather  for  deciding  the  3  dB  cutoff  frequency  for  the 
digital  filter.  Then,  the  noise  and  clutter  can  be  reduced  from  the 
frequency  domain  measured  data  using  a  digital  filter  with  proper  cutoff 
frequency.  Typical  Impulse  responses  for  a  transportation  aircraft  at 
some  different  aspect  angles  and  for  both  polarizations  are  shown  at  the 
end  of  this  section.  The  rough  geometry  of  the  aircraft  Is  sketched 
Indicating  aircraft  orientation  and  Incident  wave  polarization  for  each 
Impulse  response  waveform. 


E.  BRIEF  SUMMARY  AND  DISCUSSION 


In  this  dissertation  no  effort  Is  made  to  relate  features  of  the 
Impulse  response  waveforms  shown  to  geometrical  features  of  the 
aircraft.  To  do  so  might  compromise  the  actual  Identity  of  the  targets 
and  Is  not  essential  for  any  purpose.  Both  aircraft  show  response 
waveforms  which  ring  for  relatively  long  times.  Indicating  that  signals 
are  bouncing  between  various  structures.  Given  the  actual  size  and 
shape  of  the  aircraft,  most  of  the  dominant  scattering  mechanisms 
can  be  Identified.  For  our  purpose,  It  Is  also  clear  that  an 
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Figure  6- 
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.  Backscattered  Impulse  response  waveform  of  commercial 
transportation  aircraft  A  at  20°  from  nose-on  for 
horizontal  polarization. 
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Figure  6-25.  Back scattered  Impulse  response  waveform  of  commercial 
transportation  aircraft  A  at  90°  from  nose-on  for 
horizontal  polarization. 


E5P0NSE (VOLT) 

'0.00  300.00  600.00 


IMPULSE  RESPONSE  (VOLT) 

-600.00  -300.00  'U.OU  300.00  600.00 


examination  of  the  response  waveforms  Is  essential  for  selection  of  the 
best  digital  filter. 

The  ultimate  goal  In  the  chapter  -  extracting  the  complex  natural 
resonances  from  the  measured  complex  shape  targets  -  has  been 
accomplished  for  two  passenger  transportation  aircraft.  The 
applications  of  the  RFA  were  limited  to  1  to  3  GHz  in  order  that  the 
poles  for  many  aspect  angles  and  two  polarizations  could  be  extracted 
and  compared  easily.  The  number  of  complex  natural  resonances  extracted 
Is  somewhat  Influenced  by  the  3  dB  cut-off  frequency  of  the  digital 
filter.  It  Is  for  this  reason  that  the  Impulse  response  waveforms  for 
the  two  aircraft  were  synthesized.  The  Impulse  response  waveform  with 
zero-phase-shift  digital  filter  Is  shown  In  Figure  6-32.  The  horn 
antenna  coupling,  long  path  multiple  scattering,  target  return 
reflector,  system  noise,  and  back  wall  are  also  displayed  and  clearly 
seen.  The  result  of  filtering  the  measured  data  using  a  zero  phase 
low-pass  digital  filter  (fc  of  15  Hz)  Is  shown  In  Figure  6-33.  One  of 
the  advantages  of  using  the  digital  filter  directly  In  the  frequency 
domain  Is  an  Inverse  Fourier  transform  Is  not  required.  Therefore 
possible  distortion  due  to  the  transformation  can  be  avoided. 

The  exact  3  dB  cut-off  frequency  of  the  lowpass  filter  should  be 
different  for  various  targets,  aspect  angles,  and  polarizations.  The 
Impulse  responses  for  the  aircraft  are  valuable  In  the  determination  of 
the  3  dB  lowpass  filter. 

The  RFA  and  the  preprocessing  techniques  developed  In  this  study 
can  be  applied  to  obtain  the  complex  natural  resonances  of  a  complex 
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shape  using  measured  data.  The  complex  natural  resonances  obtained 
from  the  measured  data  are  relatively  aspect  and  polarization 
Independent. 
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Figure  6-32.  The  Impulse  response  of  aircraft  A  at  nose-on 

for  horizontal  polarization.  The  whole  range  Is 
shown  here. 
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Figure  6-33.  The  Impulse  response  of  aircraft  A  at  0°  aspect  angle 

and  horizontal  polarization.  Note  the  noise  and  clutter 
were  filtered  out  by  a  lowpass  zero-phase-shift  digital 
filter  (fc  of  15  Hz). 


CHAPTER  VII 


CONCLUSIONS 


An  analytical  model  In  the  form  of  a  rational  function 
approximation  (RFA)  for  fitting  available  complex  spectral  data  has  been 
developed.  The  available  data  which  may  be  measured  or  calculated 
correspond.  In  principle,  to  any  portion  of  the  spectrum.  The  Interest 
In  this  study  Is  In  scattered  electromagnetic  signals  and  In  the 
extraction  of  complex  natural  oscillations  (poles)  from  the  fitted 
model.*  The  methods  developed  could  have  much  wider 
applicability— to  radiated  as  well  as  to  scattered  signals,  for 
example,  and  to  seismic  and  acoustic  as  well  as  to  electromagnetic  data. 

Cther  methods  for  extracting  complex  natural  resonances  from 
data  have  been  briefly  reviewed,  and  the  differences  and  advantages  of 
the  present  technique  explored.  It  Is  also  demonstrated  that  rational 
function  models  can  be  used  to  supplement  and  complete  an  analytical 


*  The  rational  function  approximation  developed  here  Is  also  an 
excellent  tool  for  complex  data  Interpolation.  The  application  of 
the  RFA  to  data  Interpolation  has  been  made  for  thin  wires  and 
thick  wires.  The  results  have  been  shown  to  be  very  accurate  when 
compared  with  the  calculated  data.  A  typical  example  can  be  seen  In 
Figure  5-12,  where  the  first  resonances  Is  not  shown  In  the  original 
data,  but  Is  plotted  In  the  RFA  (dotted  line). 


model  (the  K-pulse  [6])  useful  for  relatively  simple  objects  at  high 
frequencies. 

The  RFA  developed  here  has  three  main  features:  low-pass  filtering 
to  combat  noise  and  clutter  effects,  windowing  procedures  to  eliminate 
so-called  curve  fitting  poles,  and  a  sum  operator  which  permits  most  of 
the  available  data  to  be  used  and  Is  usually  more  effective  than  least 
squared  error  techniques  (at  least  against  an  additive  Gaussian  white 
noise).  There  Is,  admittedly,  little  original  In  the  generation  of  a 
system  of  linear,  simultaneous  equations  for  fitting  complex,  multiple 
frequency  data.  Yet  the  method  of  solution  of  the  main  data  points  and 
In  the  preprocessing  of  the  data  (fitting,  sum  operator)  before  pole 
extraction  are  new.  Although  the  final  step,  l.e.,  pole  extraction  via 
of  the  zeros  of  the  denominator  of  the  rational  function  approxlmant  Is 
routine.  It  requires  using  double  precision  on  a  digital  computer  to 
Illuminate  near  pole-zero  cancellations. 

Application  of  the  RFA  has  been  demonstrated  for  a  variety  of 
simple  and  complex  objects,  and,  where  possible,  the  results  (complex 
natural  resonances)  have  been  compared  to  those  obtained  by  other 
methods.  Included  here  are  spheres,  discs,  loops,  thin  wires,  and 
finite  and  semi -Infinite  circular  waveguides.  Numerous  tests  of 
applications  to  a  conducting  spherical  scatterer  have  been  made.  It 
must  be  stressed,  however,  that  these  methods  have  been  tested  against 
true  poles  with  true  residues  (but  finite  In  number)  In  an  additive 
Gaussian  white  noise  model.  Tests  of  a  parameter  extraction  procedure 
where  equal  excitation  of  parameters  Is  assumed,  albeit  In  noise,  are 


relatively  meaningless.  The  RFA  has  also  been  applied  to  measured 
wl de-band  scattering  data  obtained  using  a  new  reflectivity  facility 
being  developed  at  the  ElectroSclence  Laboratory  [42].  The  targets  In 
this  case  are  of  complex  geometry  and  of  sufficient  tactical  Importance 
to  warrant  classification  (security)  problems.  Stressed,  therefore.  Is 
the  stability  of  results  of  RFA,  and  neglected  Is  the  precise 
Identification  of  the  scatterers  for  these  cases.  Excellent  pole 
extraction  has  been  demonstrated  for  both  simple  and  complex  geometries. 
Certain  far-reaching  conclusions  can  be  drawn  from  the  results  of  the 
research  reported  In  this  study. 

1.  Complex  natural  resonances  can  be  successfully  extracted  from 
additive  noise  contaminated  data  with  si gnal -to-nol se  ratios  as 
low  as  15.0  dB  or  lower. 

2.  For  si gnal -to-nol se  ratios  significantly  lower  than  15  dB,  It 
does  not  appear  that  reliable  pole  extraction  methods  are 
feasible.  Unless,  of  course,  some  new  breakthrough  Is  made. 

3.  From  2  above  It  has  been  concluded  that  pole  extraction 
procedures  from  full  scale,  field-measured  data  are  generally 
not  feasible  and  should  not  be  attempted.  Roughly, 

si gnal -to-nol se  ratios  of  13  to  15  <ffi  are  necessary  for 
successful  pole  extraction.  This  does  mean  that  pole 
extraction  from  controlled  data,  l.e.,  from  laboratory 
measurements  Is  feasible.  The  new  ElectroSclence  Laboratory 
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reflectivity  measurement  facility  or  other  careful 
measurements  easily  satisfy  this  requirement. 

4.  While  this  study  does  not  explore  target  recognition.  It  must 
be  stressed  that  pole  extraction  from  real  time  radar  data 

Is  not  necessary  for  target  recognition  using 
prediction-correlation  processing.  The  nearest  neighbor 
decision  rule  [43]  using  the  complex  natural  resonances  can 
also  be  used.  The  rule  Is  to  Identify  a  set  of  extracted 
poles  (Pf>,  as  a  member  of  the  pole  set  of  the  known  targets 
In  a  library,  Cr,  to  which  Its  nearest  neighbor  belongs. 
Therefore,  noise  and  clutter  problems  do  not  negate  this 
form  of  target  recognition  at  least  as  far  as  pole 
extraction  Is  concerned. 

5.  Parameter  extraction  for  the  complex  natural  resonances  Is 
not,  at  this  stage,  an  exact  science  using  the  procedures 
of  this  study.  Experience,  In  processing  and  In  realizing 
what  Is  and  what  Is  not  reasonable  In  terms  of  electromagnetic 
reverberations,  Is  Indispensable.  A  completely  foolproof 
automatic  procedure  cannot  be  devised  at  this  time. 

Using  both  RFA  and  asymptotic  estimates,  new  and  significant 
results  have  been  obtained  for  the  complex  natural  resonances  of  the 
thin  circular  loop.  Of  particular  Importance  Is  the  fact  that  poles 
from  RFA  which  are  most  precise  at  lower  frequencies  show  a  region  of 
overlap  with  poles  obtained  from  asymptotic  methods.  Therefore,  from  a 
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combination  of  the  two  methods  relatively  lengthy  {>30  or  40)  pole-pair 
strings  can  be  obtained.  Such  long  pole  strings  are  precisely  what  are 
needed  to  obtain  a  K-pulse  waveform  using  one  particular  approach  [44]. 

A  thln-wlre  circular  loop  has  also  been  used  to  test  RFA  for 
complex  natural  resonances  extraction.  In  this  case,  comparisons  are 
made  with  existing  poles  obtained  from  an  Integral  equation  formulation 
and  numerical  search;  good  agreement  Is  obtained.  A  K-pulse  for  the 
circular  loop  using  just  those  poles  corresponding  to  creeping  wave  type 
modes  Is  obtained,  and  the  backscatter  response  of  the  loop  to  this 
K-pulse  for  edge-on  Incidence  has  been  derived.  Derivation  of 
K-pul se-response  waveform  pairs  for  a  particular  target  Is  beyond  the 
scope  of  this  study  but  these  Initial  results  do  show  the  utility  of 
RFA.  In  the  course  of  this  study,  new  canonical  response  waveforms 
have  been  shown  for  the  thln-wlre  loop  and  for  cross-polarized  responses 
of  the  circular  disc. 

As  noted  above,  RFA  has  been  used  to  extracted  complex  natural 
resonances  of  several  relatively  complex  target  geometries.  The 
broadband  scattering  data  for  these  applications  were  measured  on  the 
new  compact  reflectivity  measurement  facility  at  the  ElectroSclence 
Laboratory.  In  this  dissertation  a  number  of  band-limited  Impulse 
response  waveforms,  generated  from  these  measured  data,  have  been  shown. 
These  waveforms  vividly  demonstrate  the  tremendous  diagnostic  potential 
of  time  domain  response  waveforms.  The  response  from  various 
geometrical  features  of  the  targets  as  well  as  possible  reverberations 
between  various  features  are  easily  seen.  This  aspect  of  the  response 


waveforms  has  not  been  stressed  here  for  obvious  reasons  but  one  must 
perceive  the  future  significance  of  what  might  be  termed  "scattering 
reflectometry." 

The  most  Immediate  application  of  the  tools  developed  In  this 
dissertation  is  to  extract  complex  target  poles  for  target  recognition 
using  prediction-correlation  [45].  In  this  regard  It  is  noted  that  the 
poles  of  complex  target  geometries  as  extracted  using  RFA  do  not  show 
precise  excitation  Invariance,  particularly  with  respect  to  the  damping 
of  a  given  pole.  The  reason— weak  excitation  of  a  particular 
resonances  at  a  given  aspect  or  polarization— Is  obvious.  As  a 
practical  matter,  however.  It  Is  clear  that  a  single  complex  geometry 
target  will  have  to  be  treated  as  several  targets  (somewhat  different 
pole  damping)  In  prediction-correlation  processing.  Research  on  this 
aspect  of  the  problem  using  the  RFA  tools  developed  here  has  already 
been  Initiated. 
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APPENDIX  A 


A  ZERO  PHASE  SHIFT  LOWPASS  DIGITAL  FILTER 


The  goal  In  this  Appendix  Is  to  design  a  10th  order  zero-phase-shift 
lowpass  digital  filter.  The  first  step  here  Is  to  design  a  5th  order 
lowpass  Butterworth  digital  filter  in  terms  of  a  5th  order  analog  filter 
[46].  Then,  two  5th  order  filters,  H(z)  and  H*(z)  (H*(z)  is  the  complex 
conjugate  of  H(z))  are  cascaded  to  make  a  10th  order  zero-phase-shift 
digital  filter.  A  block  diagram  of  5th  order  digital  filter  Is  shown  In 
Figure  A-l.  The  transfer  function  of  the  5th  order  Butterworth  analog 
filter  [46]  with  a  3  dB  cutoff  frequency  at  1  radian/second  Is  given  by 


G(S)  - _ ^ 

BQ  +  BjS  +  B2s2  +  B3s^  +  B4s4  +  Bgsb 


(A-l) 


where 


Aq  =  B0  =  1  , 

(A-2.1) 

Bi  »  3.236068 

(A-2.2) 

B2  *  5.236068 

(A-2.3) 

B3  *  5.236068 

(A-2.4) 

B4  *  3.236068 

(A-2.5) 

*-  •-  "V 


Y. 


r-.. 

V- 

L 


and 

85  »  1 

If  the  transfer  function  of  the  digital  filter  is  written  as 


H(z) 


-1  -2  -3  -4  -5 

a0  +  aiz  1  +  a2z  +  a3z  +  +  A5Z 


1  +  bjZ"1  +  b2z'z  +  b3z-3  +  b4z-*  +  a5z_s 

the  following  relationships  are  true  [45]. 

A  «  B0  +  BiC  +  B2C2  +  B3C3  +  B4C4  +B5C5  , 

®o  ~ 

ai  *  5A0/A 
a2  3  IOAq/A 

33  3  10Ao/A 

34  *  5A0/A 
a5  *  Aq/A 

bi  3  (5B0  +  3B]C  +  B2C2  -  B3C3  -  3B4C4  -  5B5C5)/A  , 

3  (10Bo  +  2BiC  +  2B2C2  -  2B3C3  +  2B4C4  +  10B5C5)/A 
t>3  3  (10Bo  +  2BiC  +  2B2C2  +  2B3C3  +  2B4C4  -  10B5C5)/A 
b4  =  (5B0  -  3BjC  +  B2C2  +  B3C3  -  3B4C4  +  5B5C5)/A  , 
b5  3  (B0  -  BjC  +  B2C2  -  B3C3  +  B4C4  -  B5C5)/A  , 


and 


C  3  x  cot(  "  l£  ) 
r  7TZ 


The  remaining  undefined  parameters  are: 
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(A-2.6) 


(A-3) 


(A-4.1) 

(A-4.2) 

(A-4.3) 

(A-4.4) 

(A-4.5) 

(A-4.6) 

(A-4.7) 

(A-4.8) 

(A-4.9) 

(A-4.10) 

(A-4.11) 

(A-4.12) 

(A-4.13) 


r.  *  -  <*.  ** 
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xr  Is  the  3  dB  cutoff  frequency  for  the  analog  filter  (assuming 
that  xr  Is  1  rad/sec  for  the  analog  filter  given  by  Equation  (A-l)),  f0 
Is  the  folding  frequency,  l.e.,  one  half  of  the  sampling  frequency,  fc 
Is  the  3  dB  cutoff  frequency  for  the  digital  filter. 

Assuming  that  the  sampling  frequency  Is  fixed  at  500  Hz,  l.e., 
f0  of  250  Hz,  the  only  variable  that  still  has  to  be  determined  Is  the  3 
dB  cut-off  frequency  fc.  Then  all  of  the  unknown  coefficients  In 
Equation  (A. 3)  can  be  obtained  using  Equations  (A-4.1)  to  (A-4.13).  The 
Implementation  for  the  transfer  function  of  the  fifth  order  digital 
filter  Is  shown  In  Figure  A-l.  The  implementation  of  H*(z)  Is  to  Input 
the  data  backward  through  the  filter  H(z)  and  reverse  the  output 
sequence. 

In  the  application  of  the  10th  order  zero-phase-shift  digital 
filter  to  the  noisy  data,  the  3  dB  bandwidth  of  the  filter  Is  chosen  to 
cover  the  target  size  and  transient  response.  The  transient  response 
waveforms  are  aspect  and  target  dependent.  However,  a  plot  of  the 
Impulse  response  may  help  to  decode  the  effective  transient  length. 
Normally,  4  to  5  TAU  (target  length  transit  time)  Is  a  good  choice  of 
the  3  dB  bandwidth  for  most  targets. 

An  example  of  the  determination  of  the  3  dB  cut-off  frequency  fc 
for  the  calculated  data  for  the  sphere  (Aka  of  0.02)  Is  shown  here. 

Since  the  data  for  the  sphere  are  sampled  In  the  frequency  domain,  the 
corresponding  Impulse  response  Is  periodic  with  a  period  tp  given  by 


where  Af  Is  the  frequency  Increment  between  successive  samples,  TAU  is 
2a/c,  a  is  the  radius  of  the  sphere,  and  c  is  the  speed  of  light  in  a 
free  space.  The  corresponding  3  dB  cutoff  frequency  fc  (in  terms  of 
target  length  transit  time)  for  the  digital  filter  is 

f  «  157.25  BW  (TAU).  (A-6) 

c 

where  BW  is  the  desired  3  dB  beamwidth  for  the  digital  filter  in  Hz. 


APPENDIX  B 


A  SUM  OPERATOR 

As  mentioned  In  Chapter  III,  an  overdeterml ned  system  of  linear 
equations  Is  obtained  after  some  manipulation  of  the  rational  function 
approximation  (Equation  (3.2)).  Assuming  the  RFA  has  a  system  order  of 
for  (M,N)  as  (3,2),  a  system  of  linear  equations  can  be  written  as 

AX  =  B  ,  (B-l) 

where  A  is  an  18x6  data  matrix  (assuming  that  9  complex  data  points  are 
used)  and  X  and  B  are  6  row  column  vectors,  l.e.,  6  unknown 
coefficients.  Thus 
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where  xj  Is  jkja,  are  the  real  part  of  the  known  Input  data  and  3i 
aro  the  Imaginary  part  of  the  Input  data.  Given  the  noisy  data  at  the 
ith  frequency 

E1  =  (ot1R  +  W  +  j(01I  +  *1IJ  .  <8*5 

where  jIjr  and  tjj  are  two  Gausslanly  distributed  pseudo  random 
noise  processes  with  zero  mean  and  standard  deviation  o.  ojr  and  Bfi 
are  the  real  part  and  Imaginary  part  of  the  noiseless  data. 

The  matrix  equation  for  additive  noise  Is 


“°1  -  *1R 

-JSl  “  J«.1 1 

■N. 

“°2  *  ^1 R 

-j&2  "  j*1I 
-«3  “  *1R 
—J  03  *  J*1I 


(B-9) 


-ag  -  HR  | 

-jBg  -  Hi  j 

"18  x  1 


In  order  to  reduce  the  additive  noise  In  each  element  of  the  matrix  A' 
and  B',  a  sum  operation  Is  used.  Each  element  of  matrix  A* 
and  matrix  B'  In  the  same  column  of  rows  1,  3  and  5*  are  added  and  then 
divided  by  3.  Thus,  a  new  row  Is  obtained.  The  same  procedure  Is 
applied  to  rows  7,  9,  and  11,  rows  13,  15  and  17,  rows  2,  4, and  6, 

.  .  .  and  rows  14,  16  and  18.  The  following  new  matrix  equation  results 
from  the  above  procedure 

A"X  *  B"  ,  (B-10) 


*  Choosing  the  odd  rows  to  be  a  group  and  the  even  rows  to  be  a  group, 
then  the  summation  will  be  limited  to  the  real  part  domain  and 
Imaginary  part  domain  respectively. 
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where  A”  Is  a  6x6  data  matrix  and  X  and  B"  are  6x1  column  vectors  In 
this  example.  The  matrix  equation  In  Equation  (B-10)  Is  exactly 
determined.  Thus,  this  matrix  equation  can  be  solved  In  the  exact 
sense.  The  variance  of  the  noise  Is  reduced  by  a  factor  of  about  3  In 
this  example.  Similarly,  the  sum  operation  can  be  applied  to  a  higher 
order  system  that  utilizes  more  data  points.  The  basic  reason  for  the 
noise  reduction  Is  an  averaging  as  shown  In  Equations  (3-18)  and 
(3-19). 

It  should  be  noted  that  there  Is  a  limitation  on  the  number  of 
samples  (N)  used  for  the  sum  operation,  because  the  noise  In  each 
element  of  the  matrix  A'  has  been  magnified  by  a  factor  of  x”  (see 
Equation  (B-7 ) ) ,  where  X^  *  jk^a  and  n  Is  an  Integer.  For  small  N  and 
small  frequency  Increments  (x^  -  *s  true  that 

Xi  *  X2  ■  X3  =»  .  .  .  ■  X||  (B— 11 ) 

so  that  the  reduction  factor  Is  N.  However,  for  larger  N,  the 
assumption  Is  no  longer  even  approximately  true.  Therefore,  the  noise 
reduction  factor  Is  less  than  N.  Considering  the  sum  operation  In  the 
above  section  of  this  Appendix,  It  Is  possible  to  define  a  sum 
operator,  S,  In  general  for  a  matrix  equation.  Given  a  matrix  equation 

CX  «  D  ,  (B-12) 


where 


’  Cn  Ci2 
C21  C22 

Cml  Cm2 


C13  ... 
C23 

Cm3  • • • 


cin“j 
C2n  | 
I 

I 

Cmn  I 

I 


mxn 


(B-13) 


n  x  1 


and 


(B-14) 


m  x  1 


( B— 15 ) 


The  sum  operator  S  Is  defined  as  an  average  of  I  elements  In  the  same 
column  but  In  I  different  rows.  Basically,  there  Is  no  limitation 
In  the  selection  of  those  I  rows  needed  In  group  summation.  However,  In 
the  application  of  a  sum  operator  to  a  matrix  equation  that  results  from 
a  RFA,  the  I  group  of  rows  are  the  ones  nearest  to  the  new  desired  row. 
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Moreover,  the  S  operator  can  be  applied  repeatedly.  Thus,  In  general , 
the  new  matrix  equation  using  one  S  operator  Is 

C'X  =  D*  ,  (B-l 

where 


r  sen 

SC12 

SC13  ... 

SCin  | 

|  sc2i 

1  • 

1  • 

SC22 

SC23  •  •  • 

s^2n  J 

s  • 

I  scmi 

scm2 

SCm3  •  •  • 

SCmn  1 

mxn 

and 


Note  again  here  that  the  sum  operator  Is  manipulated  In  the  matrix 
equation;  and  the  sum  operation  In  the  matrix  equation  Is  basically 
similar  to  row  echelon  algorithms  used  In  the  Gauss  elimination,  l.e., 


any  row  can  be  multiplied  by  a  constant  and  can  be  added  to  other  rows 


The  sum  operator  Itself  does  not  filter  out  any  of  the  signal*;  however, 
the  noise  Is  Indeed  reduced  by  using  the  sum  operator. 


*  When  the  sum  operator  and  the  RFA  are  applied  together  to  a  set  of 
noiseless  data,  the  result  Is  exactly  the  same  as  that  using  only  the 
RFA  to  the  noiseless  data. 
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APPENDIX  C 

A  COMPARISON  OF  THE  RAYLEIGH-STEVENSON  EXPANSION 

AND  THE  EXPANSION  FROM  RFA  FOR  A  DISC  AT  BROADSIDE 

In  this  Appendix,  a  list  of  the  unknown  coefficients  for  the 
rational  function  approximation  with  system  order  for  (M,N)  as  (14,13) 
to  the  backscattered  data  for  the  disc  at  broadside  excitation  Is  shown 
In  Table  C-l.  The  correspondl ng  poles,  zeros  and  residues  are  listed  In 
Table  C-2.  A  synthetic  division  of  the  rational  function  approximation 
(Equation  C-l))  using  the  coefficients  listed  In  Table  C-l  agrees 
closely  with  the  low  ferequency  power  series  expansion  (Equation  (C-2)) 
for  the  disc  at  broadside  given  by  Boersma  [47]  In  powers  of  X(jka)  up 
to  the  10th  order.  The  rational  function  yields 

E  (X)  =  .8486751X2  +  0.000437X3  -  0.4542859X4  +  .2433949X5 
n 

-  .1253862X6  -  .2085668X7  +  .0445261X8  +  .0943268X9 

-  0.0782452X10  +  .  .  .  .  (C-l) 

Compared  to  Boersma* s  result  of 
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I 


E  (X)  -  .84882636X2  +  OX3  -  .45270739X4  +  .24016873X5 
n 

+  .12934497X6  -  .21134848X7  +  .04495919X8 
+  .09588391X9  -  0.07943693X10  +  .  .  .  .  (C-2) 

Thus,  It  seems  that  the  RFA  Is  a  fairly  accurate  model  for  the  spectrum 
data  In  the  low  frequency  range. 
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TABLE  C-l 


A  LIST  OF  THE  COEFFICIENTS  OF  RFA  FITTED  TO  THE 
CALCULATED  BACKSCATTERING  DATA  FOR  THE  DISC  AT 
BROADSIDE  FOR  (M,N)  as  (14,13)  (lea  of  0.2(0.2)15.2) 


bl 

b2 

b3 

b4 

b5 

b6 

b7 

b8 

b9 

blO 

bll 

bl2 

bl3 

bl4 

•0 

•X 

a2 

a3 

«4 

«5 

a6 

•7 

a8 

a9 

alO 

all 

al2 

al3 


1.320147094030905 
1.357353646716751 
0.7226181666170916 
0.2757995040340761 
6.38468448943 84773B-02 
1 .39351666015423  998-02 
1 .7769630024111 417B-03 
2.53784 99633009358B-04 
1 .94593 541 02258490 E-05 
1.91575823 42 4014378-06 
8.67693132108662718-08 
5.89409653590826868-09 
1.29048252327416858-10 
5.66062481501240188-12 
0.8486743429354422 
1.120820652271342 
0.6982300296585792 
0.2575860119100928 
6.43968169652300958-02 
1.34124768927233328-02 
1.807 90771418965898-03 
2.47308357698554298-04 
1.98347667762588648-05 
1.87915772504621698-06 
8.84599433155288128-08 
5.80549765419457698-09 
1 .31554607644330038-10 
5.59395688772242778-12 


0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

0.00000008*00 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 

180.0000 
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TABLE  C-2 


A  LIST  OF  POLES,  ZEROS,  AND  RESIOUES*  WHICH 
CORRESPOND  TO  THE  RFA  COEFFICIENTS  IN  TABLE  C-l. 


(-4.4236 85331170900 ,-19.2850137051249) 
(-4.42360770X8895,19.2001341472798) 
(-1.269447207450867,-14.35744380950928) 
(-1.269446492195129,14.35743713378906) 
(-1 .242804050445557 ,-11 U 99048995 9716 ■> 
(•4.242005242538452,11.18808131373779) 
(-1.16 30043 943 45 947,-7. 981 972217559814) 
(-1.163003087043762,7.981966972351074) 
(-1 .033528023988342 ,-4.66654777526055) 
(-1 .033529281616211 ,4.66654777536055) 
(-1.551777601242065,-1.400591535564237) 
(-1.551777720451355,1.480591654777527) 


(-2.148778915405273,-7.35757056961972498-20) 


yprit 

(-4.550921680966912,19.51004429733050) 
(-4.550924052305061 ,-19.5300451007  1551 ) 
(-1.270Z3 33 92715454,14.41 97 4067647 988) 
(-1.270227313041687,-14.41974353790283) 
(-1.242350578308105,11.28768062591553) 
(-1.242354512214661,-11.28767204284668} 
(-1.164548754692078,8.105841636657715) 
(-1.16 45493 5073B25,-8.1 0584 83123 77 930) 
(-1 .0443 925857543 B ,4 .872451305389404) 
(-1.044391870498657,-4.872449874877930) 
(-1.635569214820062,-1.699864625930786) 
(-1.635569095611572,1 .699864149093628) 
(•0.4107 41941 08309,1 .230299234310 09) 
(-0.4907 4900150 29907 ^-1.23 02991151 80969) 


(-1.357 4670471 03 023 28-02,-2. 60 27321 8317303 36B-03) 
(-1.35747390760517118-02,2.60282512999685318-03) 
(-4. 241128953 0884 85 78-03,5. 84 836 49520162132B-04) 
(-4.24125260966873468-03,-5.85270482401078898-04) 
(-7.59766206923668288-03,1.56365080405460248  03) 
(-7.59700B3 584718568-03  ,-1.56315936475315848-03) 
(-1 .45227988528158568-02,3.76341249073819158-03) 
(-1.45236056 821 205838  02,-3.763565952965306 28-03) 
(-3 .8069529450704954802 ,1 .4498705899993531802) 
(-3.9068193050171831802,-1.44987181550969628-02) 
(-8.4212858934880469803,-8.2351453037297361802) 
(-8.4212342761281710803 ,8.2351287952638226802) 
(•0.407601655181451,0.3953218341053265) 
(-0.40760240 9S37QQ3  ,*0.3  ^216537022632) 


*  The  units  are  In  ka,  where  k  Is  the  wavenumber  and  a  Is  the  disc 
radius. 


APPENDIX  D 


THE  COMPLEX  NATURAL  RESONANCES  OF  THE  DISC 


Keller's  scalar  diffraction  coefficients  [11]  for  the  edge  at 
normal  Incidence  are  given  by 

-j  »/4 

D  *  (*,*')  «  -e  . .  Tsec  (JH[)  ?  sec  (J*JL)]  ,  (D-l-1) 

J  2v2nk  «  « 

Where  Ds  Is  edge  diffraction  coefficient  associated  with  the  soft 
(Dlrlchlet)  boundary  conditions,  D^,  is  the  edge  diffraction  associated 
with  the  hard  (Neumann)  boundary  condition,  if)'  is  the  incident  wave 
angle  and  *  Is  the  diffracted  wave  angle  as  shown  In  Figure  D-l. 

It  should  be  noted  that  Keller's  formula  Is  good  If  the  field  point 
is  not  close  to  a  shadow  or  reflection  boundary. 

As  the  field  point  Is  close  to  the  shadow  or  reflection  boundary 
(transition  regions),  the  Uniform  Theory  of  Diffraction  (UTD)  derived  by 
Kouyoumjlan  and  Pathak  [11]  should  be  used  instead  of  Keller's  formula. 
The  UTD  edge  diffraction  coefficients  for  curved  edge  or  straight  edge 
In  a  plane  screen  at  normal  Incidence  are 


00 


Figure  0-1.  The  coordinate  system  for  the  Incident  wave  diffracted 
by  a  half  plane. 


i  K 


1 1 

■j  *v 

!  _ 

i  ■“ 

I  r. 

i  I 

I  R 


-jk/4 

D  *  e _  {  F[2kLcos2U-V)/2]  .  F[2kLcos2(\|HV  )/2]  }  (d-1.2) 

?  z/m  "  cos  (*  -  ri  +  cos  cEPD  ’ 
h  2  2 

where 

F(x)  =  2j  /x  e"jx  /  e"jx2  dx  ,  (0-1.3) 

✓x 

A  Fresnel  integral,  and 

i  s  s(pe  +  s)  p}  p|  (D-1.4) 

pJ(p!  +  «)(pJ+  «)  ’ 

a  distance  parameter,  where  p|,  p£  are  the  principal  radii  of  curvature 
of  the  incident  wavefront  at  the  point  of  reflection,  p^  is  the  radius 
of  curvature  of  the  incident  waveform  at  the  diffraction  point  taken  in 

the  plane  containing  the  incident  ray  and  the  edge,  and  s  Is  the 
distance  from  the  diffraction  point  to  the  observation  point. 

The  equivalent  edge  currents  concept  was  proposed  in  [29]  to  solve 
the  field  in  the  caustic  direction  using  the  available  diffraction 
coefficients.  The  equivalent  electric  and  magnetic  edge  current  Ii  and 
Ml,  respectively  are  positioned  on  the  circular  edge  of  the  disc.  They 
are  given  by  [29]. 

_ _  A 

IlU')  »  -  Y0  M  DS(<M’)  (£*•♦')  ,  (0-2) 

MiU1)  =  -  Z0  ,  (0-3) 

3F 
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L 


where  ♦Is  the  unit  tangent  vector  to  the  circular  edge  of  the  disc, 
and  ♦  are  Interpreted  as  the  Incident  and  diffracted  ray  angles  as  shown 
In  Figure  0-2  and  0-3.  It  follows  that 
JjsI  -jkR  2it  - 

H  ( 0=0)  =  Jka  y  e  /  (xsln*'  -  ycos^'lM.d^'  ,  (0-4) 

4  it  o  R  0  1 


where 


M  =  H^lo.  (♦' 


[Sit  .  */4 

Vi r 


(D— 5 ) 


Note  In  Equation  (D.l-b)  should  be  used  for  Equation  (0.5)  due  to  the 
field  point  at  the  reflection  boundary.  Similar  first  order  scattering 
electric  field  using  the  equivalent  electric  current  concept  can  be 
obtained.  Furthermore,  the  doubly  and  triply  diffracted  field  (only 
exist  for  hard  case  In  this  example)  can  be  formulated. 


Hbs2(e*o)  a  jka  y 
T7  o 


2  it  A 

j  (x  sin*'  -  y  cos*')  M, 
0  £ 


d<|> ' 


where 


M 


2 


0, 


it)  ISlT  g"j  11 / 4 

7  V"TT 


f 


(0-6) 


(D-7) 


H1 

H2 


♦  'H  •  s1n$* 


['  ♦  8  0) 


/  p 

TaTFZaT 


-2jka 


(0-8) 


and 


p  Is  -  a  . 
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Figure  0-3.  The  singly,  doubly  and  triply  diffracted  mechanisms  of  the 
the  back scattering  from  the  circular  metallic  disc. 
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Ke+  -jkR  2ir  * 

Hdsx( e»0)  ■  Y  e  /  (x  sin*'  -  y  cos*')  M  d* 


where 


Mt  =  M2  M3  +  M4  +  •  •  • 

=  M2  (1  +  Y  +  Y^  +  •  •  •  ) 


*  M2/1-y 


(D-16) 


Y  ■  D.  (*'  -  0,  *  -  0)  • 


p  •  e”^ka| 


^a'Cp+TaT 


( p=-a) 


,-e'3,/4  .  J 

y^k  tit 

.  .  eJ  T/4  .  e"2jka 
2/jflca 


(D— 17 ) 


Consequently, 


i  k/4  -2ika 

1  -  Y  -  1  +  * _ e  3 

2/iflca 

j  n/2  2» 

-  1  +  e _ e 

2/iF 

*  1  +  j  exp{-2z)/(4irz)*^2  , 


(D-18) 


where  z  31  jka. 


*  Keller's  formula. 
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Equation  (0-18)  Is  the  charasterlstlc  equation  for  the  edge 
diffraction  node  for  the  disc  at  broadside  excitation.  The  26  pairs  of 
poles  have  been  extracted  using  a  zero  searching  routine  for  this 
equation.  They  are  listed  In  Table  D-l.  Note  the  first  3  pole-pairs 
listed  are  fro*  the  RFA  to  the  backscatterlng  data  for  the  disc  for 
broadside  excitation  because  the  GTD  Is  not  accurate  In  the 
low-frequencies.  Also  the  Induced  current  on  the  rim  due  to  the  singly 
diffracted  field  (see  Figure  D-3)  corresponds  to  the  lowest  frequency 
pole-pairs  of  the  creeping  wave  mode  because  the  variation  of  the  current 
along  the  rim  Is  exactly  one  period.  For  the  off-axis  Illumination,  the 
higher  order  diffracted  field  Is  due  to  the  multiple  diffractions  at  Qi 
and  Q2  (see  Figure  0-4).  Several  efforts  have  been  made  to  find  a  closed 
form  for  the  resonances  similar  to  the  one  derived  for  the  broadside 
Incidence  but  without  success.  Once  again  Keller's  edge  diffraction 
coefficients  are  used  here.  And  the  following  formula  Is  Introduced  for 
the  calculation  of  the  caustic  distance. 

A  A  A 

I  »  1  -  "e  *  -  s*  ,  (0-19) 

p  p["  *e 

where  p*  Is  the  radius  of  curvature  of  the  Incident  waveform  at  Qc  taken 
e  «  t 

In  the  plane  containing  the  Incident  ray  and  the  edge;  ne  Is  the 
associated  unit  normal  vector  to  the  edge  directed  away  from  the  center 

A 

of  the  curvature,  the  unit  vector  s'  Is  in  the  direction  of  Incidence  at 


L 
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TABLE  D-l 


A  LIST  OF  26  POLE*-PAIRS  FOR  THE  A  DISC  FOR  BROADSIDE 
EXCITATION  (BACKSCATTER) 


e  Number 

Real  part 

(+/-)Imai 

1 

-0.49069 

1.23103 

2 

-1.04579 

4.87347 

3 

-1.16349 

8.10608 

4 

-1.24179 

11.36106 

5 

-1.30243 

14.50748 

6 

-1.35120 

17.65236 

7 

-1.39201 

20.79634 

8 

-1.42709 

23.93976 

9 

-1.45784 

27.08279 

10 

-1.48523 

30.22556 

11 

-1.50991 

33.36812 

12 

-1.53238 

36.51053 

13 

-1.55299 

39.65282 

14 

-1.57203 

42.79502 

15 

-1.58972 

45.93715 

16 

-1.60625 

49.07920 

17 

-1.62175 

52.22122 

18 

-1.63634 

55.36318 

19 

-1.65013 

58.50512 

20 

-1.66320 

61.64701 

21 

-1.67562 

64.78889 

22 

-1.68746 

67.93073 

23 

-1.69875 

71.07256 

24 

-1.70956 

74.21437 

25 

-1.71992 

77.35616 

26 

-1.72987 

80.49794 

*  The  unit  of  poles  Is  ka,  where  k  Is  the  propagation  constant 
a  Is  the  radius  of  the  disc. 


252 


IV 


the  edge  and  s  Is  the  direction  of  diffraction,  ae  Is  the  radius  of 
curvature  of  the  edge  at  Qe. 

The  following  caustic  distances  are  obtained  using  Equation  (D-19) 
the  multiple  diffraction  at  Qj  and  Q2  shown  In  Figure  D-4. 


po-2  “ 
p0-2-l 


-a 

T+sTnF 

a  -a[l+2s1ne] 
■"TF3sTnS  “ 


p0-2-l-2 


-a[l+4s1ne] 

1+sslne 


3 


p 


p  0-2-l-2-la 


-a[l+6s1ne] 

l+7s1n0 


(D-20) 
(D— 21 ) 
(D-22) 
(0-23) 


p0-2-l-2-l-2 


-a[l+2Ns1n6] 

TT(ZRTlTsTnfl- 


(D-24) 


Although  the  diffraction  coefficients  are  the  same  In  each  order, 
the  caustic  distance  p  associated  with  each  diffraction  Is  not  the  same. 

Consequently,  the  factor  0.  •  p  ,  Is  different  In  each  order  of 

”  2a(p+2aT 

diffraction.  Therefore,  there  appears  to  be  no  way  to  form  a 
geometrical  series  similar  to  that  for  the  broadside  case. 
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APPENDIX  E 


THE  POLES  FOR  A  CIRCULAR  METALLIC  DISC  AT  EDGE-ON 

INCIDENCE 


The  analytic  form  for  the  creeping  wave  around  the  circumference 
of  a  disc  at  edge-on  Incidence  Is  given  by  Senior  [34]  as  follows: 

SCW  2  (1  -  I)"1/2  e"10*068ir(ka)2/3  e"«J2ka+j3ir/4  (E-l) 

TT 

•  exp{-1irka(l  -  1  (ka)“2/3  e"jn/3)}  , 

where  Scw  Is  the  scattering  amplitude  defined  In  [34],  y  Is  a  trial 
value,  for  example,  y  may  equal  0.8*PI.  In  terms  of  sCH,  the  normalized 
scattering  amplitude  Is  given  as 


eCW!  „  2  $cw/ka  # 


(E-2) 


or 


ECW1  ,  2  fi  .  y]-l/2  e' 10 * 068,1  (ka) 2/3  e'J2ka+J3lf/4 
n  lea  1  v 

•  exp{-1*ka(l  -  |  <ka )“2/3  e“J,r/3)}  , 


(E-3) 


or  substituting  u  for  Jka 


(E-4) 


EEW1  *  2(1-I)'1/2  u"1/3e“2y  ej5*/6  exp[-*y  +  »  y1/3]  . 

If  the  phase  reference  point  is  shifted  from  the  front  end  of  the  disc 
to  the  center  of  the  disc 

E™  *  2(1  -1)”1/2  u"1/3  e**-  exp[-iry  +  »  y1/3]  .  <E-5) 

If  we  also  consider  the  creeping  wave  due  to  an  additional  time  around 
the  disc,  then, 

EnW  "  2U-ir1/2  v'1/Ze*~  exp [ -itu  +  »  y1/3] 

•  exp[-2iry  +  rryl/3]  .  (E-6) 

Obviously,  successive  contributions  need  an  additional  phase  form  of 
P{u)=exp[-2iry  +  wpl/33  .  (E-7) 

The  complex  natural  resonances  are  obtained  by  Imposing  the  phase 
coherence  at  the  starting  points  after  one  or  more  revolutions. 

In  order  to  obtain  a  pole  string  for  the  creeping  wave  at  edge  on 
Incidence,  set 

-2wy  +  iry1/3  »  ±j2n*  •  (E-8) 

then 

w1/3  -  2y  *  ±j2n  .  (E-9) 

Let 
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n 


z  a  ka  , 


(E-10) 


z  a  Y3 


(E-U) 


Then, 


jl/3Zl/3  -  2jZ  a  ±j2n 


j1/3Y  -  2jY3  =  ±j2n 


(E-12) 


(E-13) 


Since 


jl/3  C  Ji  ±  2mrl/v 
J  =»  e  7 


(E— 14) 


The  principle  values  of  (j)1/3  are  obtained  from  (E-14)  at  na0,l,2 

eJS"/6,  eJ9’/6 

1*/6  <5ir/6 

or  eJ  ,  eJ  t  -j 

Any  of  the  three  principle  values  can  satisfy  (E-14).  We  choose 


(j)T--J  • 


(E-15) 


Then, 


-JY  -  2jY3  «  ±J2n 


(E-16) 


C 


The  zeros  of  Equation  (E-17)  are  related  to  the  poles  of  the  disc 
due  to  the  creeping  wave  mode  by  Equation  (E-8).  The  zeros  of  Equation 
(E-17)  can  be  found  by  using  the  DPOLYRT  Routine  (a  zero  searching 
routine  for  a  polynomial)  for  any  Integer  n.  A  pair  of  complex 
conjugate  zeros  and  one  real  zero  are  obtained.  However,  only  the 
complex  conjugate  poles  are  meaningful  when  they  are  substituted  into 
Equation  (E-10)  and  (E-ll).  Pure  Imaginary  resonances  could  Imply 
Internal  resonances.  A  list  of  the  first  thirty  pole-pairs  using 
Equation  (E-17)  at  n  *  1,  2,  ...  30  are  listed  In  Table  5-6. 


APPENDIX  F 


A  COMPARISON  OF  EXACT  AND  APPROXIMATE  BACKSCATTERED  FIELDS 
FOR  THE  CIRCULAR  DISK  FOR  EDGE-ON  INCIDENCE 

The  back scattered  field  for  the  disc  at  edge-on  Incidence  using 
exact  calculated  data  from  Hodge  [32]  and  Senior's  Geometrical  theory 
of  diffraction  (GTD)  result  [34]  are  compared  In  amplitude  and  phase  In 
Figures  F-l  and  F-2  respectively.  The  GTD  result  Is  not  good  at  low 
frequencies  due  to  the  asymptotic  nature  of  the  GTD  approximation. 
However,  the  amplitude  plots  at  high  frequencies,  as  shown  In  Figure 
F-l,  agree  well  except  that  there  Is  a  phase  lag  In  the  GTD  result. 
Figure  F-2  Is  a  comparison  of  the  phase  plots  both  from  GTD  and  the 
exact  result.  There  Is  a  phase  difference  between  the  calculated  data 
and  Senior's  approximation.  The  comparison  made  here  suggests  that 
additional  Improvement  could  be  made  In  the  asymptotic  model. 


APPENDIX  G 


A  REPRODUCTION  OF  ANDREJEWSKI 'S  NEAR  FIELD  DATA  OF  THE  DISC 

The  normalized  magnetic  scattered  near  field  for  the  disc  is 
reproduced  using  Andrejewski 's  formula.  From  Andrejewski,  the 
normalized  magnetic  near  field  (H)  Is  calculated  as 

H  8  1  t  1  W0!°»Y:°)*?,a,y)  -  Wil.Y.O)  ,  (G-l) 

yT  1'ot0,Y)-H|'2lO»Y)  0 

where  y  is  ka,  and  the  parameters  are  listed  in  the  following  Table 
(Table  G-l).  Using  Equation  (G-l)  and  Table  G-l,  14  data  points  (ka  of 
0.5(0. 5)4.  and  ka  of  4.(1.), 10.)  are  reproduced.  The  data  points  are 
then  approximated  by  a  rational  function  in  Chapter  V. 


TABLE  G-l 

*0,  <J>2  LISTED  IN  ANDREJEWSKI'S  DISSERTATION. 


411.  ,  w.«;  r;  w ,  X  (<;  <»,  <u. 


r 

-fi-w.it, r.o) 

K<4;r,0)  | 

ra 

1,024  ♦  1  0,5169 

0,105 

-  1  0,206 

j 

1,097  ♦  1  0,6230 

1 

0,4(9 

-  1  0,364 

1.5 

1,220  ♦  1  0,9057 

1.5 

1.25 

-  1  0,137 

2 

1,361  ♦  i  1,154 

2 

2,71 

♦  1  0,494 

*.» 

1,594  ♦  1  1,304 

2.5 

5,15 

♦  1  1,35 

5 

1,(91  *  1  1,425 

(.(5 

♦  1  2,00 

5.5 

1,73  ♦  1  1,53 

15,0 

♦  1  2.(5 

4 

1,963  ♦  1  1,644 

1(,1 

♦  1  1,09 

» 

2,0(7  ♦  1  1,(94 

20,4 

-  1  0,910 

6 

2,204  ♦  1  2,091 

37.5 

-  1  4,05 

7 

2,451  ♦  1  2.254 

45,0 

-  1  5,(0 

0 

2,555  ♦  i  2,491 

5(,5 

♦  1  0,S7( 

1 

2,759  ♦  1  2,592 

SiH 

70.3 

♦  i  5,81 

10 

2,(75  ♦  i  2,727 

104 

♦  »  (.19 

I 


I 


-0,0605  -  1  0,0000 
-0,2091  -  1  0,001# 

-0,5503  -  l  0,0206 
-0.5641  -  1  0.0259 
-U,79*1  -  1  0,0696 
-0,0169  -  1  0,2224 
-0,7652  -  1  0,99 JO 
-0,5790  -  1  0,5499 
-0,1998  -  1  0,6199 
0,4157  -  1  0,4692 
0,7701  ♦  1  0,0025 
0,4665  ♦  i  0,5742 
-0,2327  ♦  1  0,5095 
-0,7274  ♦  1  0,0955 
-0,576  -  1  0,4907 


rr 

'k  ci; jr) 

ESI 

0,1000  -  1  0,52(5 

1 

0,5971  -  1  0,(9(2 

1.5 

0,0631  -  1  1.162 

2 

1,445  “  1  1,750 

2.3 

2,075  -  1  2,465 

3 

2,72  -  t  3,251 

3,5 

5,40  -  1  3,(7 

4 

4,1(5  -  1  4,797 

5 

6,020  -  1  6,501 

6 

7,900  -  1  0,419 

7 

10,05  -  110,65 

0 
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APPENDIX  H 


SOME  TRANSIENT  RESPONSES 

In  this  Appendix,  several  transient  response  waveforms  are 
presented.  The  Impulse  response  and  the  ramp  response  of  the  loop  at 
edge-on  Incidence  are  presented  In  Figure  H-l  and  Figure  H-2.  The  ramp 
response  of  the  loop  for  broadside*  excitation  Is  shown  In  Figure  H-3. 
All  of  these  plots  are  obtained  using  Fourier  synthesis  of  the 
calculated  spectral  data  for  the  loop  at  edge  on  and  broadside 
excitation.  All  of  the  plots  shown  In  this  Appendix  are  limited  to  30 
transit  times  for  the  loop  diameter.  The  actual  waveform  durations 
are  somewhat  longer  than  those.  Singly  stated  that  circular  loop  (like 
the  wire)  Is  a  high  Q  scatterer. 

Figure  H-4  shows  ramp  response  waveforms  for  the  circular  disc  at 
several  aspect  angles  (cross  polarization)  where  10  harmonics  have  been 
used  In  the  synthesis  of  the  ramp  response. 


*  The  Impulse  response  of  the  loop  for  broadside  excitation 
has  also  been  shown  In  Chapter  V. 
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*  If  n  Is  N/2,  W(n)  Is  1/2. 
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APPENDIX  I 

K-PULSE  WAVEFORM  FOR  THE  THIN  CIRCULAR  LOOP 

The  K-pulse  waveform  for  an  object,  as  defined  by  Kennaugh  [6],  Is 
a  single  tine-limited  excitation  waveform  which  will  elicit  unique  tine 
United  response  waveforms  from  the  object  at  arbitrary  aspects.  One 
method  for  generating  the  K-pulse  Involves  a  weighted  product  of  the 
complex  natural  resonances  of  the  object  [44],  l.e.. 


E1(*)  -  n  [1  -  e 

*  j-1 


"e<,“pj)3  . 


d-l) 


where  e  Is  an  arbitrary  positive  number,  and  Pj  are  the  poles  of  the 
object.  Then,  El(s)  has  a  finite  duration  Inverse  Laplace 
transformation  of  length  (n  +  l)e.  Having  formulated  the  K-pulse  Input 
spectrum,  the  K-pulse  Input  waveform  Is  generated  by  the  Fourier 
synthesis  of  the  spectrum.  Note  the  weighting  function. 


.  Sin  2vn/N  .  (N/2)2 

25n7R -  (n/2)7^~ 

where  N  Is  the  number  of  sampled  data.  Is  used  to  avoid  Gibbs 


(1-2) 


phenomenon.  Since  the  K-pulse  and  response  waveforms  must  be  time 
limited,  the  K-pulse  actually  should  have  exactly  maximum  time  length  of 
(n  +  1 )e  as  specified,  these  waveforms  are  an  excellent  test  of  complex 
natural  resonances  locations  (poles)  deduced  using  approximate  methods. 
That  Is,  the  generated  K-pulse  and  response  waveforms  will  Indeed  be 
time-limited  If  the  pole  locations  are  accurate.  Some  allowance  must  be 
made  of  course  for  truncated  pole  products  used  to  obtain  the  K-pulse. 

The  most  simple  application  of  these  Ideas  Involves  a  single  string 
of  poles  (one  resonant- type  mode).  If  more  than  one  type  of  resonance 
Is  possible  the  K-pulse  becomes  more  complicated.  For  the  loop  at 
edge-on  Incidence,  resonance  Is  confined  to  the  creeping  wave  mode  and 
represents  an  excellent  test. 

The  first  10  pole-pairs  listed  In  Table  5-6  (creeping  wave  poles 
for  circular  loop)  were  used  to  obtain  the  K-pulse  spectrum  (amplitude 
and  phase)  and  waveform  for  the  loop  shown  In  Figure  1-1,  1-2  and  Figure 
1-3  respectively.  Figure  1-1  (amplitude)  and  Figure  1-2  (phase)  are  the 
K-pulse  Input  spectrum.  Figure  1-3  Is  the  K-pulse  Input  waveform.  Note 
the  K-pulse  length  Is  assumed  arbitrarily  to  be  *  TAU  (loop  diameter 
transit  time).  The  ringing  shown  Is  primarily  associated  with 
truncation,  and  not  Inaccuracies  of  the  poles.*  Response  to  the  K-pulse 
shown  In  Figure  1-3  Is  shown  In  Figure  1-4.  Clearly  there  Is  some 
slight  ringing  and  the  response  waveform  Is  not  precisely  time-limited. 

*  Unpublished  results  by  Kennaugh  confirm  this  postulate.  Using  20 
pole-pairs  deleted  most  of  the  ringing. 


Frequency  <kb) 


Figure  1-1.  The  K-pulse  Input  spectrum  (amplitude  plot)  generated  from 
Equation  (1-1)  using  the  first  10  pole-pairs  of  the 
circular  loop  at  edge  of  Incidence  (TE  mode).  The  pulse 
duration  Is  confined  to  »  TAU,  where  »  ■  3.14159  and  TAU 
Is  tiie  loop  diameter  transit  time.  Note  the  spectrum 
Is  spread  over  all  the  frequency,  however,  as  shown 
above,  the  amplitude  Is  snail  for  kb  greater  than  10. 
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Figure  1-2.  The  K-pulse  Input  spectrum  (phase  plot)  generated  from 

(1-1)  using  10  pole-pairs  for  the  circular  loop  (TE  mode). 
The  pulse  duration  at  edge-on  Is  confined  to  v  TAU,  where 
*  ■  3.14159  and  TAU  Is  the  loop  diameter  transit  time. 

The  data  are  truncated  at  ka  of  10  because  the 
corresponding  amplitude  Is  small  for  kb  greater  than  10. 
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.  A  K-pulse  Input  waveform  for  a  loop  for  edge-on  Incidence 
(TE  mode).  The  ringing  Is  due  to  the  truncation  of 
poles  (10  pole-pairs  used  here).  The  duration  Is 
confined  to  *  TAU,  where  *  Is  3.14159. 
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Figure  1*4.  A  K-pulse  output  waveform  for  a  loop  for  edge-on  Incidence 
The  duration  1$  longer  than  w  TAU,  where  w  Is  3.14159. 

This  Is  due  to  both  the  truncation  of  poles  and  spectrum. 


Figure  1-4  Is  generated  using  Fourier  synthesis  of  the  K-pulse  output 
spectrum,  which  Is  the  product  of  the  K-pulse  Input  spectrum  and  the 
transfer  function  of  the  target.  All  the  spectra  are  generated  In  the 
region  up  to  kb  of  10,  where  b  Is  the  loop  radius.  The  result  shown  In 
Figure  1-4  Is  similar  to  the  result  shown  by  Kennaugh  [6]  for  a  finite 
thln-wlre  scatterer.  Also  the  K-pulse  response  Is  about  *  TAU  lengths 
as  expected.  Therefore,  the  pole  locations  found  from  RFA  or  Integral 
equation  formulation  and  numerical  search  are  reasonably  accurate.  More 
Importantly,  It  would  appear  that  It  Is  feasible  to  obtain  K-pulse 
waveforms  from  pole  locations  found  from  measured  data. 

At  this  stage  we  are  admittedly  not  In  a  position  to  discuss  the 
K-pulse  and  response  waveforms  for  the  loop  (edge-on  Incidence)  In  the 
same  fashion  as  could  be  done,  for  example,  for  the  transmission  line 
by  Kennaugh  [6].  The  cancellation  Ideas  are  less  direct,  but 
still  occur.  A  more  detailed  study  of  the  loop.  Involving  general 
aspects  and  diffraction  as  well  as  creeping  wave  resonances,  would 
appear  to  be  a  fruitful  area  for  future  study. 


